matlab实验汇总.doc
文本预览下载声明
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^
显示全部