Matlab编程实验报告.docx

1、Matlab编程实验报告(此文档为word格式，下载后您可任意编辑修改！) 中北大学Matlab编程实验报告班 级： 姓 名： 学 号： 2013 年 11 月2，实验一 Matlab运算基础一、实验目的1、熟悉启动和退出MATLAB的方法。2、熟悉MATLAB工作窗口的组成。3、掌握建立矩阵的方法，矩阵分析与处理。4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。二、实验内容1 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。z1=2*sin(85*pi180)(1+exp(2)z1 =0.2375（2），其中 x=2,1+2*i;-0.45,5x = 2

2、.0000 1.0000 + 2.0000i -0.4500 5.0000 z2=12*log(x+sqrt(1+x*x)z2 = 0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i（3）提示：利用冒号表达式生成a向量，求各点的函数值时用点乘运算。 a=-3.0:0.1:3.0a = Columns 1 through 9 -3.0000 -2.9000 -2.8000 -2.7000 -2.6000 -2.5000 -2.4000 -2.3000 -2.2000 Columns 10 through 18 -

3、2.1000 -2.0000 -1.9000 -1.8000 -1.7000 -1.6000 -1.5000 -1.4000 -1.3000 Columns 19 through 27 -1.2000 -1.1000 -1.0000 -0.9000 -0.8000 -0.7000 -0.6000 -0.5000 -0.4000 Columns 28 through 36 -0.3000 -0.2000 -0.1000 0 0.1000 0.2000 0.3000 0.4000 0.5000 Columns 37 through 45 0.6000 0.7000 0.8000 0.9000 1.

4、0000 1.1000 1.2000 1.3000 1.4000 Columns 46 through 54 1.5000 1.6000 1.7000 1.8000 1.9000 2.0000 2.1000 2.2000 2.3000 Columns 55 through 61 2.4000 2.5000 2.6000 2.7000 2.8000 2.9000 3.0000 z3=(exp(0.3.*a)-exp(-0.3.*a)2.*sin(a+0.3)z3 = Columns 1 through 9 0.4387 0.5072 0.5640 0.6089 0.6424 0.6648 0.6

5、766 0.6784 0.6709 Columns 10 through 18 0.6549 0.6313 0.6011 0.5652 0.5247 0.4805 0.4337 0.3854 0.3366 Columns 19 through 27 0.2881 0.2410 0.1962 0.1543 0.1162 0.0824 0.0535 0.0299 0.0120 Columns 28 through 36 -0.0000 -0.0060 -0.0060 0 0.0117 0.0288 0.0509 0.0775 0.1080 Columns 37 through 45 0.1418

6、0.1780 0.2159 0.2547 0.2934 0.3311 0.3669 0.3998 0.4289 Columns 46 through 54 0.4532 0.4719 0.4841 0.4891 0.4862 0.4748 0.4543 0.4243 0.3846 Columns 55 through 61 0.3350 0.2755 0.2061 0.1272 0.0392 -0.0574 -0.1619（4） ，其中t0：0.5：2.5提示：用逻辑表达式求分段函数值。t=0:0.5:2.5;z4=(t=0 & t=1 & t=2 & t A=100:999; B=find(

7、mod(A,21)=0)B = Columns 1 through 16 6 27 48 69 90 111 132 153 174 195 216 237 258 279 300 321 Columns 17 through 32 342 363 384 405 426 447 468 489 510 531 552 573 594 615 636 657 Columns 33 through 43 678 699 720 741 762 783 804 825 846 867 888 C=length(B)C = 43（2）建立一个字符串向量，删除其中的大写字母。 X=LIsa is A

8、BeAutiful Girl; Y=find(X=A&X X(Y)=X =sa is eutiful irl提示：利用find函数和空矩阵。5 由指令A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。（提示：find和sub2ind） A=rand(3,5)A = Columns 1 through 4 0. 686 0. 803 0. 736 0. 816 0. 799 0. 679 0. 663 0. 989 0. 444 0. 057 0. 531 0. 228 Column 5 0. 809 0. 208 0. 031

9、K=find(A0.5)K = 1 2 3 4 8 1012C,D=ind2sub(size(A),K)C = 1 2 3 1 2 1 3D = 1 1 1 2 3 4 46 先运行clear, format long, A=rand(3,3)，然后根据A写出两个矩阵：一个对角阵B，其相应元素由A的对角元素构成；另一个矩阵C，其对角元素全为0，而其余元素与对应的A阵元素相同。（提示：diag） clear format long A=rand(3,3)A = 0. 773 0. 168 0. 609 0. 333 0. 557 0. 421 0. 916 0. 562 0. 890 b=dia

10、g(A)b = 0. 773 0. 557 0. 890 B=diag(b)B = 0. 773 0 0 0 0. 557 0 0 0 0. 890 C=A-BC = 0 0. 168 0. 609 0. 333 0 0. 421 0. 916 0. 562 0三、实验报告及要求实验报告给出详细的实验过程、代码以及实验结果，必要时按照实验任务要求对结果进行分析。实验二 多项式及字符运算一、实验目的1、掌握多项式的表示与创建方法。2、掌握多项式的四则运算以及多项式求根方法。3、掌握字符串数组的创建及操作方法。二、实验内容1、已知，分别取和一个矩阵，求在处的值。 p=2 -1 0 3; x=2;

11、y=polyval(p,x)y = 15 x1=1 2;3 4; y1=polyval(p,x1)y1 = 4 15 48 1152、 有3个多项式试进行下列操作：求。（1） 求的根。 p1=1 2 4 0 5; p2=1 2; p3=1 2 3; R=conv(p2,p3)R = 1 4 7 6 p4=0 1 4 7 6 ;p5=p1+p4p5 = 1 3 8 7 11 r=roots(p5)r = -1.3840 + 1.8317i -1.3840 - 1.8317i -0.1160 + 1.4400i -0.1160 - 1.4400i（2） 当x取矩阵A的每一元素时，求P（x）的值。其

12、中 A=-1 1.2 -1.4;0.75 2 3.5;0 5 2.5; y=polyval(p5,A)y = 1.0e+003 * 0.0100 0.0382 0.0125 0.0223 0.0970 0.4122 0.0110 1.2460 0.1644 （3） 当以矩阵A为自变量时，求P（x）的值。其中A的值与（3）相同。X=AX = -1.0000 1.2000 -1.4000 0.7500 2.0000 3.5000 0 5.0000 2.5000 Y=X4+3*X3+8*X2+7*X1+11Y = 1.0e+003 * 0.0076 -0.1172 -0.0665 0.1437 1.

13、3900 1.1754 0.1934 1.7474 1.51983.用两种方法建立如下的二维字符串数组：第一行为“Matlab编程”，第二行为“实验环节”，最后利用该二维字符串形成一个一行字符串数组，“Matlab编程实验环节” A=Matlab编程;实验环节 A =Matlab编程实验环节 B=str2mat(Matlab编程,实验环节)B =Matlab编程实验环节 C=B(1,:),B(2,:)C =Matlab编程实验环节 三、实验报告及要求实验报告给出详细的实验代码以及实验结果，必要时按照实验任务要求对结果进行分析。实验三 Matlab符号计算一、实验目的1、熟练掌握MATLAB符号

14、表达式的创建2、熟练掌握符号表达式的代数运算3、掌握符号表达式的化简和替换4、熟练掌握符号微积分5、熟练掌握符号方程的求解二、实验内容1. 创建符号表达式：。 （用两种方式分别创建） f=sym(a*x3+b*x2+c*x+d) f = a*x3 + b*x2 + c*x + dsyms a b c d f=a*x3+b*x2+c*x+d f = a*x3 + b*x2 + c*x + d2. 创建符号矩阵。 A=sym(a*cos(x)+b*sin(y) 10+20;a*x2+b*y2+c*z2 sqrt(t*t+1) A = a*cos(x) + b*sin(y), 30 a*x2 + b

15、*y2 + c*z2, (t2 + 1)(12)3. 已知表达式，计算当时的值；计算与的复合函数，、的逆函数。f=1-sin(x)*sin(x); x=1x = 1 subs(f)ans = 0.2919 f=sym(1-sin(x)*sin(x); g=sym(2*x+1); compose(f,g) ans =1 - sin(2*x + 1)2 finverse(f)Warning: finverse(1 - sin(x)2) is not unique. ans = asin(1 - x)(12) finverse(g) ans = x2 - 12 4. 符号函数，分别对、进行微分，对趋

16、向于1求极限，并计算对的二次、三次微分，用findsym得出符号变量。syms x y a b c df=a*x3+b*y2+c*y+d; A=diff(f,x) A = 3*a*x2 B=diff(f,y) B = c + 2*b*y C=diff(f,c) C = y D=diff(f,d) D = 1 L=limit(a*x3+b*y2+c*y+d,y,1) L = a*x3 + b + c + d E=diff(f,x,2) E = 6*a*x F=diff(f,x,3) F = 6*a findsym(E)ans =a,x findsym(F)ans =a 5. 求极限。syms n

17、;L=limit(symsum(12n,n,0,inf),n,inf) L = 26. 求定积分和不定积分。syms x; f=sqrt(x)(1+x)2; int(f,x,1,inf) ans = pi4 + 12 syms x y; q=x*exp(-x*y);I=int(q) I = -(x*y + 1)(y2*exp(x*y) S=int(I) S = (x*y + 2)(y3*exp(x*y) 7. 求微分方程在满足，的解。y=dsolve(D2y+2*Dy+2*y,y(0)=1,Dy(0)=0,x) y = cos(x)exp(x) + sin(x)exp(x)8. 求解微分方程组

18、x,y=dsolve(Dx=2*x+3*y,Dy=x-2*y,x(0)=1,y(0)=2,t) x = (exp(7(12)*t)*(2*7(12) - 3)2 - (2*7(12) + 3)(2*exp(7(12)*t) + (7(12)*exp(7(12)*t)*(2*7(12) - 3)7 + (7(12)*(2*7(12) + 3)(7*exp(7(12)*t) y = (7(12)*exp(7(12)*t)*(2*7(12) - 3)14 + (7(12)*(2*7(12) + 3)(14*exp(7(12)*t) ，其中9. 解微分方程的通解。 y=dsolve(Dy+y*tanx

19、=cosx,x) y = (cosx - C11exp(tanx*x)tanx10、求解方程组x,y,z=solve(x+2*y-z=7,x+z=3,x2+3*y2=28,x,y,z) x = 154 - 37(12)4 37(12)4 + 154 y = 37(12)4 + 54 54 - 37(12)4 z = 37(12)4 - 34 - 37(12)4 - 34三、实验报告及要求实验报告给出详细的实验代码以及实验结果，必要时标注代码注释，并按照实验任务要求对结果进行分析。实验四 选择结构程序设计一、实验目的1、掌握建立和执行M文件的方法。2、掌握利用if语句实现选择结构的方法。3、掌握

20、利用switch语句实现多分支选择结构的方法。4、掌握try语句的使用。二 、实验内容1.求下列分段函数的值。 要求：用if语句实现，分别输出x= -5.0, -3.0, 1.0, 2.0, 2.5, 3.0, 5.0 时的y值。提示：x的值从键盘输入。 请输入x的值： -5x = -5y = 14请输入x的值： -3x = -3y = 11请输入x的值： 1x = 1y = 2请输入x的值： 2x = 2y = 1请输入x的值： 2.5x = 2.5000y = -0.2500请输入x的值： 3x = 3y = 5请输入x的值： 5x = 5y = 62输入一个百分制成绩，要求输出成绩等级A

21、,B,C,D,E 。其中90100分为A，8089分为B，7079分为C，6069分为D，60分以下为E。要求：（1）用switch语句实现。（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出错误信息。(1)请输入成绩:61grade = 61rate =D请输入成绩:54grade = 54rate =E请输入成绩:87grade = 87rate =B请输入成绩:78grade = 78rate =C(2)请输入成绩:103grade = 103成绩无效3、建立56矩阵，要求输出矩阵第n行元素。当n值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。请输入行数n

22、: 2n = 2 3 5 6 7 8 9请输入行数n: 6n = 6 4 6 8 9 0 3ans =Index exceeds matrix dimensions.三、实验报告及要求实验报告给出详细的实验代码以及实验结果，必要时按照实验任务要求对结果进行分析。实验五 循环结构程序设计一、实验目的1、掌握利用for语句实现循环结构的方法。2、掌握利用while语句实现循环结构的方法。3、熟悉利用向量运算来代替循环操作的方法。二、实验内容1.根据，求的近似值。当n分别取100、1 000、10 000时，结果是多少？要求：分别用循环结构和向量运算（使用sum函数）来实现。(1)请输入n值; 100n = 100s = 3.1321请输入n值; 1000n = 1000s = 3.1406请输入n值; 10000n = 10000s =3.1415(2)请输入n值; 100n = 100s = 3.1321请输入n值; 1000n = 1000s = 3.1406请输入n值; 10000n = 10000s = 3.14152.根据.，求：（1）y0)length(find(f=0)length(find(f0)（2） 正数、零、负数的个数。A = 4.