MATLAB书上的一些实验报告和习题
MATLAB实验报告
实验一MATLAB操作基础实验目的和要求：1、熟悉MATLAB的操作环境及基本操作方法。
2、掌握MATLAB的搜索路径及设置方法。
3、熟悉MATLAB帮助信息的查阅方法实验内容：1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试验用help命令能否查询到自己的工作目录。
2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。
例1-1例1-2例1-3例1-43、利用帮助功能查询inv、plot、max、round等函数的功能。
4、完成下列操作：（1）在matlab命令窗口输入以下命令：x=0:pi/10:2*pi;y=sin(x);（2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。
5、访问mathworks公司的主页，查询有关MATLAB的产品信息。
主要教学环节的组织：教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。
思考题：1、如何启动与退出MA TLAB集成环境？启动：（1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动MATLAB系统。
（2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。
（3）在桌面上建立快捷方式后。
双击快捷方式图标，启动MA TLAB。
退出：（1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。
（2）在MA TLAB命令窗口中输入exit或quit命令。
（3）单击MATLAB主窗口的关闭按钮。
2、简述MATLAB的主要功能。
MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。
3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？使用‘；’隔开4、Help命令和lookfor命令有何区别？help是显示matlab内置的帮助信息一般是help 命令；而lookfor 关键词是通过关键词查找，“关键词”就是你要找的词语5、在MA TLAB环境下，建立了一个变量fac，同时又在当前目录下建立了一个m文件fac.m，如果需要运行fac.m文件，该如何处理？实验二MATLAB矩阵及运算实验目的和要求：1、掌握MATLAB数据对象的特点及运算规则2、掌握MATLAB建立矩阵的方法及矩阵处理的方法3、掌握MATLAB分析的方法实验内容：1.求下列表达式的值（1）（2）（3）（4）2.已知A=[-1,5,-4;0,7,8;3,61,7],B=[8,3,-1;2,5,3;-3,2,0] 求下列表达式的值：（1）A+6B和A^2-B+I(其中I为单位矩阵)。
MATLAB实验报告
MATLAB实验报告一、实验目的本次 MATLAB 实验旨在深入了解和掌握 MATLAB 软件的基本操作和应用，通过实际编程和数据处理，提高解决问题的能力，培养编程思维和逻辑分析能力。
二、实验环境本次实验使用的是 MATLAB R2020a 版本，运行在 Windows 10 操作系统上。
计算机配置为英特尔酷睿 i5 处理器，8GB 内存。
三、实验内容（一）矩阵运算1、矩阵的创建使用直接输入、函数生成和从外部文件导入等方式创建矩阵。
例如，通过｀1 2 3; 4 5 6; 7 8 9` 直接输入创建一个 3 行 3 列的矩阵；使用｀ones(3,3)｀函数创建一个 3 行 3 列元素全为 1 的矩阵。
2、矩阵的基本运算包括矩阵的加减乘除、求逆、转置等。
例如，对于两个相同维度的矩阵｀A` 和｀B` ，可以进行加法运算｀C ＝ A ＋ B` 。
3、矩阵的特征值和特征向量计算通过｀eig` 函数计算矩阵的特征值和特征向量，加深对线性代数知识的理解和应用。
（二）函数编写1、自定义函数使用｀function` 关键字定义自己的函数，例如编写一个计算两个数之和的函数｀function s ＝ add(a,b) s ＝ a ＋ b; end` 。
2、函数的调用在主程序中调用自定义函数，并传递参数进行计算。
3、函数的参数传递了解值传递和引用传递的区别，以及如何根据实际需求选择合适的参数传递方式。
（三）绘图功能1、二维图形绘制使用｀plot` 函数绘制简单的折线图、曲线等，如｀x ＝ 0:01:2pi; y ＝ sin(x)； plot(x,y)｀绘制正弦曲线。
2、图形的修饰通过设置坐标轴范围、标题、标签、线条颜色和样式等属性，使图形更加清晰和美观。
3、三维图形绘制尝试使用｀mesh` 、｀surf` 等函数绘制三维图形，如绘制一个球面｀x,y,z ＝ sphere(50)； surf(x,y,z)｀。
（四）数据处理与分析1、数据的读取和写入使用｀load` 和｀save` 函数从外部文件读取数据和将数据保存到文件中。
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。
矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。
MATLAB综合性实验报告8
一、实验目的：1.学会运用Matlab解决实际问题2.进一步掌握Matlab的一些基本操作3.通过范例体会Matlab的初步建模过程4学会用做动画二．实验仪器、设备或软件：电脑、MATLAB软件三．实验内容：1.已知前两年的猪肉的产量和猪肉的价格分别为：39吨，28吨，12元/公斤，17元/公斤，根据前一年的猪肉价格影响后一年猪肉产量当年猪肉产量影响当年猪肉价格的线性关系，编写程序，利用动画原理画出前十年猪肉的产量——价格的动态图形。
（参数设置为，c1=39,c2=28,c3=36,r1=12,r2=17,k=16)。
2．在地面上建有一座圆柱形水塔,水塔内部的直径为d,并且在地面处开了一个高为H的小门.现在要对水塔内部进行维修施工,施工方案要求把一根长为l(l>d)的水管运到水塔内部.请问水塔的门高H多高时,才有可能成功地把水管搬进水塔内。
四．实验记录：1．蛛网模型：实验代码：clear%c1为产量1, c2为产量2, c3为产量3, r1为%肉价1, r2为肉价2, k为K年后产量与肉价%是否稳定c1=39;c2=28;c3=36;hold offr1=12;r2=17;k=16;a1=[c1 1;c2 1];%系数矩阵b1=[r1,r2]';%列向量a2=[r1 1;r2 1];b2=[c2,c3]';a=a1\b1;b=a2\b2;%x0(1)=c1;for n=1:30y0(n)=a(1)*x0(n)+a(2); %a(1)a(2)为矩阵元素的引用，下行类似x0(n+1)=b(1)*y0(n)+b(2);x(n)=x0(n);y(n)=x0(n+1);endplot(x,y0,'-g',y,y0,'-b')hold onm=moviein(100);for n=1:kfor j=1:30t1=x0(n)+(j-1)*(x0(n+1)-x0(n))/30;t2=x0(n)+j*(x0(n+1)-x0(n))/30;if t2<t1t=t1;t1=t2;t2=t;%这样做方便画图elseendt=t1:0.01:t2;plot(t,y0(n),'.r')%划横线t1=y0(n)+(j-1)*(y0(n+1)-y0(n))/30;t2=y0(n)+j*(y0(n+1)-y0(n))/30;if t2<t1t=t1;t1=t2;t2=t;elseendt=t1:0.01:t2;plot(x(n+1),t,'.r')%划竖线endm(:,n)=getframe;endmovie(m,20)实验结果：ans =Columns 1 through 639.0000 28.0000 36.0000 30.1818 34.4132 31.3358 Columns 7 through 1233.5739 31.9462 33.1300 32.2691 32.8952 32.4398 Columns 13 through 1832.7710 32.5302 32.7053 32.5779 32.6706 32.6032 Columns 19 through 2432.6522 32.6166 32.6425 32.6236 32.6374 32.6274 Columns 25 through 3032.6346 32.6294 32.6332 32.6304 32.6324 32.6310 Column 3132.6320实验图像：2．水塔模型：实验代码：Clearfunction x=lt523(l,d,h) k1=0;a=0l=7;d=3ymax=0;k2=0;b=pi/2;h=3;while (b-a)>10^(-2)k2=k2+1;m=0;a=0;if k1==0n=ceil(b/h)-1 elsen=ceil(b/h);endfor i=1:n+1;x(i)=a+(i-1)*(b-a)/n;y(i)=l*sin(x(i))-d*tan(x(i));endfor i=1:nif y(i)>ymaxymax=y(i);a=x(i-1);else y(i)<ymaxb=x(i);endendm=m+1;k1=k1+1;h=(b-a)/5;endabx=(a+b)/2t=0:0.01:1.35y1=l*sin(t)-d*tan(t);plot(t,y1,'-b')方法：2syms h xh=8*sin(x)-2*tan(x);d1=0;d2=pi/2;d0=(d1+d2)/2;d=0.03;a=subs(h,x,d1);b=subs(h,x,d2);c=subs(h,x,d0);for n=d1:d:d2if (((c-a)/(d0-d1))*((b-c)/(d2-d0)))<0 d1=d1+d; d2=d2-d; d0=(d1+d2)/2;a=subs(h,x,d1);b=subs(h,x,d2);c=subs(h,x,d0);endendd0hold onaxis([0,pi/2,-5,5])x=0:0.05:pi/2;plot(x,8*sin(x)-2*tan(x),'-b')实验运行结果：五、实验总结：通过做此实验，让我对MATLAB有更进一步的了解，学会怎样才能正确运用MATLAB求解实际问题，了解如何利用数学模型去解释和分析社会经济问题，特别是这个典型经济问题的求解。
MATLAB实验报告
MATLAB实验报告专业班级学号姓名实验一Matlab环境及命令窗口的使用实验目的与要求：熟练掌握MA TLAB的启动和退出；熟悉MATLAB的命令窗口；熟悉常用选单和工具栏；熟悉MA TLAB桌面的其它窗口。
实验内容：启动MATLAB，熟悉各环境窗口及作用，熟悉file、edit、View等菜单下各选项的功能；熟悉工具栏中各项的功能；打开演示程序（demo），开始其中的演示程序；退出MA TLAB（exit、quit）。
在命令窗口中，输入简单的命令操作；打开函数编辑窗，编写简单的程序并保存；程序：>> x1=sqrt(100),x2=5.48,y=9/x2x1 =10x2 =5.4800y =1.6423>>demo演示框左侧是库目录，选定demo下的图形类（Graphics）按钮，双击演示框右边的栏里的function of complex variable选项，就会出现如图中的演示框。
此例为复数Z的三次方：Z^3。
实验二Matlab数值计算实验目的与要求：熟练掌握matlab变量的使用；熟练掌握矩阵的创建；熟悉掌握matlab 的矩阵和数组的运算；熟悉matlab多项式的运算。
实验内容：给2×3阶矩阵变量a赋值（注意逗号与分号的使用），取其第2行第2个元素，取其第2行元素，取其第1列和第3列的元素；计算a/3、5*a、a+2等的值；程序：>> a=[7,9,5;5,9,4]a =7 9 55 9 4>> a=[7,9,5;5 9 4]a =7 9 55 9 4>> b=a([2],[2])b =9>> c=a([2],:)c =5 9 4>> d=a(:,[1,3])d =7 55 4>> e=a/3e =2.33333.0000 1.66671.6667 3.0000 1.3333>> 5*aans =35 45 2525 45 20>> a+2ans =9 11 77 11 62、生成2×3阶的全1矩阵x和3×3魔方矩阵y，计算a+x,x*y，x/y，y^2,2^a;a.*x, a.\x，y.^2,2.^a 的值。
MATLAB实验报告3
MATLAB实验报告3MATLAB实验报告3一、实验目的1.掌握MATLAB程序的调试方法；2.掌握MATLAB中的矩阵操作；3.熟悉MATLAB中处理图像的基本操作。
二、实验内容1.用MATLAB调试程序；2.用MATLAB进行矩阵运算；3.用MATLAB处理图像。
三、实验原理及步骤1.MATLAB程序的调试方法在MATLAB中调试程序可以采用设置断点、逐行运行、单步调试等方法。
设置断点可以在程序中的其中一行上点击左键，会出现一个红色的圆点表示断点已设置。
逐行运行可以通过点击Editor界面上的运行按钮实现。
单步调试可以通过点击断点所在行的左侧按钮实现。
2.矩阵运算在MATLAB中，对于矩阵的运算可以使用一些基本的函数，如矩阵加法、减法、乘法等。
矩阵加法可以使用"+"操作符实现，减法可以使用"-"操作符实现，乘法可以使用"*"操作符实现。
另外，MATLAB还提供了一些更复杂的矩阵运算函数，如矩阵的转置、逆等。
3.图像处理在MATLAB中，可以使用imread函数加载图像文件，使用imshow函数显示图像，使用imwrite函数保存图像。
另外，还可以使用一些图像处理函数对图像进行处理，如灰度化、二值化、平滑滤波等。
四、实验步骤1.调试程序首先，在MATLAB的Editor界面中打开要调试的程序文件。
然后，在程序的其中一行上点击左键，即设置了一个断点。
最后，点击运行按钮，程序会在断点处停下，然后可以通过单步调试和逐行运行来逐步查看程序的执行过程和变量的取值。
2.矩阵运算首先，定义两个矩阵A和B，并赋值。
然后，使用"+"操作符对两个矩阵进行相加，得到矩阵C。
最后，使用disp函数显示矩阵C的值。
3.图像处理首先，使用imread函数加载一张图像。
然后，使用imshow函数显示加载的图像。
接着，使用rgb2gray函数将彩色图像转换为灰度图像。
MATLAB实验报告
班级：通信班姓名：彭羊平学号： 222008315222033实验一：matlab的基本操作一、实验目的：1、了解MATLAB的集成环境，熟悉其基本操作。
2、了解MATLAB的基础知识，包括矩阵的建立、简单操作、逻辑操作和关系运算。
3、熟悉基本的数学函数和逻辑函数。
4、在命令窗口输入命令完成一些简单的功能，为MATLAB程序设计奠定基础。
二、实验内容：1、利用diag等函数产生下列矩阵：然后利用reshape函数将它们变换成行向量。
2、利用rand函数产生（0,1）间均匀分布的10*10随机矩阵A，然后统计A中大于等于0.6的元素的个数。
3、有一矩阵A，找出矩阵中值大于1的元素，并将它们重新排列成列向量B。
4、在一测量矩阵A(6*6)，存在有奇异值（假设大于100的值为奇异值），编程实现删去奇异值所在的行。
三、实验结果：1、程序如下：a=fliplr(diag([8,-7,2])+diag([5,3],-1))a=reshape(a,1,9)b=diag([2,5,8])+diag([4],2)+diag([7],-2) b=reshape(b,1,9)结果如下：2、程序如下：A=rand(10)A=A>0.6sum(sum(A))结果如下：3、程序如下：A=fix(5*rand(5))a=A>1;K=find(a);for k=1:length(K)B(k)=A(K(k));endB=B'结果如下：4、程序如下：A=fix(100*rand(8))+5 B=A>100;k=find(any(B'))A(k,:)=0结果如下:班级：通信班姓名：彭羊平学号： 222008315222033实验二： matlab绘图（1）一、实验目的：1、了解MATLAB图形系统和各种图形函数。
2、熟悉MATLAB的基本图形操作，具备MATLAB画图能力。
3、熟悉各种数学函数，并通过图形函数画出。
matlab实验内容答案解析
实验报告说明：matlab 课程实验需撰写8个实验报告，每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。
第一次实验内容：实验一 MATLAB 运算基础一、实验目的1.熟悉启动和退出MATLAB 的方法。
2.熟悉MATLAB 命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握MATLAB 各种表达式的书写规则以及常用函数的使用。
二、实验内容1.先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
（1）22sin 8511z e ︒=+（2）12ln(2z x =，其中2120.455i +⎡⎤=⎢⎥-⎣⎦（3）0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02a ae e z a a --=+=---提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。
（4）2220141122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪--≤<⎩，其中t ＝0：0.5：2.5 提示：用逻辑表达式求分段函数值。
2.已知12344347873657A -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，131203327B -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦求下列表达式的值：（1） A+6=B 和A-B+I(其中I 为单位矩阵)。
（2） A*B 和A.*B 。
（3） A^3和A^.3 。
（4） A/B 和B\A 。
（5）[A ，B]和[A([1,3],;);B^2] 。
3.设有矩阵A 和B12345678910111213141516171819202122232425A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦， 30161769023497041311B ⎡⎤⎢⎥-⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎣⎦（1） 求它们的乘积C 。
（2） 将矩阵C 的右下角3×2子矩阵赋给D （3） 查看MATLAB 工作空间使用情况。
4.完成下列操作：（1）求[100，999]之间能被21整除的数的个数。
matlab实验报告
实验一：MATLAB基本操作一、实验目的1、学习掌握MA TLAB语言的基本操作方法2、掌握命令窗口的使用3、熟悉MA TLAB的数据表示、基本运算和程序控制语句4、熟悉MA TLAB程序设计的基本方法二、实验内容和要求2.b.帮助命令的使用，查找SQRT函数的使用方法按Start-help-search，输入sqrtc.矩阵运算（1）已知A=[1 2;3 4];B=[5 5;7 8];求A^2*B输入A=[1 2;3 4];B=[5 5;7 8];A^2*B结果为：ans =105 115229 251(2)矩阵除法已知A=[1 2 3 ;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];求A/B,A\BA/B=1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000A\B=1.0e+016 *0.3152 -1.2609 0.9457-0.6304 2.5218 -1.89130.3152 -1.2609 0.9457（3）矩阵的转置及共轭转置已知A=[15+i,2-i,1;6*i,4,9-i];求A.',A'A=15.0000 + 1.0000i 2.0000 - 1.0000i 1.00000 + 6.0000i 4.0000 9.0000 - 1.0000iA.'=15.0000 + 1.0000i 0 + 6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000iA'=15.0000 - 1.0000i 0 - 6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i（4）使用冒号选出指定元素已知A=[1 2 3;4 5 6;7 8 9];求A中第三列前两个元素；A中所有第二行的元素A =1 2 34 5 67 8 9求A中第三列前两个元素:A([1,2],3)ans =36求A中所有第二行的元素:A(2,1:1:end)ans =4 5 6（5）方括号[]用magic函数生成一个4阶魔方矩阵，删除该矩阵的第四列magic(4)ans =16 2 3 135 11 10 89 7 6 124 14 15 1a=ans；a(:,[1,2,3])ans =16 2 35 11 109 7 64 14 153..多项式（1）求多项式p(X)=x^3-2x-4的根（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]求矩阵A的特征多项式；求矩阵多项式中未知数为20时的值；把矩阵A作为未知数代入到多项式中；(1)a=[1,0,-2,-4]x=roots(a)x =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i(2)求特征多项式系数：A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]p=poly(A)p =1.0000 -6.9000 -77.2600 -86.1300 604.5500求矩阵多项式中未知数为20时的值：c=polyval(p,20)c =7.2778e+004把矩阵A作为未知数代入到多项式中：c=polyval(p,A)c =1.0e+003 *0.3801 -0.4545 -1.9951 0.4601-1.9951 0.2093 -1.9951 -2.8880-0.4545 -4.8978 0.6046 0.43530.4353 0.0840 -0.4545 -1.16174、基本程序设计（1）编写命令文件：计算1+2+3+……+n<2000；n=0;s=0;while s<2000n=n+1;s=s+n;endn=n-1n =62（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
matlab实验报告
实验仿真一：控制系统中的MATLAB编程应用已知单位负反馈的开环传递函数为：G(s)=ωξω（1）假设自然频率ω给定，取五个不同的阻尼比ξ。
试用MATLAB编程，在同一张图中绘制其单位阶跃响应曲线，并分析阻尼比对二阶动态系统的影响。
取ω=2，ξ=0.25, 0.5, 1，2，4由开环传递函数求出闭环传递函数：G(s)=ωξωω将数据带入闭环传递函数得：G(s)=，ξ则随着ξ的值变化，绘制的曲线将不同。
写入MATLAB程序如下：clearclcnum=4;den1=[1 1 4];sys1=tf(num,den1);den2=[1 2 4];sys2=tf(num,den2);den3=[1 4 4];sys3=tf(num,den3);den4=[1 8 4];sys4=tf(num,den4);den5=[1 16 4];sys5=tf(num,den5);step(sys1,'b-',sys2,'g-',sys3,'r-',sys4,'c-',sys5,'y-',20);图1-1 阻尼比不同的单位阶跃响应曲线图阻尼比对二阶动态系统的影响：随着阻尼比ξ的增大，系统超调量逐渐减小。
当阻尼比大于1时，系统响应曲线为单调曲线，已经没有振荡了。
（2）假设阻尼比ξ给定，取五个不同的自然频率ω。
试用MATLAB编程，在同一张图中绘制其单位阶跃响应曲线，并分析自然频率对二阶动态系统的影响。
取ξ=0.5，ω=0.5，1，2，3， 4将数据带入闭环传递函数得：G(s)=ω，ωω则随着ω的值变化，绘制的曲线将不同写入MATLAB程序如下：clearclcnum1=0.25;den1=[1 0.5 0.25];sys1=tf(num1,den1);num2=1;den2=[1 1 1];sys2=tf(num2,den2);num3=4;den3=[1 2 4];sys3=tf(num3,den3);num4=9;den4=[1 3 9];sys4=tf(num4,den4);num5=16;den5=[1 4 16];sys5=tf(num5,den5);step(sys1,'b-',sys2,'g-',sys3,'r-',sys4,'c-',sys5,'y-',20);图1-2 自然频率不同的单位阶跃响应曲线图自然频率对二阶动态系统的影响：随着自然频率的增加，系统的响应速度加快，而响应曲线的峰值保持不变。
matlab实验报告(实验4)
学生实验报告开课学院及实验室: 机电学院2012年12月21日学院机电学院年级、专业、班姓名学号实验课程名称MATLAB程序设计成绩实验项目名称实验4: 数据和函数的可视化指导老师一、实验目的1、掌握MATLAB绘图的基本步骤和相关指令调用的先后顺序。
2、掌握MATLAB绘图指令的调用方法。
二、实验内容数学函数从形式上可以分为离散函数和连续函数。
MATLAB对这两种函数数据的可视化都提供了相应的指令。
仔细阅读教材【例5.1-1】的实现代码, 运行并保存结果；并改用stem函数, 画出【例5.1-1】的序列图。
仔细阅读教材【例5.1-2】的实现代码, 运行并保存结果；并分别使用描点和连折线方式, 画出连续函数y=xcosx的近似图形（采样点数自定, 要求画出的图尽量接近原连续函数的图）。
仔细阅读【例5.2-2】的实现代码, 理解plot指令画多条曲线的运用方法, 运行并保存结果；并使用plot函数和legend函数, 在同一个图形窗口上画出y=sint和y=sin(2t)在[0,2pi]区间上的图形, 并标出图例。
仔细阅读【例5.2-4】的实现代码, 理解图形标识选项的运用方法, 运行并保存结果；并修改代码, 把“sin(t)”字体改为正体, 大小改为20, “极大值”改为宋体。
阅读【例5.2-6】, 理解使用hold on指令画多幅图的方法, 运行并保存结果。
阅读【例5.2-8】, 理解使用subplot函数画多个子图的方法, 运行并保存结果。
（1）综合实验: 阅读以下关于通过绘制二阶系统阶跃响应综合演示图形标识的示例, 理解示例中所有图形标识指令的作用, 掌握各个图形标识指令的运用方法, 并在原指令上改动以实现以下功能:（2）把横坐标范围改为0至5pi, 纵坐标范围改为0至2；（3）把图中的横轴的刻度改为从0开始到4pi, 中间各点间隔为pi/2；纵轴刻度改为从0开始到1.5, 中间各点间隔为0.3；（4）把图中的α改为σ。
MATLAB原理应用实验报告第三章(符号运算)
《MATLAB原理及应用》实验报告第三章MATLAB的符号运算一．实验目的1、掌握符号对象的命名方法2、掌握符号表达式的基本运算3、掌握符号级数的求法二．实验设备计算机、MATLAB软件三．实验内容1.确定符号表达式的变量为了简化符号对象的操作和计算，MATLAB为用户提过了findsym命令。
r=findsym(S)确定符号表达式或者矩阵S中自由符号变量r=findsym(S，n)确定符号表达式或者矩阵S中靠近x最近的n个独立符号变量。
【实验3-1】使用MA TLAB的命令确定符号表达式的变量。
在MATLAB的命令窗口中输入下例内容：>> syms a x y z t确定下面简单符号表达式中的符号变量信息：>>findsym(sin(pi*t))ans =t确定下面简单符号表达式中的符号变量信息：>>findsym(x+i*y-j*z)ans =x, y, z确定下面简单符号表达式中的符号变量信息：>>findsym(a+y,1)ans =y2.符号表达式元算1.符号表达式的四则运算表达式的四则运算与数字运算一样，用+、-、/、运算符实现，其运算结果依然是一个符号表达式。
【实验3-2】在MATLAB的命令窗口中输入下例内容：>>f=sym('2*x^2+3*x-5');%定义符号表达式g=sym('x^2-x+7');f+gans =3*x^2+2*x+2ans =3*x^2+2*x+2>> f^gans =(2*x^2+3*x-5)^(x^2-x+7)3.符号表达式的提取分子和分母运算如果符号表达式是一个有理分式或可以展开为有理分式，可以可利用numden函数来提取符号表达式的分子或分母。
期一般调用格式为[n,d]=numden函数来提取符号表达式该函数提取的符号表达式s的分子和分母，分别将它们存放在n和d中。
matlab 实验报告
matlab 实验报告Matlab实验报告引言：Matlab是一种强大的数值计算和可视化软件，广泛应用于科学、工程和经济等领域。
本实验报告将介绍我在使用Matlab进行实验过程中的一些经验和结果。
实验一：矩阵运算在这个实验中，我使用Matlab进行了矩阵运算。
首先，我创建了一个3x3的矩阵A和一个3x1的矩阵B，并进行了矩阵相乘运算。
通过Matlab的矩阵乘法运算符*，我得到了一个3x1的结果矩阵C。
接着，我对矩阵C进行了转置操作，得到了一个1x3的矩阵D。
最后，我计算了矩阵C和矩阵D的点积，并将结果输出。
实验二：数据可视化在这个实验中，我使用Matlab进行了数据可视化。
我选择了一组实验数据，包括时间和温度两个变量。
首先，我将数据存储在一个矩阵中，并使用Matlab的plot函数将时间和温度之间的关系绘制成曲线图。
接着，我使用Matlab的xlabel、ylabel和title函数添加了横轴、纵轴和标题。
最后，我使用Matlab的legend函数添加了图例，以便更好地理解图表。
实验三：数值积分在这个实验中，我使用Matlab进行了数值积分。
我选择了一个函数f(x)进行积分计算。
首先，我使用Matlab的syms函数定义了符号变量x，并定义了函数f(x)。
接着，我使用Matlab的int函数对函数f(x)进行积分计算，并将结果输出。
为了验证结果的准确性，我还使用了Matlab的diff函数对积分结果进行了求导操作，并与原函数f(x)进行了比较。
实验四：信号处理在这个实验中，我使用Matlab进行了信号处理。
我选择了一个音频文件，并使用Matlab的audioread函数读取了该文件。
接着，我使用Matlab的fft函数对音频信号进行了傅里叶变换，并将结果绘制成频谱图。
为了进一步分析信号的特征，我还使用了Matlab的spectrogram函数绘制了信号的时频图。
通过对信号的频谱和时频图的观察，我可以更好地理解信号的频率和时域特性。
Matlab实验报告
实验一：Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境。
2．学习使用图形函数计算器命令funtool及其环境。
二、实验容熟悉Matlab操作环境，认识命令窗口、存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1)y=sin(x); (2) y=(1+x)^3*(2-x);求下列函数的符号积分(1)y=cos(x);(2)y=1/(1+x^2);(3)y=1/sqrt(1-x^2);(4)y=(x1)/(x+1)/(x+2)求反函数(1)y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2));代数式的化简(1)(x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2)sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6)y6=x^2/23．两个函数之间的操作求和(1)s in(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5乘积(1)e xp(-x)*sin(x) (2) sin(x)*x商(1)s in(x)/cos(x); (2) x/(1+x^2); (3) 1/(x-1)/(x-2);求复合函数(1)y=exp(u) u=sin(x) (2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x) (4) y=sinh(u) u=-x实验二：MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。
matlab实验一实验报告
matlab实验一实验报告实验一：Matlab实验报告引言：Matlab是一种强大的数学软件工具，广泛应用于科学计算、数据分析和工程设计等领域。
本实验旨在通过使用Matlab解决实际问题，探索其功能和应用。
一、实验目的本次实验的主要目的是熟悉Matlab的基本操作和常用函数，了解其在科学计算中的应用。
二、实验内容1. 数值计算在Matlab中，我们可以进行各种数值计算，包括基本的加减乘除运算，以及更复杂的矩阵运算和方程求解。
通过编写相应的代码，我们可以实现这些功能。
例如，我们可以使用Matlab计算两个矩阵的乘积，并输出结果。
代码如下：```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A * B;disp(C);```2. 数据可视化Matlab还提供了强大的数据可视化功能，可以将数据以图表的形式展示出来，更直观地观察数据的规律和趋势。
例如，我们可以使用Matlab绘制一个简单的折线图，来展示某个物体在不同时间下的位置变化。
代码如下：```matlabt = 0:0.1:10;x = sin(t);plot(t, x);xlabel('Time');ylabel('Position');title('Position vs. Time');```3. 图像处理Matlab还可以进行图像处理，包括图像的读取、处理和保存等操作。
我们可以通过Matlab对图像进行增强、滤波、分割等处理，以及进行图像的压缩和重建。
例如，我们可以使用Matlab读取一张图片，并对其进行灰度化处理。
代码如下：```matlabimg = imread('image.jpg');gray_img = rgb2gray(img);imshow(gray_img);```三、实验结果与分析在本次实验中，我们成功完成了数值计算、数据可视化和图像处理等任务。
matlab实验报告实验二
matlab实验报告实验二Matlab实验报告实验二引言Matlab是一种功能强大的数学软件，广泛应用于科学研究和工程实践中。
在实验二中，我们将探索Matlab的图像处理功能，并通过实际案例来展示其应用。
图像处理基础图像处理是指对图像进行数字化处理的过程，其目的是改善图像质量、提取有用信息或实现特定的应用需求。
在Matlab中，我们可以利用各种函数和工具箱来实现图像处理的各种任务，如图像增强、滤波、分割和特征提取等。
实验步骤1. 图像读取与显示在Matlab中，我们可以使用imread函数读取图像文件，并使用imshow函数将图像显示在屏幕上。
例如，我们可以读取一张名为"lena.jpg"的图像，并显示出来：```matlabimg = imread('lena.jpg');imshow(img);```2. 图像灰度化图像灰度化是将彩色图像转换为灰度图像的过程。
在Matlab中，我们可以使用rgb2gray函数将彩色图像转换为灰度图像。
例如，我们可以将上一步读取的图像转换为灰度图像：```matlabgray_img = rgb2gray(img);imshow(gray_img);```3. 图像二值化图像二值化是将灰度图像转换为二值图像的过程，其中只包含黑色和白色两种颜色。
在Matlab中，我们可以使用imbinarize函数将灰度图像二值化。
例如，我们可以将上一步得到的灰度图像二值化：```matlabbinary_img = imbinarize(gray_img);imshow(binary_img);```4. 图像平滑图像平滑是指去除图像中的噪声或细节，使得图像更加平滑和清晰。
在Matlab 中，我们可以使用imfilter函数对图像进行平滑处理。
例如，我们可以对上一步得到的二值图像进行平滑处理：```matlabsmooth_img = imfilter(binary_img, fspecial('average'));imshow(smooth_img);```5. 图像边缘检测图像边缘检测是指提取图像中物体边缘的过程，常用于目标检测和图像分割。
MATLAB实验报告
实验一基本操作和简单语句输入一、实验环境计算机MATLAB软件二、实验目的1、熟悉MATLAB的命令窗口。
2、掌握MATLAB的一些基本操作，能够进行一般的数值计算。
3、实现语句的重调和修改。
三、实验内容与步骤1、启动MATLAB2、观察MATLAB窗口的各个组成部分（1）了解菜单栏各窗口项的功能，用鼠标打开MA TLAB的各个菜单，在状态栏里显示当前鼠标所指的菜单项的含义。
（2）用鼠标指向常用工具栏的每个工具按钮，了解各工具按钮的含义。
3、命令窗口的打开和关闭（1）查看菜单窗口中有哪些菜单命令。
（2）在命令窗口中输入命令a=3；b=4；y=a*b+a/b，然后回车，查看命令显示结果。
>> a=3;b=4;y=a*b+a/b,y =12.7500（3）利用MATLAB中编辑命令时常用的按键功能，调用上一个语句，对它进行修改（如把分号改成逗号，看运行结果），并把运行结果复制到word中保存。
>> a=3,b=4,y=a*b+a/b,a =3b =4y =12.7500（4）关闭命令窗口。
（5）打开命令窗口。
4、使用MATLAB帮助熟悉MATLAB的帮助系统，通过帮助系统了解有关内容。
5、在命令窗口中输入demo，将出现MA TLAB的演示窗，通过演示窗，对MATLAB 的功能进行一次浏览。
四、练习1、计算y=x^3+(x-0.98)^2/(x+1.25)^3-5*(x+1/x)，x=2，x=3时的值。
>> x=2;y=x^3+(x-0.98)^2/(x+1.25)^3-5*(x+1/x)y =-4.4697>> x=3;y=x^3+(x-0.98)^2/(x+1.25)^3-5*(x+1/x)y =10.38652、计算cos(pi/3)-sqrt(9-sqrt(2))>> cos(pi/3)-sqrt(9-sqrt(2))ans =-2.25423、已知：a=3,A=4,b=a^2,B=b^2-1,c=a+A-2*B,C=a+B+2*c，求：C>> a=3;A=4;b=a^2;B=b^2-1;c=a+A-2*B;C=a+B+2*cC =-2234、复数z1=1+3*i,z2=1+2*i,z3=2*exp((pi/6)*i),求z=(z1*z2)/z3>> z1=1+3*i;z2=1+2*i;z3=2*exp((pi/6)*i);z=(z1*z2)/z3z =-0.9151 + 3.4151i实验二矩阵和数组的操作一、实验环境计算机MATLAB软件二、实验目的1、掌握矩阵和数组的一般操作，包括创建、保存、修改和调用等。
Matlab实验报告
实验一：MATLAB语言平台与SIMULINK工具箱一、实验目的：1．学习了解MATLAB语言环境；2．练习MATLAB命令；4．掌握SIMULINK工具箱的应用。
二、实验要求：1．实验内容（1）学习了解MATLAB语言环境；（2）练习MATLAB命令的基本操作；（3）掌握SIMULINK仿真环境的使用方法;(4) 掌握线性系统仿真常用基本模块的用法2．要求按照实验文档的要求与步骤完成实验，撰写实验报告。
三、实验基本知识：1.熟悉MATLAB环境：MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。
2.掌握MATLAB常用命令3.了解SIMULINK模块库中各子模块基本功能四、实验内容：1.学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）help eyeEYE Identity matrix.EYE(N) is the N-by-N identity matrix.EYE(M,N) or EYE([M,N]) is an M-by-N matrix with 1's onthe diagonal and zeros elsewhere.EYE(SIZE(A)) is the same size as A.EYE with no arguments is the scalar 1.EYE(M,N,CLASSNAME) or EYE([M,N],CLASSNAME) is an M-by-N matrix with 1'sof class CLASSNAME on the diagonal and zeros elsewhere.Example:x = eye(2,3,'int8');See also speye, ones, zeros, rand, randn.Reference page in Help browserdoc eyeA=eye(3)A=1 0 00 1 00 0 12．任意创建3个变量，分别为数值、矩阵、字符串,保存到mydat.mat文件，再把my_data.mat文件中的变量读取到MatLab工作空间内；查找已创建变量的信息，删除变量；学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
matlab实验报告
matlab实验报告实验名称熟悉matlab环境课程名称数值分析学生姓名乔幸帅学号 331201000018专业信号与信息处理一、实验目的本次实验的目的是在理论基础的学习前提下熟悉matlab的基本操作，学会matlab的基本绘图以及函数的编写，为学习数值分析这门课奠定基础，并为以后在其他学习和工程使用中打下基础。
二、实验内容1.熟悉matlab的三种指令输入方式2.学会简单的绘图命令3.用matlab编程学会创建函数三、实验步骤3.1 matlab的三种指令输入方式3.1.1 在指令框中直接输入指令并查看结果，如图所示：3.1.2 创建M文件写入指令并运行○1点击matlab的FILE按键选中NEW-Script创建一个如图所示：○2在M文件中输入想要运行的程序并保存○3点击debug中的Run按键运行程序显示结果如图：012345678910matlab Notebook的启动点击matlab左下角的start按钮，点击里面的notebook选项即可新建一个类似word的文档，在里面输入程序，并在Notebook下拉菜单选择Define Input Cell，继续选择Ealuate Input Cell，即得结果如图所示3．2 matlab的简单绘图绘制图像的最基本函数是plot（），其基本调用格式是plot（x,y）如果要绘制多个图像可以使用subplot函数，即在同一个窗口绘制多个子图。
用subplot函数花多个子图。
t=0:pi/100:2*pi;y1=sin(t);y2=cos(t);y3=sin(t)+cos(t);y4=sin(t).*cos(t);subplot(2,2,1)plot(t,y1)xlabel('\itt\rm/s');ylabel('y_1=sin(t)');subplot(2,2,2)plot(t,y2)xlabel('\itt\rm/s'); ylabel('y_2=cos(t)'); subplot(2,2,3) plot(t,y3) gridxlabel('\itt\rm/s'); ylabel('y_3=sin(t)+cos(t)'); subplot(2,2,4) plot(t,y4) gridxlabel('\itt\rm/s'); ylabel('y_4=sin(t)*cos(t)'); 程序运行结果如图所示2468-1-0.500.51t /sy 1=s i n (t )02468-1-0.500.51t /sy 2=c o s (t )2468-2-1012t /sy 3=s i n (t )+c o s (t )02468-0.50.5t /sy 4=s i n (t )*c o s (t )3．3 用 matlab 编写函数MATLAB 的 M 函数是由 function 语句引导的，其基本格式如下： function [返回变量列表] = 函数名 (输入变量列表)[例子]：假设我们想生成一个 n*m 阶的 Hilbert 矩阵, 它的第 i 行第 j 列的元素值为 1/(i+j-1)。
MATLAB实验报告(猎狗追兔子的问题)
MATLAB实验报告电气14班张程21104011202012年4月12日星期四一．实验目的1.学会用MATLAB软件求解微分方程的初值问题。
2.学会根据实际问题建立简单微分方程数学模型。
3.了解级计算机数据仿真、数据模拟的基本方法。
二．实验题目有一只猎狗在B处发现了一只兔子在正东北方距离它200米的地方O处，此时兔子开始以8米每秒的速度向正西北方向距离为120米的洞口A全速跑去，假设猎狗在追赶兔子时始终朝着兔子的方向全速奔跑。
（1）问猎狗能追上兔子的最小速度是多少？（2）选取猎狗的速度分别为15、18米每秒，计算猎狗追上兔子是所跑过的路程和所用的时间。
（3）画出猎狗追赶兔子奔跑的曲线图。
三．实验过程（1）将所有路径转化入第一象限，从而转化为教材中缉私艇和走私船的问题。
兔子起始位置为（0，0），奔跑方向为Y轴正方向，猎狗起始位置为（200，0）。
兔子刚好被追上时跑的距离Y=CR/1-R²。
（R为兔子与猎狗速度之比）时间T=CR/A(1-R²)=BC/(B²-A²)。
当兔子进洞的时候刚好被追到，这种情况下猎狗所需速度最小。
将C=200,Y=120带入方程并利用MATLAB求解：所以最小速度为17.08m/s。
(2)程序如下：运行结果：t =18.2000 s =273.0000将程序中b改为18运行结果为：t =13.6000 s =244.8000 兔子跑过距离分别为145.6000、108.8000，所以第一次兔子已经进洞，猎狗追不上。
（3）下图分别为猎狗速度为18m/s、15m/s时的模拟图。
四．反思总结1.实际问题的解决可以通过matlab等工具进行很好地计算和模拟，关键在于找到适当的方法。
2.实验过程中总会有这样或者那样的错误，在不断的更正中发现对程序的了解更深了一步。
3.同一个问题可以有多个解答方法，只要擅于分析和发现，在不断的总结中就可以有所收获。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告10姓名： 学号： 时间：2012年5月5日实验名称：用MATLAB 进行有约束优化问题的计算机求解； 实验目的：学会如何运用MATLAB 进行优化问题求解； 实验内容：1、试求解线性规划问题 min -3x 1+4 x 2-2x 3+5x 4⎪⎪⎩⎪⎪⎨⎧-≥-≥---≤+-+-=-+-无约束43,2,1432143214321,1232142224..x x x x x x x x x x x x x x t s x 2、某车间有甲、乙两台机床，可用于加工三种工件。
假定这两台车床的可用台时数分别为600和700，三种工件的数量分别为300、500和400，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表。
问怎样分配车1、求解如下f=[-3 4 -2 5]';A=[1 1 -1 2;-2 3 1 1];B=[14; 2];Ae=[4 -1 2 -1];Be=[-2];xm=[-1;-1;-1; ]; ff=optimset;ff.TOLX=1e-15;ff.TolFun=1e-20; TolFun=1e-20;[x,f_opt,key,c]=linprog(f,A,B,Ae,Be,xm,[],[],ff)Exiting: Maximum number of iterations exceeded; increase options.MaxIter.x =-1 5/2 -1 -13/2f_opt =-35/2key =0 c =iterations: 85algorithm: 'large-scale: interior point' cgiterations: 0message: [1x73 char]2、解：设甲车床加工工件1、工件2、工件3的数量分别为321,,x x x ，乙车床上加工工件1、工件2、工件3的数量为分别为654,,x x x ，则：6543218121110913min x x x x x x +++++ ⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=≥≤++≤++=+=+=+6,...,2,1,07003.12.15.06001.14.0400500300..654321635241i x x x x x x x x x x x x x t s x i解答：f=[13 9 10 11 12 8]';A=[0.4 1.1 1 0 0 0;0 0 0 0.5 1.2 1.3]; B=[600;700];Ae=[1 0 0 1 0 0;0 1 0 0 1 0;0 0 1 0 0 1]; Be=[300 500 400]; xm=[0 0 0 0 0 0]; ff.TolX=1e-15; ff.TolFun=1e-20; ff.TolCon=1e-20;[x,f_opt,key,c]=linprog(f,A,B,Ae,Be,xm,[],[],ff)Optimization terminated. x =1/664613997892457940000000000000000000 5001/225184372151181690000000 3001/11100951164476451000000000000 400f_opt =11000key =1 c =iterations: 7algorithm: 'large-scale: interior point' cgiterations: 0message: 'Optimization terminated.'因此，甲车床加工工件1、工件2、工件3的数量分别为0，500，0，乙车床上加工工件1、工件2、工件3的数量为分别为300，0，400，此时加工费用最低，为11000实验报告11姓名： 学号： 时间：2012年5月15日实验名称：二次型规划和非线性规划的最优化问题的计算机求解实验目的：使用MATLAB 软件解决二次型规划和非线性规划最优化问题 实验内容：P211~21211、试求解下面的二次型规划问题。
2122212136442min x x x x x x --+-x ⎪⎩⎪⎨⎧≥≤+≤+0943..2,12121x x x x x t s求解过程：首先，原始问题的目标函数为：2122212136442)(x x x x x x x f --+-=]364[],80;84[--=-=T f H程序如下：f=[-6 -3];H=[4 -5;-3 8];OPT=optimset;rgeScale='offo';A=[1,1;4,1];B=[3;9];Ae=[];Be=[];LB=zeros(2,1); [x,f_opt]=quadprog(H,f,A,B,Ae,Be,LB,[],[],OPT) 运行结果如下: x =1.9500 1.0500f_opt =-11.0250故，当05.1,95.121==x x 时，0250.11min -=f 。
12、试求解下面的非线性规划问题。
)12424(min 22122211++++x x x x x e xx ⎪⎪⎩⎪⎪⎨⎧≤≤--≥≥++-≤+10,10105.10..2121212121x x x x x x x x x x t s 程序如下：1、function y=opt_fun1(x)y=exp(x(1))*(4*x(1)*x(1)+2*x(2)*x(2)+4*x(1)*x(2)-2*x(2)+1);2、function [c,ceq]=opt_con1(x)ceq=[];c = [x(1)+x(2);x(1)*x(2)-x(1)-x(2)+1.5;-x(1)*x(2)-10];运行：ff=optimset; rgeScale='off'; ff.Display='iter'; ff.TolFun=1e-30; ff.TolX=1e-15; ff.TolCon=1e-20;x0=[1;1]; xm=[-10;-10]; xM=[10;10]; A=[]; B=[]; Aeq=[]; Beq=[];[x,f_opt,c,d]=fmincon('opt_fun1',x0,A,B,Aeq,Beq,xm,xM,'op t_con1',ff);x,f_opt,d.funcCount 运行结果如下： x =1.1865-1.7104f_opt =25.5022ans =201故，当7104.1,1865.121-==x x 时，5022.25min =y 。
15、试求解下面的0~1规划问题。
5432131075min x x x x x ++++x ⎪⎪⎩⎪⎪⎨⎧≤≤≤--+-≥+--+-≥-++-10122022362245.54325432154321i x x x x x x x x x x x x x x x t s 程序如下：f=[5 7 10 3 1]; A=[-1 1 -5 -1 4; 2 -6 3 2 -2; 0 -2 2 -1 -1]; B=[-2;0;1];xm=[0;0;0]; xM=[1;1;1]; x=bintprog(f,A,B,[],[])' 运行结果如下： x =0 1 1 0 0实验报告 12姓名： 学号： 时间：2012年5月17日实验名称：线性代数问题的计算机求解 实验目的：运用MATLAB 求解线性代数问题实验内容：18.抽查某地区3所小学五年级男学生的身高由下表给出，问该地区这3所小学五年级男学1.操作工之间的差异是否显著？2．机器之间的差异是否显著？3.交互作用是否显著（a=0.05）?18. A=[128.1 134.1 133.1 138.9 140.8 127.4;150.3 147.9 136.8 126 150.7 155.8;140.6 143.1 144.5 143.7 148.5 146.4]A =128.1000 134.1000 133.1000 138.9000 140.8000 127.4000150.3000 147.9000 136.8000 126.0000 150.7000 155.8000140.6000 143.1000 144.5000 143.7000 148.5000 146.4000>> m=mean(A),[p,tb1,stats]=anova1(A)m =139.6667 141.7000 138.1333 136.2000 146.6667 143.2000p =0.7838tb1 ='Source' 'SS' 'df' 'MS' 'F' 'Prob>F''Columns' [ 211.3361] [ 5] [42.2672] [0.4813] [0.7838]'Error' [1.0538e+003] [12] [87.8167] [] []'Total' [1.2651e+003] [17] [] [] []stats =gnames: [6x1 char]n: [3 3 3 3 3 3]source: 'anova1'means: [139.6667 141.7000 138.1333 136.2000 146.6667 143.2000] df: 12s: 9.3711可知道，P=0.7838>a=0.05,故应该接受假设，身高没有显著影响。
19.>>A=[15,15,17,17,17,17,15,17,16,18,20,22;19,19,16,15,15,15,18,17,16,15,16,17;16,18,21,19,22,22, 18,18,18,17,17,17];>> m=mean(A),[p,tab,stats]=anova1(A)m =16.6667 17.3333 18.0000 17.0000 18.0000 18.0000 17.0000 17.3333 16.6667 16.6667 17.6667 18.6667p =0.9920tab ='Source' 'SS' 'df' 'MS' 'F' 'Prob>F''Columns' [ 14.0833] [11] [1.2803] [0.2352] [0.9920]'Error' [130.6667] [24] [5.4444] [] [] 'Total' [144.7500] [35] [] [] []stats =gnames: [12x2 char]n: [3 3 3 3 3 3 3 3 3 3 3 3] source: 'anova1'means: [16.6667 17.3333 18 17 18 18 17 17.3333 16.6667 16.6667 17.6667 18.6667] df: 24 s: 2.3333由于得出的概率值α>=9920.0p ，其中α=0.05，所以操作工之间的差异不显著。