Matlab实用教程分析.docVIP

Matlab实用教程 2.2 数组 2.2.1 建立数组 （1）逐个元素生成法 x=[1 2 3;2 3 4;2 4 5] （2）冒号法 x=1:1:5 （3）特殊方法 x=linspace(1,5,5) x=logspace(0,2,5) 在与之间插入5-2个数。 x=logspace(1,pi,5) 如果中间数为，在与之间插入5-2个数；如果中间数小于2，则返回值为。 2.2 引用数组 （1）一维数组 已知一维数组x=[1 2 3 4 5]，求x(1:3)。 （2）二维数组 二维数组x=[1 2 3;2 3 4;2 4 5]，求x(2,2:3)。 2.3 运算符 2.3.1 算术运算 定义 加 减 乘 右除 左除 乘方 操作符 + - * / \ ^ 点运算 .* ./ .\ .^ 2.3.2 关系运算 定义 等于 不等于 大于 大于等于 小于 小于等于 操作符 == ~= = = 2.3.3 逻辑运算符 定义 逻辑与 逻辑或 逻辑非 逻辑异或 操作符 | ~ xor Matlab语言中的逻辑运算规则 输入 与 或 异或 非 a b ab a|b xor(a,b) ~a 0 0 0 0 0 1 0 1 0 1 1 1 1 0 0 1 1 0 1 1 1 1 0 0 所有非0数均被认为是真，而0为假；在逻辑判断结果中，判断为真时输出1，判断为假时输出0。例： A=[1 1;0 1]; b=0; B=[0 1;0 0]; C1=Ab C2=A|b C3=xor(A,B) 2.4 m文件 例 分别绘制花瓣图案，，，。 clear; theta=-pi:0.01:pi; rho(1,:)=2*sin(5*theta).^2; rho(2,:)=cos(10*theta).^3; rho(3,:)=sin(theta).^2; rho(4,:)=5*cos(3.5*theta).^3; for i=1:4 subplot(2,2,i); polar(theta,rho(i,:)) end 【例】定义一个名为f.M的函数文件，并调用该函数。 function z=f0226(x,y) z=sin(x.^2)+y.^2-2*x-5; a=[1 3]; b=[2 2]; c=f0226(a,b) 2.5 程序设计 2.5.1 顺序结构 【例】顺序结构 disp(the begin of the program) disp(the first line) disp(the end of the progrogram) 2.5.2 分支结构 （1）条件语句 【例】 a=1; if a==0; disp(nihao1); else disp(nihao2); end 【例】 a=0; if a==0; disp(nihao1); else if a==1; disp(nihao2); else end end （2）开关语句 定义函数 function y=f0233(x) switch x case {1,2} y=x*0.1; case {3,4} y=x*0.4; otherwise y=x*0.5; end 在命令窗输入 a=2.5; b=f0233(a) 2.5.3 循环结构 【例】 for(i=1:3) for (j=1:4) A(i,j)=1.5; end end A 【例】 i=0; while(i==5) A=1.5; end; A 3.6 求解线性方程组 线性方程组的求解可以分为两类：一类是方程组求唯一解或求特解，另一类是方程组求无穷解或通解，可以通过方程组系数矩阵的秩来判断： 若系数矩阵的秩，（为方程组中未知量的个数），则有唯一解。 若系数矩阵的秩，则可能有无穷多解。 线性方程组的无穷解=对应齐次方程组的通解+非齐次方程组的一个特解。 3.6.1 线性方程组唯一解或特解的解法 （1）矩阵的除法 【例】 求方程组 A=[5 6 0 0;1 5 6 0;0 1 5 6;0 0 1 5]; b=[1 0 0 0]; R_A=rank(A) %求系数矩阵的秩 X=A\b 另一种解法 C=[A,b]; R=rref(C) 3.6.2 齐次线性方程组的通解的解法 Matlab中函数null用来求解零空间，即满足的空间，实际上求出解空间的一组基（基础解系），它可用于获得齐次线性方程组的通解。 【例】 求解方程组的通解。 A=[1 2 2 1;2 1 -2 -2;1 -1 -4 -3]; format rat