《实验报告Matlab程序设计.docx》由会员分享，可在线阅读，更多相关《实验报告Matlab程序设计.docx（6页珍藏版）》请在三一文库上搜索。
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(请输入一个三位数字：);x=mod(a,10);y=mod(a,100)-x;z=a-y-x;e=100*x+y+z
3、/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=convert(x,y);disp(复数,num2str(x),+,num2str(y),i 的指数为：,num2str(a),对数为：,num2str(b),正弦为
4、: ,num2str(c),余弦为：,num2str(d)【例4】利用函数文件，实现直角坐标(x, y)与极坐标(p, 9)之间的转换。已知极坐标的矢y22 arctan径、极角分别为：x y ,x函数文件：文件名：transform程序：function r,theta=transform(x,y)r=sqrt(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 为
5、:,num2str(theta)【例5】if语句有几种表现形式？(1) if表达式语句组Aend(2) if表达式语句组Aelse语句组Bend(3) if表达式1语句组Aelseif表达式2语句组Belse语句组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=80&x=70&x=60&x=0&x=59dispC为 E)elsedisp(成绩有误)end,111y 12 -2 .2【例7】已知 23 n ,当n=100时，求y的值。令s0=O,则对于任意i,有si=si-1+/(i*i),i=1,2,，n要求：采用for循环语句实现 程序：y=0;n=100;for i=1:ny=y+1iA2;endy