数学建模MATLAB实验报告.docxVIP

学生实验报告

萨荼Satu

实验项目名称：MATLAB程序设计实验学时：

同组学生姓名：实验地点：

实验日期：实验成绩：

批改教师：批改时间：

一、实验目的和要求

1.掌握MATLAB软件的安装、配置和介绍等内容。

2.掌握初步使用MATLAB方法。

3.掌握数据和函数绘图命令。

4.熟悉MATLAB软件的主要界面、查阅帮助文档等内容。

二、实验仪器和设备

主机要求：安装WindowsXP、7、8、10等操作系统的客户端和安装MATLAB2014a以上版本的软件。

三、实验过程

3.1M文件与M函数

例1.1用M文件画出衰减震荡曲线y=e-t3cos3t以及它的包络线

t=linspace(0,4*pi,1000);%设置t的范围

y=exp(-t/3).*cos(3*t);

y0=exp(-t/3);

figure;

plot(t,y,b-,LineWidth,1.5);???%衰减振荡曲线

holdon;

plot(t,y0,r--,LineWidth,1.2);??%上包络线

plot(t,-y0,r--,LineWidth,1.2);?%下包络线

holdoff;

xlabel(t);

ylabel(y);

title(衰减振荡曲线及其包络线);

legend(y=exp(-t/3)cos(3t),y_0=exp(-t/3),-y_0=-exp(-t/3));

gridon;

例1.2M函数文件示例，随机方阵A，同时求A,A2，A-1

functionrand_matrix_operations()

??n=4;

??A=rand(n);

??disp(随机方阵A=);

??disp(A);

??%行列式

??detA=det(A);

??disp(行列式|A|=);

??disp(detA);

??%平方

??Asquared=A^2;

??disp(A^2=);

??disp(Asquared);

??%逆矩阵

??invA=inv(A);

??disp(A^{-1}=);

??disp(invA);

??%转置

??AT=A;

??disp(A=);

??disp(AT);

end

rand_matrix_operations

3.2程序结构

例1.3一个典型的顺序文件：

disp(请看执行结果：)

disp(thebeginoftheprogram)

disp(thefirstline)

disp(thesecondline)

disp(thethirdline)

disp(theendoftheprogram)

例1.4利用for循环求1-250的整数之和。

functionsum_1_to_250()

??total=0;

??fori=1:250

????total=total+i;

??end

??disp([1到250的整数之和为:,num2str(total)]);

end

例1.5利用for循环找出250-500之间的所有素数。

functionprimes_250_500()

??prime_list=[];

??forn=250:500

????ifisprime(n)

??????prime_list=[prime_list,n];

????end

??end

??disp(250到500之间的所有素数为:);

??disp(prime_list);

end

例1.6利用while循环计算1！+2！+3！+...+20!的值。

functionsum_factorial_1_to_20()

??n=1;?

??sum_val=0;????

??fact=1;?????

??whilen=20

????fact=fact*n;?

????sum_val=sum_val+fact;

????n=n+