Matlab作业电气工程专业试题及答案.docx
- 文档编号:8971377
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:18
- 大小:180.09KB
Matlab作业电气工程专业试题及答案.docx
《Matlab作业电气工程专业试题及答案.docx》由会员分享,可在线阅读,更多相关《Matlab作业电气工程专业试题及答案.docx(18页珍藏版)》请在冰点文库上搜索。
江苏大学
京 江 学 院
《电气工具软件训练Ⅲ》课程设计报告
设计题目:
MATLAB
专业班级:
学生姓名:
学生学号:
指导老师:
完成日期:
10/16
江苏大学·电气信息工程学院
—MATLAB课程设计的目的和要求
1.MATLAB软件功能简介
MATLAB的名称源自MatrixLaboratory,1984年由美国
Mathworks公司推向市场。
它是一种科学计算软件,专门以矩阵的形式处理数据。
MATLAB将高性能的数值计算和可视化集成在一起,并提供了大量的内置函数,从而被广泛的应用于科学计算、控制系统和信息处理等领域的分析、仿真和设计工作。
MATLAB软件包括五大通用功能,数值计算功能(Nemeric)、符号运算功能(Symbolic)、数据可视化功能(Graphic)、数字图形文字统一处理功能(Notebook)和建模仿真可视化功能(Simulink)。
其中,符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLAB命令窗口。
该软件有三大特点,一是功能强大;二是界面友善、语言自然;三是开放性强。
目前,Mathworks公司已推出30多个应用工具箱。
MATLAB在线性代数、矩阵分析、数值及优化、数值统计和随机信号分析、电路与系统、系统动力学、次那好和图像处理、控制理论分析和系统设计、过程控制、建模和仿真、通信系统以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。
2.MATLAB课程设计的目的
本次课程设计主要是为了使学生了解MATLAB软件的基本知
识,熟悉MATLAB的上机环境,掌握MATLAB数值运算、程序设计、二维/三维绘图、符号运算、Simulink仿真等相关知识,并初步具备将一般数学问题转化为对应的计算机进行处理的能力,以便为今后进一步的学习打下坚定基础。
二MATLAB课程内容
1MATLAB语言基础
实验目的:
基本掌握MATLAB向量、矩阵、数组的生成及其基本运算(区分数组运算和矩阵运算)、常用的数学函数。
了解字符串的操作。
实验内容:
①创建以下矩阵:
A为初值为1,终值为12,元素数目为6的行向量;
æ2
4
3
6ö
ç1
ç
5
3
2÷
ç2
7
7
9÷
ç1
2
3
5÷
ø
B= ÷;C为B的三次方矩阵;D由B和C横向拼接并去除第一列、
è
最后一列和第一行元素而成;E由B和C纵向拼接而成;F抽取E的3~5行和第2~3列元素生成;G由F经变形为3×2的矩阵而得;H由B和C数组相乘运算而成,同时将H(1,1)和H(2,1)分别变为π的平方和立方,H(2,2)=arccos
(2),H(3,3)=H(1,1)+H(2,1)。
源程序:
A=linspace(1,12,6)B=[2,4,3,6;1,5,3,2;2,7,7,9;1,2,3,5];C=B^3
D1=[B,C];
D=D1(2:
4,2:
7)E=[B;C]F=E(3:
5,2:
3)
G=reshape(F,3,2)H=B*C
H(1,1)=pi^2;
H(2,1)=pi^3;H(2,2)=acos
(2);H(3,3)=H(1,1)+H(2,1)
运行结果:
A
C
=1.0000
=
3.2000
5.4000
7.6000
9.8000
12.0000
292
938
873
1140
227
751
681
868
518
1675
1555
2021
227
724
681
895
D
=
5
3
2
227
751
681
7
7
9
518
1675
1555
2
3
5
227
724
681
E
=
2
4
3
6
1
5
3
2
2
7
7
9
1
2
3
5
292
938
873
1140
227
751
681
868
518
1675
1555
2021
227
724
681
895
F
=
7 7
2 3
938 873
G
=
7 7
2 3
H
938 873
=
4408 14249
13221
17185
3435 11166
10305
13333
7842 25374
23527
30558
H
3435 11085
=
1.0e+004*
10305
13414
0.0010 1.4249
1.3221
1.7185
0.0031 0
+
0.0001i
1.0305
1.3333
0.7842 2.5374
0.0041
3.0558
0.3435 1.1085
1.0305
1.3414
②
(1)用矩阵除法求下列方程组的解x=[x1;x2;x3];
ì6x1+3x2+4x3=3
ï-2x+5x+7x=-4
í 1
ï8x-x
2
-3x
3
=-7
î 1 2 3
(2)求矩阵的秩(rank函数);
(3)求矩阵的特征值与特征向量(eig函数);
(4)系数矩阵的3次幂与开方;
(5)系数矩阵的指数运算和数组对数运算;
(6)系数矩阵a(1,2)、a(1,3)、a(2,2)、a(2,3)的元素不变,其余元素变为零。
(7)提取系数矩阵主对角线上的元素,并依次相加赋予b。
源程序:
A=[6,3,4;-2,5,7;8,-1,-3];B=[3;-4;-7];
X=inv(A)*Brank(A)[v,d]=eig(A)A1=A^3
A2=sqrtm(A)A3=expm(A)A4=logm(A)A(:
1)=0;
A(3,:
)=0
A=[6,3,4;-2,5,7;8,-1,-3];
b=trace(A)运行结果:
X=
1.0200
-14.0000
9.7200
ans=3v=
0.8013
-0.1094
-0.1606
0.3638
-0.6564
0.8669
0.4749
0.7464
-0.4719
d=
9.7326
0
0
0
-3.2928
0
0
0
1.5602
578
298
352
228
156
202
384
154
156
2.2447
+
0.2706i
0.6974
-
0.1400i
0.9422
-
0.3494i
-0.5815
+
1.6244i
2.1005
-
0.8405i
1.7620
-
2.0970i
1.9719
-
1.8471i
-0.3017
+
0.9557i
0.0236
+
2.3845i
A3=
A1=
A2
1.0e+004
1.0653
*
0.5415
0.6323
0.4830
0.2465
0.2876
0.6316
0.3206
0.3745
1.7129
+
0.4686i
0.5305
-
0.2425i
0.5429
-
0.6049i
1.1938
+
2.8123i
0.3658
-
1.4552i
-0.5514
-
3.6305i
-0.0748
-
3.1978i
0.7419
+
1.6546i
1.8333
+
4.1282i
A=
0
3
4
0
5
7
0
0
0
b=
A4=
8
2MATLAB数值运算
实验目的:
掌握MATLAB的数值运算及其运算中所用到的函数,掌握结构数组的操作。
实验内容:
①已知多项式a(x)=x2+2x+3;b(x)=4x2+5x+6
(1)求多项式a(x)和多项式b(x)的乘法运算结果,并在命令窗口中显示该多项式c;
(2)求多项式c的根及其微分;
源程序:
p1=[1,2,3];
p2=[4,5,6];
y1=poly2sym(p1);y2=poly2sym(p2);p=conv(p1,p2);
c=poly2sym(p)roots(p)c1=diff(c)
运行结果:
c=4*x^4+13*x^3+28*x^2+27*x+18ans=
-1.0000+1.4142i
-1.0000-1.4142i
-0.6250+1.0533i
-0.6250-1.0533i
c1=16*x^3+39*x^2+56*x+27
(s2+1)(s+3)(s+1)
②求
s3+2s+1
的“商”及“余”多项式并在命令窗口中显示该多项式。
源程序:
root=[-i,i,-3,-1];p=poly(root);disp(poly2sym(p))d=[1,0,2,1];
c2=poly2sym(d)[q,r]=deconv(p,d);c3=poly2sym(q)c4=poly2sym(r)
运行结果:
c3=x+4
c4=2*x^2-5*x-1
③
(1)计算当x=2,x=3时,错误!
未找到引用源。
(x-0.98)2 1
=3+ - + 的
f(x) x
值;
(x+1.25)3
5(x )
x
9-2
(2)计算错误!
未找到引用源。
cos60+arccos(p)- 的值
æ2
4
3
6ö
ç1
ç
5
3
2÷
÷
ç2
7
7
9÷
ç1
2
3
5÷
ø
(3)A= ,B=A2+3,C=A-2B,,求:
C
è
源程序:
x=2
y1=x^3+(x-0.98)^2/(x+1.25)^3-5*(x+1/x)x=3
y1=x^3+(x-0.98)^2/(x+1.25)^3-5*(x+1/x)
y2=cos(60)+acos(pi)-sqrt(9-sqrt
(2))A=[2,4,3,6;1,5,3,2;2,7,7,9;1,2,3,5]B=A^2+3;
C=A-2*B
运行结果:
x=2
y1=-4.4697
x=3
y1=10.3865
y2=-3.7066+
1.8115i
A=
2 4
3
6
1 5
3
2
2 7
7
9
1 2
3
5
C=
-44
-124
-117
-154
-35
-109
-93
-110
-72
-219
-205
-265
-35
-94
-93
-125
3MATLAB符号运算
实验目的:
掌握符号变量和符号表达式的创建,掌握MATLAB的symbol工具箱的一些基本应用。
实验内容:
①已知
f=(ax2+bx+c-3)3-a(cx2+4bx-1),按照自变量x和自变量a,对表
达式f分别进行降幂排列(同幂合并)。
源程序:
symsabcx
f=sym('(a*x^2+b*x+c-3)^3-a*(c*x^2+4*b*x-1)')collect(f)
collect(f,a)
运行结果:
f=(a*x^2+b*x+c-3)^3-a*(c*x^2+4*b*x-1)
ans=a^3*x^6+3*b*a^2*x^5+((c-3)*a^2+2*b^2*a+a*(2*(c-3)*a+b^2))*x^4+(4*
(c-3)*b*a+b*(2*(c-3)*a+b^2))*x^3+((c-3)*(2*(c-3)*a+b^2)+2*b^2*(c-3)+a
*(c-3)^2-a*c)*x^2+(3*(c-3)^2*b-4*b*a)*x+(c-3)^3+aans=a^3*x^6+3*(b*x+c-3)*x^4*a^2+(3*(b*x+c-3)^2*x^2-c*x^2-4*b*x+1)*a+(b*x+c-3)^3
②已知f1=1/(a-b),f2=2a/(a+b),f3=(a+1)(b-1)(a-b),分别求f1和f2的符号和、
f1和f3的符号积、f1和f3的符号商。
源程序:
symsabf1=sym('1/(a-b)');
f2=sym('2*a/(a+b)');f3=sym('(a+1)*(b-1)*(a-b)');f1+f2
f1*f3f1/f3
运行结果:
ans=1/(a-b)+2*a/(a+b)ans=(a+1)*(b-1)
ans=1/(a-b)^2/(a+1)/(b-1)
③对下列表达式进行符号运算
(1)已知数学表达式y(x)=(ex+x)(x+2),将其展开。
(2)已知数学表达式y(x)=a3-1,对其进行因式分解。
(3)已知数学表达式y(x)=
x+3+
x(x+1)
x-1
x2(x+2)
,对其进行通分。
(4)已知数学表达式y(x)=2cos2x-sin2x,对其进行化简。
源程序:
symsxay1=sym('(e^x+x)*(x+2)');y2=sym('a^3-1');
y3=sym('(x+3)/(x*(x+1))+(x-1)/(x^2*(x+2))');
y4=sym('2*(cos(x))^2-(sin(x))^2');s1=expand(y1)
s2=factor(y2)[a,b]=numden(y3)s3=simple(y4)
运行结果:
s1=e^x*x+2*e^x+x^2+2*x
s2=(a-1)*(a^2+a+1)a=x^3+6*x^2+6*x-1b=x^2*(x+1)*(x+2)s3=3*cos(x)^2-1
④已知数学表达式f(x)=axn+bt+c,对其进行如下的符号替换:
1)a=sint,b=lnz,c=de2t的符号变量替换。
2)n=3,c=л的符号常量替换。
3)c=1:
2:
5替换。
17/16
ç24÷
4)c=æ13ö
è ø
的数组矩阵替换。
源程序:
symsaxbtcnsintlnzde2tf=sym('a*x^n+b*t+c')f1=subs(f,'a',sint);
f2=subs(f1,'b',lnz);
subs(f2,'c',de2t)
f3=subs(f,'n',3);
subs(f3,'c',pi)
subs(f,'c',1:
2:
5)
subs(f,'c',[1,3;2,4])
运行结果:
f=a*x^n+b*t+c
ans=sint*x^n+lnz*t+de2tans=a*x^3+b*t+pi
ans=[a*x^n+b*t+1,a*x^n+b*t+3,a*x^n+b*t+5]ans=
[a*x^n+b*t+1,a*x^n+b*t+3][a*x^n+b*t+2,a*x^n+b*t+4]
⑤已知符号表达式f
函数f(g(x))。
源程序:
symsx
=1-sin2x,g=2x+1,计算x=0.5时,f的值;计算复合
f=sym('1-(sin(x))^2');
g=sym('2*x+1');
f1=subs(f,'x',0.5)f2=compose(f,g)
运行结果:
f1=0.7702
f2=1-sin(2*x+1)^2
⑥求lim
x®2x
x2-1
2
。
-3x+1
源程序:
symsx
f=sym('(x^2-1)/(x^2-3*x+1)');
limit(f,x,2)运行结果:
ans=-3
ex+xsinx
⑦求函数f(x)=cos2x-sin2x的积分;求函数g(x)= 的导数。
源程序:
symsx
f=sym('cos(2*x)-sin(2*x)');
g=sym('sqrt(e^x+x*sin(x))');f1=int(f)
f2=diff(g)运行结果:
f1=1/2*sin(2*x)+1/2*cos(2*x)
f2=1/2/(e^x+x*sin(x))^(1/2)*(e^x*log(e)+sin(x)+x*cos(x))
ò
p
⑧计算定积分6(sinx+2)dx
0
源程序:
symsx
f=sym('sin(x)+2');
s1=int(f,x,0,pi/6)运行结果:
s1=-1/2*3^(1/2)+1/3*pi+1
⑨求下列线性代数方程组的解。
ìx+y+z=10
ï
í3x+2y+z=14
î
ï2x+3y-z=1
源程序:
symsxyz
f1=sym('x+y+z=10')
f2=sym('3*x+2*y+z=14')f3=sym('2*x+3*y-z=1')[s1,s2,s3]=solve(f1,f2,f3)
运行结果:
f1=x+y+z=10
f2=3*x+2*y+z=14
f3=2*x+3*y-z=1s1=1
s2=2
s3=7
⑩求解当y(0)=2,z(0)=7时,微分方程组的解。
ï
í
ìdy-z=sinxdx
dz
ï
îdx
+y=1+x
源程序:
symsxzy
[s1,s2]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2','z(0)=7','x')
运行结果:
s1=cos(x)+6*sin(x)+1/2*sin(x)*x+x+1
s2=-3/2*sin(x)+6*cos(x)+1+1/2*cos(x)*x
4MATLAB程序设计
实验目的:
掌握MATLAB程序设计的主要方法,熟练编写MATLAB函数。
实验内容:
①用π/4≈1-1/3+1/5-1/7+…公式求π的近似值,直到最后一项的绝对值小于10-6
为止,试编写其M脚本文件。
源程序:
方法一:
n=1;
sum1=0;
while(1/(2*n-1))>1.0000e-006sum1=sum1+4*(-1)^(n+1)*(1/(2*n-1));n=n+1;
endsum1
方法二:
sum1=0;
fork=1:
10^200
sum1=sum1+4*(-1)^(k+1)*(1/(2*k-1));
if(1/(2*k-1))<10^(-6)
break;end;
endsum1
运行结果:
sum1=3.1416
②分别用for和while结构计算1+22+33+…+100100的运行程序。
源程序:
While:
n=1;Sum2=0;
whilen<=100Sum2=sum2+n^n;n=n+1;
endSum2For:
Sum2=0;
forn=1:
100
Sum2=sum2+n^n;
endSum2
执行:
Sum2=1.0037e+200
5MATLAB绘图
实验目的:
掌握MATLAB二维图形绘制,掌握图形属性的设置和图形修饰;掌握图像文件的读取和显示。
实验内容:
①绘制图形
图1
px
其中x初值为0,终值为12.5,步长为0.1,第一个曲线y1为sin( ),第二
6
px
个曲线y2为cos( )。
写出图1的绘制源程序。
按照以下的步骤进行
(1)产生曲线
6
的数据(共有3组数据:
x,y1,y2);
(2)选择合适的线形、标记、颜色(正弦曲线为红色,余弦曲线为紫色);(3)添加图例及文字说明信息;(4)添加坐标轴说明与图标题。
源程序:
x=0:
0.1:
12.5;
y1=sin(pi/6*x);
y2=cos(pi/6*x);plot(x,y1,'r',x,y2,'--m')
title('两曲线比较图')x=[2/3*pi;6/3*pi];y=[0.42;0];
s=['sin(pi/6*x)';'cos(pi/6*x)'];text(x,y,s)legend('sinx','cosx')
xlabel('时间t')
ylabel('幅值')gridon
运行结果:
②在同一窗口不同坐标系里分别绘出y1=sinx,y2=cosx,y3=cinh(x),y4=cosh(x)4
个图形,其中x为以0为初值,2p为终值,元素数目为30的行向量。
源程序:
x=linspace(0,2*pi,
30);
y1=sin(x);y2=cos(x);y3=sinh(x);y4=cosh(x);subplot(2,2,1);
plot(x,y1');
subplot(2,2,2);
plot(x,y2');
subplot(2,2,3);
plot(x,y3');
subplot(2,2,4);
plot(x,y4');
运行结果:
③绘制一个三维曲线其中x=cosp,y=sinp,z=p,p为以0为初值,20p为终值,步长0.1p的行向量。
源程序:
p=0:
0.1*pi:
20*pi;x=cos(p);
y=sin(p);z=p;plot3(x,y,z);
运行结果:
三 课程设计体会
经过一周的学习,我对mutlab有了初步的认识和能够完成一些简单的程序。
在学习过程中,老师认
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 作业 电气工程 专业 试题 答案