实验报告Matlab程序设计

1、学号1051020103班级数计1103姓名李东指导教师孙桂利实验题目实验二Matlab程序设计实验目的：1、掌握交互式操作命令，input、pause、disp、subplot、hold on；3、掌握建立和执行M文件的方法；4、掌握定义和调用MATLAB函数文件的方法，以及function命令的使用；6、了解MATLAB程序调试的基本方法；8、掌握利用语句if、for语句实现循环结构的方法。实验课时：6小时【例11求次方程ax2 bx c 0的根。采用input函数提示分别输入方程系数a、b、c 的2.值，令d b 4ac,则万程的根,.dx1b2a又2、. d2a 。采用disp函数显示

2、两根的值。实验内容:命令文件： %求二元一次方程a*xA2+b*x+c=0的根a=input('a=');b=input('b=');c=input('c=');d=b*b-4*a*c;x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a) disp('x1=',num2str(x(1),', x2=',num2str(x(2) %num2str表示数字化为字符串。【例2】将一个三位数分别拆成分别个位，十位，百位。然后交换顺序输出。命令文件：%一个三位数交换顺序输出a=input('请输

3、入一个三位数字：');x=mod(a,10);y=mod(a,100)-x;z=a-y-x;e=100*x+y+z/100;disp(e)【例3】定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数 文件。函数文件：文件名：convert程序：functionz,d,s,c=convert(a,b)z=exp(a+b*i);d=log(a+b*i);s=sin(a+b*i);c=cos(a+b*i);函数调用：程序：clearx=input('请输入复数实部的值');y=input('请输入复数虚部的值');a,b,c,d=co

4、nvert(x,y);disp('复数',num2str(x),'+',num2str(y),'i 的指数为：',num2str(a),',对数为：,num2str(b),',正弦为: ',num2str(c),',余弦为：',num2str(d)【例4】利用函数文件，实现直角坐标(x, y)与极坐标(p, 9)之间的转换。已知极坐标的矢y22 arctan径、极角分别为："x y ,x函数文件：文件名：transform程序：function r,theta=transform(x,y)r=sq

5、rt(x*x+y*y);theta=atan(y/x);函数调用：程序：x=input('Please input x=:');y=input('Please input y=:');r,theta=transform(x,y);disp('矢径 r 为:',num2str(r),'','极角 theta 为:',num2str(theta)【例5】if语句有几种表现形式？(1) if表达式语句组Aend(2) if表达式语句组Aelse语句组Bend(3) if表达式1语句组Aelseif表达式2语句组Belse

6、语句组Cend【例6】输入一个百分制成绩，要求输出成绩等级A,B,C,D,E其中90100分为A, 8089分为B, 7079分为C, 6069分为D, 60分以下为E。要求：1、采用if语句实现。2、输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出错误信息。程序：x=input('请输入成绩')if x>=90&x<=100dispC为 A')elseif x>=80&x<=89dispC为 B')elseif x>=70&x<=79dispC为 C')elseif x>=60&x<=69dispC为 D')elseif x>=0&x<=59dispC为 E')elsedisp(成绩有误')end,111y 12 -2 .2【例7】已知 23 n ,