实验报告matlab程序设计

《实验报告matlab程序设计》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

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

2、值。命令文件：%求二元一次方程a*x^2+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

3、);y=mod(a,100)-x;z=a-y-x;e=100*x+y+z/100;disp(e)【例3】定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。函数文件：文件名：convert程序：function[z,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);d

4、isp(['复数',num2str(x),'+',num2str(y),'i的指数为:',num2str(a),'，对数为：',num2str(b),'，正弦为：',num2str(c),'，余弦为：',num2str(d)])【例4】利用函数文件，实现直角坐标（x，y）与极坐标（ρ，θ）之间的转换。已知极坐标的矢径、极角分别为：，函数文件：文件名：transform程序：function[r,theta]=transform(x,y)r=sqrt(x*x+y*y);theta=atan(y/x);函数调用：程序：x=input('Pleas

5、einputx=:');y=input('Pleaseinputy=:');[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语句组Cend【例6】输入一个百分制成绩，要求输出成绩等级A,B,C,D,E。其中90～100分为A，80～89分为B，70～79分为C，6

6、0～69分为D，60分以下为E。要求：1、采用if语句实现。2、输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出错误信息。程序：x=input('请输入成绩')ifx>=90&x<=100disp('等级为A')elseifx>=80&x<=89disp('等级为B')elseifx>=70&x<=79disp('等级为C')elseifx>=60&x<=69disp('等级为D')elseifx>=0&x<=59disp('等级为E')elsedisp('成绩有误')end【例7】已知，当n=100时，求y的值。令s0=0，则对于

7、任意i，有si=si-1+1/（i*i),i=1,2,…,n。要求：采用for循环语句实现。程序：y=0;n=100;fori=1:ny=y+1/i^2;endy