Matlab基本介绍重点.ppt

什么是数学实验？;《数学实验》课程 ;先看两个例子; 例2 导弹打击过程仿真 zuijiwent; MATLAB程序： ? k=0; A=[0,0]; %导弹初始位置 B=[0,100]; %飞机初始位置 v=1; dt=1; %离散时间改变量 d=100; %相距距离 while d0.5 plot(A(1),A(2),’r+’); %画导弹位置 hold on plot(B(1),B(2),’b*’); %画飞机位置 pause(0.2); k=k+1; B=B+[v*dt,0]; %飞机移动位置 e=B-A; %导弹指向飞机向量 d=norm(e); e0=e/d; %取向量方向（单位化） A=A+2.0*v*dt*e0; %导弹追击位置 fprintf(k=%.0f A(%.2f,%.2f) B(%.2f,100) d=%.2f\n,k ,A(1),A(2),B(1),d); end ;课程内容 ;《数学实验》课程的考核;MATLAB命令简单介绍 ;运行MATLAB创建一个或多个窗口 ; 工具箱及命令查询 help topics;MATLAB中基本代数运算符;变量及数组输入 ;构造矩阵的方法------q3;求方程组Ax=B的解-----------q1.m;A=[1 2 0;2 5 -1;4 10 -1]; %矩阵输入（ A为3阶矩阵） B=[366;804;351] %列矩阵输入 D=A’ %求A的转置 H=det(A) %求A的行列式 r=rank(A) %计算矩阵A的秩 ni=inv(A) %求A的逆 x=A\B %Ax=B方程组的解 y=inv(A)*B %与x相同 tezhengzhi=eig(A) %求A的特征值 [m,n]=eig(A) %求A的特征值与特征向量 B=poly(A) %求A的特征多项式 round(poly(A)) %round用来对数值取整 rref(A) %计算矩阵A的行最简形;求向量组a1＝（l -2 2 3），a2＝（-2 4 -1 3），a3＝（-1 2 0 3）,a4＝(O 6 2 3)，a5＝（2 -6 3 4）的秩和一个极大线性无关组；并求由该向量组为列向量的矩阵的标准阶梯形。 解 a1=[1 -2 2 3]; a2=[-2 4 -1 3]; a3=[-1 2 0 3]; a4=[0 6 2 3]; a5=[2 -6 3 4]; A=[a1;a2;a3;a4;a5]; format rat %以有理格式输出 b=rank(A), B=rref(A); 例 求解方程组 的基础解系。 解 在Matlab编辑器中建立M文件:f.m A=[1,2,2,1;2,1,-2,-2;1,-1,-4,-3]; format rat %指定有理式格式输出。 b=rref(A) B=null(A,r) %求解空间的基。 syms k1 k2 %定义符号常数。 X=k1*B(:,1)+k2*B(:,2) %写出方程组的通解。;MATLAB中的数值函数运算 ;MATLAB中符号函数的表示; MATLAB中符号运算－－求极限 -------------q0.m; MATLAB符号计算----求导数 ---q4;MATLAB符号计算---求积分-----q5; 级数求和与函数展开－－q6; 代数方程求根----------q7;二维平面图形------q9;图形的比较显示----在同一窗口绘制多图方法1--------q10;图形的比较显示----在同一窗口绘制多图