MATLAB大作业MALAB大作业.doc

贵州大学实验报告 姓名学号实验组实验时间指导教师成绩实验项目名称MATLAB大作业实验目的通过综合考察，让同学们能对学到的知识综合运用，更好地灵活地解决问题实验原理利用matlab提供的一些函数实现数据分析和数值计算。实验仪器运行MATLAB7.0或以上的计算机一台实验内容 实验数据试编写名为test01.m的MATLAB函数，用以计算下述的值： 要绘制t关于函数f(t)的图形，其中t的取值范围为，间距为。 （提示：注意要产生一系列的点，这里可考虑t的输入是向量形式，可以利用find函数找出所需限定值的元素的位置，对其按需要赋值后，再进行绘图；其次，另外一种思路，也可考虑使用循环的形式来实现） test01.m代码如下： function y=test01(t) if sin(t)sin(pi/4) y=sin(pi/4); elseif sin(t)sin(-pi/4) y=sin(-pi/4); else y=sin(t); end 调用test01函数计算值 绘制t关于函数f(t)的图形，其中t的取值范围为，间距为 代码如下： t=-6*pi:pi/10:6*pi; y=sin(t); y(find(ysin(pi/4)))=sin(pi/4);%f(t)=sin(pi/4),ysin(pi/4) y(find(ysin(-pi/4)))=sin(-pi/4);%f(t)=sin(pi/4),ysin(pi/4) plot(t,y) %绘制图形 title(t from -6{\pi} to 6{\pi}) %添加图形标题 xlabel(t) %x轴标签 ylabel(f(t)) %y轴标签 实验结果： 编写函数，在同一窗口的4个子图中利用plot等语句绘制y=at2图像，其中a=[1 2 5 10],t QUOTE ∈（-2,5） 范围[-2，5]。 代码如下： a=[1,2,5,10]; t=-2:0.01:5; %生成t值抽样值 y1=a(1)*t.^2; %计算y=t2值 y2=a(2)*t.^2; %计算y=2t2值 y3=a(3)*t.^2; %计算y=5t2值 y4=a(4)*t.^2; %计算y=10t2值 subplot(2,2,1) %分割为2*2四个子窗口，活动窗口为1 plot(t,y1) % 绘制y=t2图象 title(y=t^2) %添加图象标题 subplot(2,2,2) plot(t,y2) %绘制y=2t2图象 title(y=2t^2) subplot(2,2,3) plot(t,y3) %绘制y=5t2图象 title(y=5t^2) subplot(2,2,4) plot(t,y4) %绘制y=10t2图象 title(y=10t^2) 实验结果： 求函数在区间[-5,5]上的最大值和最小值。 ㈠ 方法一： 代码如下: x=-5:0.01:5; %生成x的抽样值 y=x.^3.+2.*x.^2.-3; %计算抽样x对应的函数值 MAX=max(max(y)) %找到最大值 MIN=min(min(y)) %找到最小值 实验结果 ㈡ 方法二： 对于方法一，虽然大概能找到最大值和最小值，但那都只是近似最大值和最小值，因为对x 的选取都是离散的，事先并不敢保证取最大值和最小值的x就在抽样的队列当中（当然本题比较特殊，取最大值和最小值的x在抽样的x队列中），下面利用更一般的解法求解连续可微函数在区间上的最大值和最小值问题 原理：对于连续可微函数而言，最大值和最小值只会出现在函数的极值点处，因此本题还可以利用fminbnd函数来求得最小值和最大值（求最大值时需要一定的转换）。另外，由于直接利用fminbnd(@fname,x1,x2,options)函数计算，则只能针对开区间（-5,5）,而本题给出一闭区间[-5，5]，因而对得到的极值点都还要和函数在x=-5和x=-5处的函数值比较，只有这样才能得到真正的最大值和最小值。 另外，由于matlab中没有提供专门的函数求最大值，但我们可以求-f（x）上最小值就是f（x）上的最大值了，本题就采用这种转换求得最大值 利用fminbnd函数计算最大值和最小值代码如下： format long %将结果显示长度设置为long型，便于对结果对比 y1=(-5)^3+2*(-5)^2-3; %计算在x=-5处的函数值 y2=(5)^3+2*(5)^2-3; %计算在x=5算的函数值 [x,mi]=fminbnd(@(x)(x.^3.+2.*x.^2.-3),-5,5) %计算极小值 [x,ma]=fminbnd(@(x)