课程设计报告MATLAB江苏大学

江苏大学《电气工程工具软件培训》课程设计报告设计题目：MATLAB工具软件专业班级：电气1201：学生学号：：完成日期：2013.7.3江苏大学•电气信息工程学院（a组指导老师:黄永红；b组指导老师：刘辉；c组指导老师:王博）一MATLAB课程设计的目的和要求MATLAB软件功能简介MATLAB的名称源自MatrixLaboratory,1984年由美国Mathworks公司推向市场。它是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而被广泛的应用于科学计算、控制系统和信息处理等领域的分析、仿真和设计工作。MATLAB软件包括五大通用功能，数值计算功能（Nemeric）、符号运算功能（Symbolic）、数据可视化功能（Graphic）＞数字图形文字统一处理功能（Notebook）和建模仿真可视化功能（Simulink）。其中，符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLAB命令窗口。该软件有三大特点，一是功能强大；二是界面友善、语言自然；三是开放性强。目前，Mathworks公司已推出30多个应用工具箱。MATLAB在线性代数、矩阵分析、数值及优化、数值统计和随机信号分析、电路与系统、系统动力学、次那好和图像处理、控制理论分析和系统设计、过程控制、建模和仿真、通信系统以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。MATLAB课程设计的目的本次课程设计主要是为了使学生了解MATLAB软件的基本知识，熟悉MATLAB的上机环境，掌握MATLAB数值运算、程序设计、二维/三维绘图、符号运算、Simulink仿真等相关知识，并初步具备将一般数学问题转化为对应的计算机进行处理的能力，以便为今后进一步的学习打下坚定基础。二MATLAB课程内容1MATLAB语言基础实验目的：基本掌握MATLAB向量、矩阵、数组的生成及其基本运算(区分数组运算和矩阵运算)、常用的数学函数。了解字符串的操作。实验内容：①创建以下矩阵：A为初值为1，终值为12，元素数目为6的行向量；(2436)1532B=2779；C为B的三次万矩阵；D由B和C横向拼接并去除第一列、"1235)最后一列和第一行元素而成；E由B和C纵向拼接而成；F抽取E的3〜5行和第2〜3列元素生成；G由F经变形为3x4的矩阵而得；H由B和C数组相乘运算而成，同时将H(1,1)和H(2,1)分别变为n的平方和立方，H(2,2)=arccos(2)，H(3,3)=H(1,1)+H(2,1)。源程序：>>A=[1,3,2,5,6,12]>>B=[2,4,3,6;1,5,3,2;2,7,7,9;1,2,3,5]>>C=B*B*B>>M=[BC]>>D=M(2:4,2:7)>>E=[B;C]>>F=E(3:5,2:3)>>G=repmat(F,1,2)>>H=B*C>>H(1,1)=piA2>>H(2,1)=piA3>>H(2,2)=acos⑵>>H(3,3)=H(1,1)+H(2,1)运行结果：A=1B=3256122436153227791235C=2929388731140227751681868518167515552021227724681895M=Columns1through624362929381532227751277951816751235227724Columns7through8

