matlab实验汇总….docVIP

2014电科上机题 第一次上机：1、熟悉桌面平台及帮助系统 p5 2、向量、矩阵及运算 p30 3、用左除和右除求解线性方程组 p36,37 第二次上机：1、矩阵的特殊操作 p50 2、多项式运算 p59 3、多项式拟合 p61 4、设原始数据为x， 在11个点上测得的y值如下： X=[-2.0 -1.6 -1.2 -0.8 -0.4 0 0.4 0.8 1.2 1.6 2.0]; Y=[2.8 2.96 2.54 3.44 3.56 5.4 6.0 8.4 9.5 13.3 15]; 采用2阶和10阶多项式拟合，并画图比较。 a=polyfit(x,y,2); x1=-2.0:0.01:2.0; y1=a(1)*x1.^2+a(2)*x1+a(3); b=polyfit(x,y,10); x2=-2.0:0.01:2.0; y2=b(1)*x2.^10+b(2)*x2.^9+b(3)*x2.^8+b(4)*x2.^7+b(5)*x2.^6+b(6)*x2.^5+b(7)*x2.^4+b(8)*x2.^3++b(9)*x2.^2+b(10)*x2+b(11); plot(x,y,b*,x1,y1,r-,x2,y2,g-) legend(二阶,十阶) 5、设a1=[2 4 6 8]为系统分母系数向量，b1=[3 6 9]为系统分子系数向量，求此系统的频率响应并画出频率特性。令频率数组w取线性间隔(0-10)。 w=0:0.01:10; s=j*w; A=3*s.^2+6*s+9; B=2*s.^3+4*s.^2+6*s+8; H=A./B; subplot(1,2,1) plot(w,abs(H)) title(幅值--w) subplot(1,2,2) plot(w,angle(H)) Title(角度--w) 第三次上机：1、符号表达式的生成p64；符号矩阵的创立p68；符号矩阵的简化p73。 2、符号代数方程求解。P79 3、求解： 2cos(x)+y^(1/2)-ln(z)=7 2^x+2y-8z=-1 x+y-cosh(z)=0 先建立.m文件： function y=fc(x); y(1)=2*cos(x(1))+x(2)^(1/2)-log(x(3))-7; y(2)=2^x(1)+2*x(2)-8*x(3)+1; y(3)=x(1)+x(2)-cosh(x(3)); y=[y(1) y(2) y(3)]; 第四次上机：1、符号求解 d3u/dx^3=u, u(0)=1,u’(0)=-1,u’’(0)=pi 命令窗口执行： dsolve(D3u=u) ans = C1*exp(t)+C2*exp(-1/2*t)*sin(1/2*3^(1/2)*t)+C3*exp(-1/2*t)*cos(1/2*3^(1/2)*t) y=dsolve(D3u=u, u(0)=1,Du(0)=-1,D2u(0)=pi) y = 1/3*pi*exp(t)-1/3*(1+pi)*3^(1/2)*exp(-1/2*t)*sin(1/2*3^(1/2)*t)+(-1/3*pi+1)*exp(-1/2*t)*cos(1/2*3^(1/2)*t) simplify(y) ans = 1/3*pi*exp(t)-1/3*3^(1/2)*exp(-1/2*t)*sin(1/2*3^(1/2)*t)-1/3*3^(1/2)*exp(-1/2*t)*sin(1/2*3^(1/2)*t)*pi-1/3*exp(-1/2*t)*cos(1/2*3^(1/2)*t)*pi+exp(-1/2*t)*cos(1/2*3^(1/2)*t) 2、分别用solve和fsolve函数求方程组 　　　 x^2*y^2 - 2*x - 1 = 0 　　　 x^2 - y^2 - 1 = 0 Solve求解： [x,y]=solve(x^2*y^2-2*x-1=0,x^2-y^2-1=0) x = 1/2+1/2*5^(1/2) 1/2+1/2*5^(1/2) 1/2-