资源预览内容
第1页 / 共24页
第2页 / 共24页
第3页 / 共24页
第4页 / 共24页
第5页 / 共24页
第6页 / 共24页
第7页 / 共24页
第8页 / 共24页
第9页 / 共24页
第10页 / 共24页
亲,该文档总共24页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
第二章 数值数组及其运算,秦伟刚电气工程与自动化学院4D506,例子:,绘制函数 在 时的曲线。,x=0:0.1:1y=x.*exp(-x)plot(x,y),xlabel(x),ylabel(y),title(y=x*exp(-x),一维数组的子数组寻访,rand(state,0) %产生均匀分布的随机数x=rand(1,5) %生成1行5列的随机数x(3) x(1 2 5) x(1:3) x(3:end) x(3:-1:1) x(find(x0.5),一维数组的子数组赋值,x(3) = 0 x(1 4)=1 1,二维数组的创建,在MATLAB环境下,用下面三条指令创建二维数组C 。a=2.7358; b=33/79;C=1,2*a+i*b, b*sqrt(a);sin(pi/4),a+5*b,3.5+i,二维数组的创建,复数数组的另一种输入方式。 M_r=1,2,3;4,5,6,M_i=11,12,13;14,15,16CN=M_r+i*M_i,二维数组的子数组寻访和赋值,A=zeros(2,4) %全零数组%赋值A(:)=1:8%赋值s=2 3 5;A(s)Sa=10 20 30A(s)=Sa,执行数组运算的常用函数,执行数组运算的常用函数,执行数组运算的常用函数,执行数组运算的常用函数,【例】演示pow2的数组运算性质。A=1:4;5:8 A = 1 2 3 4 5 6 7 8 pow2(A) ans = 2 4 8 16 32 64 128 256,数组运算和矩阵运算的区别,两种不同转置的比较,例:clear;A=zeros(2,3);A(:)=1:6;A=A*(1+i)A_A=A.A_M=A,多项式表达方式的约定,MATLAB约定降幂多项式P(x)=a1xn+a2xn-1+Anx+an+1用以下系数行向量表示: P=a1 a2 an an+1,多项式行向量的创建方法,1.直接输入法 多项式系数应以降幂次序进行,如果多项式中缺某幂次项,那么认为该幂次项的系数为零。2.利用指令:P=poly(AR),产生多项式系数向量 若AR是方阵,则多项式P就是该方阵的特征多项式。 若AR=ar1 ar2 arn,则AR的元素被认为是多项式P的根。所得的多项式应满足 (x-ar1)(x-arn)=a0xn+a1xn-1+an-1X+an,多项式行向量的创建方法,【例 】求3阶方阵A的特征多项式。 A=11 12 13;14 15 16;17 18 19; PA=poly(A) %A的特征多项式 PPA=poly2str(PA,s) %显示多项式,多项式运算函数,多项式运算,【例】求 的“商”及“余”多项式。p1=conv(1,0,2,conv(1,4,1,1);p2=1 0 1 1;q,r=deconv(p1,p2);cq=商多项式为 ; cr=余多项式为 ;disp(cq,poly2str(q,s),disp(cr,poly2str(r,s),多项式运算,【例 】部分分式展开。 a=1,3,4,2,7,2;%分母多项式系数向量 b=3,2,5,4,6; %分子多项式系数向量 r,s,k=residue(b,a) 注意:当分母阶数高于分子时,计算结果中的k将是空阵。,标准数组生成函数,数组操作函数,“非数”和“空”数组,非数NaN 【例】非数的产生a=0/0,b=0*log(0),c=inf-inf “空”数组【例】 “空”数组的产生a=,b=ones(2,0),c=zeros(2,0),d=eye(2,0),关系操作,例,【例】关系运算示例。A=1:9,B=10-A,r0=(A0) %判断哪些元素不大于0 L2=A0 %逻辑非优先于关系运算 L3=AL4=A-2&A1 %判断哪些元素大于-2小于1,关 %系运算优先于逻辑与,
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号