- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
matlab指导书
实验一 数值运算实验
1 实验目的
(1)学习MATLAB语言的基本矩阵运算
(2)学习MATLAB语言的点运算
(3)学习复杂运算
2 实验内容
在下面的实验操作中,认真记录每项操作的作用和目的。
基本矩阵运算
创建数值矩阵
a=[1 2 3;4 5 6;7 8 9];
观察
a
a(3,2)
a(:,1)
键入
t=0:10;
u=0:0.1:10;
观察矩阵变量t、u的值。
键入
a(:,3)=[2;3;4];
a
观察矩阵a的变化。
键入
b=[1 1+2i;3+4i 3];
观察复数矩阵。
创建特殊矩阵
键入
a=ones(3,3);
b=zeros(2,2);
c=eye(4);
magic(4);
观察特殊矩阵。
练习矩阵运算
a=[0 1 0;0 0 1;-6 -11 -6];
b=[1 2;3 4;5 6];
c=[1 1 0;0 1 1];
作矩阵乘运算
v1=c*a
v2=a*b
v3=c*a*b
v4=b*c
v5=c*b
矩阵乘方运算
a^2
a^(1/2)
矩阵加减运算
a1=a+b*c
a2=c*b-a(1:2,1:2)
a3=a(1:2,2:3)+c*b
矩阵右除
ar=c/a
矩阵左除
al=a\b
练习矩阵特征运算
完成以下特征运算
inv(a),tril(a),rank(a)
MATLAB语言的点运算
练习点乘与点除
a1=[1 2;3 4];
a2=0.2*a1;
观察
[a1 a2]
[a1.*a2 a1./a2]
由点运算完成标量函数运算与作图
正余弦函数的点运算
t=0:2*pi/180:2*pi;
y1=sin(t);y2=cos(t);
y=y1.*y2;
plot(t,[y’ y1’ y2’]);
复变函数的点运算。
w=0.1:0.1:2;
g1=(1+0.5*w*i)/(1-0.5*w*i);
g1
g2=(1+0.5*w*i)./(1-0.5*w*i);
g2
plot(g2);xlabel(‘real g2(w)’);ylabel(‘inag g2(w)’)
axis(‘square’);
多项式运算
建立多项式向量
ap=[1 2 2 1];
b=[-1 -2 -3];bp=poly(b)
练习多项式乘与求根
p=conv(ap,bp)
roots(p)
练习多项式运算
a=[1 2 3 4];b=[1 -1];
c=a+[zeros(1,length(a)-length(b)) b];
poly2str(c,’x’)
polyvalm(a,3)
代数方程组
1)恰定方程组
A=[1 2;2 3];b=[8;13];
方法1:逆矩阵求解
X=inv(a)*b
方法2:矩阵左除求解
X=a\b
3 实验报告要求
按上述步骤完成数值运算实验,并按练习实例做好记录。
实验二 MATLAB语言绘图
1 实验目的
(1)学习MATLAB的各种二维绘图
(2)学习MATLAB的三维绘图
(3)MATLAB的绘图修饰
2 实验内容
(1)基本二维绘图
1)向量绘图
x=0:2*pi/100:2*pi;
y1=sin(2*x);y2=cos(2*x);
plot(x,y1)
plot(x,y2)
plot(x,y1,x,y2)
plot(x,y1);hold on;
plot(x,y2);hold off;
plot(x’,[y1’ y2’])
plot(x,y1,’c:’,x,y2,’ro’)
figure(1);plot(x,y1);
figure(2);plot(x,y2);
subplot(2,2,1);plot(x,y1);
subplot(2,2,2);plot(x,y2);
subplot(2,2,3);plot(x,y1,x,y1+y2);
subplot(2,2,4);plot(x,y2,x,y1-y2);
2)函数绘图
fplot(‘sin’,[0 4*pi])
f=’sin(x)’;fplot(f,[0 4*pi])
fplot(‘sin(1/x)’,[0.01 0.1],1e-3)
fplot(‘[tan(x),sin(x),cos(x)]’,[-2*pi,2*pi,-2*pi,2*pi])
3)符号函数快捷绘图
syms x
f=exp(-0.5*x)*sin(x);
ezplot(f,[0,10])
(2)多种二维绘图
1)半对数绘图
w=logspace(-1,1);
g=20*log10(1./(1+2*w*i));
p=angle(1./(1+2*w*i))*180/pi;
subplot(2,2,1);semilogx(w,g);grid;
subplot(2,1,2);semilogx(w,p);grid;
2)极坐标绘图
t=0:
文档评论(0)