MATLAB基础上机练习题.docx

MATLAB基础练习上机例题讲解矩阵的表示与运算例1：设，，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。解： A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; C1=A*BC1 = 19 -82 30 12 27 3 -38 54 29 C2=A*BC2 = -15 16 -24 36 63 -1793 -105 22 6 117 -60 19 46 84 -10 C3=A.*BC3 = 5 16 24 -26 -18 -12 -15 72 -2 -21 108 -56 inv(C1)ans = 0.0062 0.0400 -0.0106 -0.0046 0.0169 0.0030 0.0168 0.0209 0.0150 inv(C2)Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 8.997019e-019.ans = 1.0e+015 * -0.9553 -0.2391 -0.1997 0.2700 0.9667 0.2420 0.2021 -0.2732 -0.4473 -0.1120 -0.0935 0.1264 -1.1259 -0.2818 -0.2353 0.3182 inv(C3)??? Error using == invMatrix must be square.方程组求解例2：求下列联立方程的解解： a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; b=[4;4;9;4]; c=a\bc = 5.2226 4.4570 1.4718 1.5994多项式的表示与运算例3： 已经两个多项式的系数分别是：[1 2 3 4]和[1 4 9 16]，请求这两个多项式的乘积，及商和余数。解：: p1=[1 2 3 4]; p2=[1 4 9 6]; c=conv(p1,p2)c = 1 6 20 40 55 54 24 [q,r]=deconv(p1,p2)q = 1r = 0 -2 -6 -2函数的表达与函数值计算例4：在matlab的命令窗口里完成如下计算，其中t的值分别取-1,0,1，表达式如下： 解： t=[-1 0 1]t = -1 0 1 y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3*t)+pi/3)y = 1.0e+004 * 2.0897 + 1.2065i 0.0001 0.0000零极点计算例5：已知：一个多项式的系数向量是p=[1 -6-72 -27]，求这个多项式的根。解： p=[1 -6 -72 -27]p = 1 -6 -72 -27 r=roots(p)r = 12.1229 -5.7345 -0.3884例6：用fminsearch函数求如下函数的极小值：解： f f= inline(‘100*(x(2)-x(1)^2)^2+(1-x(1))^2’,’x’); x0= [-1.2 1]; [sx,sfval] = fminsearch(ff,x0);sx = 1.00001.0000sfval = 8.1777e-010符号运算例7：用符号运算求解：解：例8：用符号运算求解：解：绘图基础例9：设 ，把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线。解： x=linspace(0,2*pi,101); y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); plot(x,y,r)数据拟合例10：已知一组测量值M语言编程解： t=[1 2 3 4 5 6 7 8 9 10]; y=[15.0 39.5 66.0 85.5 89.0 67.5 12.0 -86.4 -236.9 -448.4]; p1=polyfit(t,y,2)p1 = 0.0411 -0.6383 5.3937 p2=polyfit(t,y,3)p2 = -0.0041 0.1085 -0.9493 5.7443 y3=polyval(p1,x); y4=polyval(p2,x); plot(t,y,O,x,y3,r:,