TOC\o"1-5"\h\z873114068186815552021681895D=5322277516817795181675155523522772468124361532277912352929388731140227751681868518167515552021227724681895TOC\o"1-5"\h\z7723938873G=7777232393887393887344081424913221171853435111661030513333784225374235273055834351108510305134141.0e+004*0.00101.42491.32211.71850.34351.11661.03051.33330.78422.53742.35273.05580.34351.10851.03051.34141.0e+004*0.00100.00310.78420.00100.00310.78420.34351.42491.11662.53741.10851.32211.71851.03051.33332.35273.05581.03051.34141.0e+004*0.00101.42491.32211.71850.00310+0.0001i1.03051.33330.78422.53740.00413.05580.34351.10851.03051.3414H=1.0e+004*0.00101.42491.32211.71850.00310+0.0001i1.03051.33330.78422.53740.00413.05580.34351.10851.03051.3414②(1)用矩阵除法求下列方程组的解村[明；知灼；6x+3x+4x=3〈一2x+5x+7x=-48x一x一3x=-7V123⑵求矩阵的秩(rank函数)；求矩阵的特征值与特征向量(eig函数)；系数矩阵的3次幕与开方；系数矩阵的指数运算和数组对数运算；系数矩阵a(1,2)、a(1,3)、a(2,2)、a(2,3)的元素不变，其余元素变为零。提取系数矩阵主对角线上的元素，并依次相加赋予b。源程序：>>a=[6,3,4;-2,5,7;8,-1,-7]>>b=[3;-4;-7]>>x=a\b>>c=rank(x)»lambda=eig(a)»d=aA3»e=sqrtm(a)»f=expm(a)»g=log(a)»a(l,l)=O;»a(2,l)=0;»a(3,l)=0;»a(3,2)=0;»a(3,3)=0»a=[6,3>4;-2,5,7;8,-l,-7]»b=a(l,l)+a(22)+a(3>3)运行结果：a=TOC\o"1-5"\h\z634-2578-1-7b=3x=0.8196-3.97942.5052c=lambda=-7.74878.95192.7968d=TOC\o"1-5"\h\z450314332418437450438-360

2.44950+2.44950+1.4142i2.8284e=2.23612.64582.4144+0.2615i-0.2367+0.9088i1.3810-1.5804if=0.6223-0.0987i2.4144+0.2615i-0.2367+0.9088i1.3810-1.5804if=0.6223-0.0987i2.0722-0.3431i-0.0883+0.5966i0.7573-0.4741i1.1524-1.6476i0.1778+2.8652i1.0e+003*5.26543.28822.76211.80101.14950.95902.52931.57441.3238g=1609/898713/6491588/2291+355/113i1603/9964319/207702731/1970+355/113i1475/7581475/758+355/113ia=034057000a=634-2578-1-7b=42MATLAB数值运算实验目的：掌握MATLAB的数值运算及其运算中所用到的函数，掌握结构数组的操作。实验内容：①已知多项式a(x)=x2+2x+3；b(x)=4x2+5x+6求多项式a(x)和多项式b(x)的乘法运算结果，并在命令窗口中显示该多项式c；求多项式c的根及其微分；源程序：>>p1=[1,2,3];>>p2=[4,5,6];>>c=conv(p1,p2)>>c=poly2sym(c)>>r=roots(c)>>q=polyder(c)运行结果：c=413282718c=4*xA4+13*xA3+28*xA2+27*x+18r=-1.0000+1.4142i-1.0000-1.4142i-0.6250+1.0533i-0.6250-1.0533iq=16395627②求(睛+1)((+3)((+D的"商，，及“余,，多项式并在命令窗口中显示该多项式。s3+2s+1源程序：>>formatratp1=conv([1,0,1],conv([1,3],[1,1]));p2=[1,0,2,1];»formatrat»pl=conv([l,0,U»conv([l,3],[l,l]));»p2=[l,0,2,U；»[q,r]=deconv(pl,p2);»cq=，商多项式为"《*=涂多项式为"»disp([cq,poly2str(q,*s,)]),disp([cr,poly2str(r,'s')])运行结果：商多项式为s+4余多项式为2sA2-5s-l③⑴计算当x=2,x=3时，y(x)=X3+(^~0-98)2_5(x+-)的值;3+1.25)3x计算cos60+arccos(7i)-V9-^2的值‘2436、a=1532,b=A2+3,C=A-2B,，求：C2779J235?源程序;»symsx»f=xA3+(x-0.98)A2/(x+1.25)A3-5*(x+Vx)»fl=subs(f,,2,)»answ=vpa(fl,6)»f2=subs(f,,3,)»answ=vpa(f2,6)»A=[2,4,3,6;l,5,3,2;2,7,7,9;l,2,3,5]»B=AA2+3»c=A-2*B运行结果：f=xA3+(x-49/50)A2/(x+5/4)A3-5*x-5/xn=⑵A3+((2)-49/50)A2/(⑵+5/4)人3-5*⑵-5/(2)answ=

-4.46969f2=(3)A3+((3)-49/50)A2/((3)+5/4)A3-5*(3)-5/(3)answ=10.3865(3)2121B=)A=45723373629523183718645711348604810648805613765c=-44-35-72-35-124-109-219-94-117-93-205-93-154-110-265-1253MATLAB符号运算实验目的：掌握符号变量和符号表达式的创建，掌握MATLAB的symbol工具箱的一些基本应用。实验内容：①已知f=(ax2+bx+c-3)3-a(cx2+4bx-1)，按照自变量x和自变量a，对表达式f分别进行降幕排列(同幕合并)。源程序：>>symsxabc%定义符号变量>>f=(a*xA2+b*x+c*3)A3-a*(c*xA2+4*b*x-1)>>fx=collect(f)%对f按x的降幂排列>>fa=collect(f,a)%对f按a的降幕排列运行结果：f=(a*xA2+b*x+3*c)A3-a*(c*xA2+4*b*x-1)fx=aA3*xA6+3*b*aA2*xA5+(3*c*aA2+2*bA2*a+a*(6*c*a+bA2))*xA4+(12*c*b*a+b*(6*c*a+bA2))*xA3+(3*c*(6*c*a+bA2)+6*bA2*c+9*a*cA2-c*a)*xA2+(27*cA2*b-4*b*a)*x+27*cA3+afa=aA3*xA6+3*(b*x+3*c)*xA4*aA2+(3*(b*x+3*c)A2*xA2-c*xA2-4*b*x+1)*a+(b*x+3*c)A3②已知f1=1/(a-b),f2=2a/(a+b),f3=(a+1)(b-1)(a-b)，分别求fl和f2的符号和、fl和f3的符号积、f1和f3的符号商。源程序：>>symsab;>>f1=1/(a-b)>>f2=2*a/(a+b)>>f3=(a+1)*(b-1)*(a-b)>>h1=f1+f2>>h2=f1*f3>>h3=f1/f3运行结果：f1=1/(a-b)f2=2*a/(a+b)f3=(a+1)*(b-1)*(a-b)h1=1/(a-b)+2*a/(a+b)h2=(a+1)*(b-1)h3=1/(a-b)A2/(a+1)/(b-1)③对下列表达式进行符号运算⑴已知数学表达式y(x)=(ex+x)(x+2),将其展开。已知数学表达式y(x)=a3-1,对其进行因式分解。⑶已知数学表达式心)=上二+工^，对其进行通分。x(x+1)x2(x+2)⑷已知数学表达式y(x)=2cos2x-sin2x，对其进行化简。源程序：>>symsx;>>y=(exp(x)+x)*(x+2);>>y=expand(y)>>symsa;>>y=a人3-1;>>y=factor(y)⑶>>y=sym('(x+3)/x*(x+1)+(x-1)/xA2(x+2)*)>>[n,d]=numden(y)>>symsx;>>y=2*cos(x)A2-sin(x)A2;>>simple(y)运行结果：y=exp(x)*x+2*exp(x)+xA2+2*xy=(a-1)*(aA2+a-1)y=(x+3)/x*(x+1)+(x-1)/xA2(x+2)4*xA2+xA3+4*x-1d=xA2simplify:3*cos(x)人2-1radsimp:2*cos(x)人2-sin(x)人2combine(trig):3/2*cos(2*x)+1/2factor:2*cos(x)人2-sin(x)人2expand:2*cos(x)人2-sin(x)人2combine:3/2*cos(2*x)+1/2convert(exp):2*(1/2*exp(i*x)+1/2/exp(i*x))A2+1/4*(exp(i*x)-1/exp(i*x))A2convert(sincos):2*cos(x)A2-sin(x)A2convert(tan):2*(1-tan(1/2*x)A2)A2/(1+tan(1/2*x)A2)A2-4*tan(1/2*x)A2/(1+tan(1/2*x)A2)A2collect(x):2*cos(x)人2-sin(x)人2ans=3*cos(x)人2-1④已知数学表达式f(x)=axn+bt+c,对其进行如下的符号替换:a=sint,b=lnz,c=de2t的符号变量替换。n=3,c=^的符号常量替换。c=1:2:5替换。(1a、…......c=的数组矩阵替换。124)源程序：(1)>>a=sym('sint')>>b=sym(*lnz*)>>c=sym(*de2t*)>>symsabxt;>>f=a*xAn+b*t+c;>>n=sym('3')>>c=sym(*pi*)>>symsabcxtn;>>f=a*xAn+b*t+c;>>c=sym('1:2:5')c=sym('[1,2;3,4]‘)运行结果：(1)a=sintb=lnzc=de2t(2)n=pi(3)c=1:2:5(4)c=[1,2][3,4]已知符号表达式f=1-sin2x,g=2x+1,计算x=0.5时，f的值；计算复合函数fg(x))。源程序：>>symsx;>>f=1-sin(x)A2>>g=2*x+1;>>f1=subs(f,*0.5,)>>fg=compose(f,g)运行结果：f=1-sin(x)A2f1=1-sin((0.5))A2fg=1-sin(2*x+1)A2求lim—X2__\—。xr2X2—3x+1源程序：>>symsx;>>ans=limit((xA2-1)/(xA2-3*x+1),x,2)运行结果：ans=-3求函数f(x)=cos2x-sin2x的积分；求函数g(x)=ex+xsinx的导数。源程序：»symsxyzab»S=cos(2*x)-sin(2*x)；»int(S)»symsxytuvzab»S=sqrt(exp(x)+x*sin(x))；»diff(S)运行结果：ans=V2*sin(2*x)+l/2*cos(2*x)(2)ans=l/2/(exp(x)+x*sin(x))A(l/2)*(exp(x)+sin(x)+x*cos(x))计算定积分f6(sinx+2)dx0源程序：»symsxyzab»S=sin(x)+2；»int(S,0,pi/6)运行结果：ans=-1/2*3^(1/2)+l/3*pi+l求下列线性代数方程组的解。x+y+z=10<3x+2y+z=142x+3y-z=1源程序：»[x,y,z]=solve('x+y+zTO','3*x+2*y+zT4','2*x+3*y-zT')运行结果：

7⑩求解当y(0)=2,z(0)=7时，微分方程组的解。'dy.——-z=sinxdxdz1、dx源程序：—+y=1+x>>[y,z]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2,x(0)=7','x')；Warning:Explicitsolutioncouldnot、dx源程序：>InC:¥MATLAB6p5¥toolbox¥symbolic¥dsolve.matline326>>[y,z]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2,z(0)=7','x')；>>pretty(y)>>pretty(z)运行结果：cos(x)+6sin(x)+1/2xsin(x)+1+x-3/2sin(x)+6cos(x)+1+1/2cos(x)x4MATLAB程序设计实验目的：掌握MATLAB程序设计的主要方法，熟练编写MATLAB函数。实验内容：①用n74-1-1/3+1/5-1/7+^公式求n的近似值，直到最后一项的绝对值小于10-6为止，试编写其M脚本文件。源程序：>>t=1；pi=0；n=1；s=1；>>whileabs(t)>1e-7pi=pi+t；n=n+2；s=-s；t=s/n；end>>>>pi=pi*4；>>fprintf('pi=%f/n',pi)运行结果：pi=3.141592/n②分别用for和while结构计算1+22+33+・・・+100100的运行程序。源程序：>>clearsum=0；fori=1：100sum=sum+i"i；end>>sum>>clearall；sum=0；i=1；whilei<=100sum=sum+「i；i=i+1；end>>sum执行：sum=1.0037e+200sum=1.0037e+2005MATLAB绘图实验目的：掌握MATLAB二维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。