北科,北京科技大学,数学实验,MATLAB第二次作业
数学实验练习二 参考答案
数学实验二实验内容：学习matlab的m文件编写和函数的编写，体会matlab编程特点，掌握matlab 的编程基本方法。
要求：一．学习ppt教案的例题代码，能正确的输入、运行代码；二．写出如下各段代码的作用，将以下各段循环执行的代码，改为不需要循环的矩阵和数组运行，并使用tic,toc测试不同代码的执行时间：%程序1,文件名：ex2_2_1.mticdx = pi/30;nx = 1 + 2*pi/dx;for i = 1:nxx(i) = (i-1)*dx;y(i) = sin(3*x(i));endtoc以上程序实现将[0,2*pi]间隔pi/30分成60等分，x和y分别为61个元素的数组，y为计算sin(3x)的值。
以上程序可以使用简单的matlab数组计算实现：x2=0:pi/30:2*pi;y2=sin(3*x2);大家可以比较一下，x1和x2完全相同，y和y2也完全相同。
%程序2，文件名：ex2_2_2.mticA=round(2+rand(50,60)*6); 生成一个在[2,8]上均匀分布的50*60随机数组[X,Y]=size(A); 求出其大小；X=50,Y=60minA=A(1,1); 设最小值为矩阵A的第1行1列的元素for i=1:Xfor j=1:Yif A(i,j)<minAminA=A(i,j);minX=i;minY=j;endendend 以上程序按行、列搜索矩阵A的最小值，若当前值A(I,j)小，则将最小值设为当前值；[minA ,minX,minY] 输出矩阵最小值minA及矩阵最小值所在的行minX、列minY。
toc上述程序可以使用find函数及min函数实现；此时只需： minA=min(A(:));[minX,minY]=find(A ’==minA,1);%注意此处需将矩阵A 转置，因为matlab 中是按列优先搜索的，而题目的程序是按行有限搜索。
Matlab第二次实验报告
Matlab第二次实验报告实验目的：1、了解plot函数和subplot函数的基本用法和matlab绘图的基本原理。
2、了解图形的属性设置。
比如画图的颜色，画图采用的线性标识符等。
二：实验基本知识1：1. 单窗口单曲线绘图；2. 单窗口多曲线绘图；3. 单窗口多曲线分图绘图；4. 多窗口绘图；5.可任意设置颜色与线型；6.图形加注功能；7.fplot——绘制函数图函数；8.ezplot——符号函数的简易绘图函数9：subplot—将画图区域分块函数。
实验内容：将高数课本后的18个图画出上机练习程序第1——6图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1);x=-50:1:50;y=a.*x.^3;subplot(3,2,1);plot(x,y,'r');title('y=a*x^3');xlabel('X');ylabel('Y');x=0:1:50;y=sqrt(a.*x.^3);subplot(3,2,2);plot(x,y,x,-y,'r');title('y^2=a*x^3');xlabel('X');ylabel('Y');x=-3:0.1:3;y=exp(-x.^2);subplot(3,2,3);plot(x,y,'r');title('y=e^x^2');xlabel('X');ylabel('Y');x=-3:0.1:3;y=8*a^3./(x.^2+4*a^2);subplot(3,2,4);plot(x,y,'r');title('y=8*a^3./(x.^2+4*a^2)'); xlabel('X');ylabel('Y');x=0:0.01:5;y=sqrt(x.^3./(a*2.-x)); subplot(3,2,5);plot(x,y,x,-y,'r');title('y^2*(2a-x)=x^3'); xlabel('X');ylabel('Y');%t=-1:0.01:5;%x=3*a.*t./(1+t.^3);%y=3*a*t.^2./(1+t.^3);subplot(3,2,6);ezplot('x.^3+y.^3-3*3*x.*y'); %plot(x,y);title('x^3+y^3-3axy=0'); xlabel('X');ylabel('Y');输入a=3：显示结果第7——10图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0:0.1:2*pi;x=a*cos(t).^3;y=a*sin(t).^3;subplot(2,2,1);plot(x,y,'b')title('x^1.5+y^1.5=a^1.5'); xlabel('X');ylabel('Y');t=-2*pi:0.1:2*pi;x=a.*(t-sin(t));y=a.*(1-cos(t));subplot(2,2,2);plot(x,y,'b')title('°ÚÏß');xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*(1-cos(t));subplot(2,2,3);plot(p.*cos(t),p.*sin(t),'b'); title('p=a(1-cos£¨t£©)'); xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*t;subplot(2,2,4);plot(p.*cos(t),p.*sin(t),'b'); title('p=at');xlabel('X');ylabel('Y');输入a=1，第11题图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º') figure(1);t=-2*pi:pi/100:2*pi;p=exp(a.*t);plot(p.*cos(t),p.*sin(t)); title('p=e^at');xlabel('X');ylabel('Y');输入a=0.1，显示结果：clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0.2*pi:pi/100:100*pi;p=a./t;plot(p.*cos(t),p.*sin(t));title('p=at');xlabel('X');ylabel('Y');输入a=1，显示结果：第13——14题图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º');figure(1);t=0:pi/100:2*pi;p=(a^2.*sin(2*t)).^0.5;subplot(1,2,1);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*sin(2t)');xlabel('X');p=(a^2.*cos(2*t)).^0.5;subplot(1,2,2);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*cos(2t)');xlabel('X');ylabel('Y');输入a=1，显示结果：第15-16题图：clear;clc;a=input('ÇëÊäÈëa:');figure(1)t=0:pi/200:pi;p=a.*cos(3*t);subplot(1,2,1);plot(p.*cos(t),p.*sin(t),'r--');title('p=a*cos(3t)');xlabel('X');p=a.*sin(3*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'r--'); title('p=a*sin(3t)');xlabel('X');ylabel('Y');输入a=1，显示结果：第17——18题图：clear;clc;a=input('ÇëÊäÈëa:');figure(1);t=-pi:pi/200:pi;subplot(1,2,1);p=a.*sin(2*t);plot(p.*cos(t),p.*sin(t),'k-.'); title('p=a*sin(2t)');xlabel('X');ylabel('Y');p=a.*cos(2*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'k-.');title('p=a*cos(2t)');xlabel('X');ylabel('Y');输入a=1，显示结果：实验心得：这次实验主要是学习matlab软件的二维绘图功能：matlab软件还有丰富的图形修饰功能，如改变线条的形式和颜色；除此之外其还可以以多种形式如在一个图形中显示多个函数图形，一个figure中显示多个坐标系，还有一个程序中就可以用多个figure一起来表示不同的函数……在画图过程中你可能因为角度的区间取值不一样，画出来的图形可能跟书本上给出的图形有很大的差异，但没关系，只要你思路是对的，区间可以自己慢慢改动知道跟书本上给出的图形想接近。
MATLAB第二次上机实验报告
电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析学生姓名：学号：指导教师：一、实验名称实验二 线性方程组求解和函数的数值逼近二、实验目的通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。
实验涉及的核心知识点：病态方程求解、矩阵分解和方程组求解、Lagrange 插值。
实验重点与难点：算法设计和MATLAB 编程三、实验内容1. 对高阶多项式()()()()()2011220k p x x x x x k ==---=-∏编程求下面方程的解()190p x x ε+=并绘图演示方程的解与扰动量ε的关系。
2. 对220n =，生成对应的Hilbert 矩阵，计算矩阵的条件数；通过先确定解获得常向量b 的方法，确定方程组()n H x b =最后，用矩阵分解方法求解方程组，并分析计算结果。
3. 对函数()21125f x x =+ []1,1x ∈-的Chebyshev 点()()21cos 21k k x n π⎛⎫-= ⎪ ⎪+⎝⎭，1,2,,1k n =+编程进行Lagrange 插值，并分析插值结果。
四、实验数据及结果分析1. 对高阶多项式()()()()()2011220k p x x x x x k ==---=-∏编程求下面方程的解()190p x x ε+=并绘图演示方程的解与扰动量ε的关系。
p=[1,-1]; for i=2:20 n=[1,-i];p=conv(p,n); % 求多项式乘积 endm=zeros(1,21); % m 的最高次幂为20，有21项 hold on x=1:20;d=[-1,0,0.1,0.5,1]; for i=1:5delt=d(i); m(2)=delt;y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); endtitle('方程p(x)=0的解与扰动量delt 的关系')legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1')24681012141618200102030405060方程p(x)=0的解与扰动量delt 的关系delt=-1delt=0delt=0.1delt=0.5delt=12.对220n =，生成对应的Hilbert 矩阵，计算矩阵的条件数；通过先确定解获得常向量b 的方法，确定方程组()n H x b =最后，用矩阵分解方法求解方程组，并分析计算结果。
北京科技大学matlab作业1
《数学实验》报告
2016年3月
一、【实验目的】
1、熟悉matlab的基本命令和基本操作。
2、掌握matlab命令文件和函数文件的建立、运行和保存
二、【实验任务】
P16第4题
编写函数文件，计算，并求出当k=20时表达式的值。
P27第2题
矩阵A=，B=，计算A*B,A.*B,并比较两者的区别。
三、【实验程序】
P16第4题
function zr01=jiecheng(n)
zr01=0;
for t=1:n
zr01=zr01+prod(1:t);
end
jiecheng(20)
P27第2题
>>A=[123;456;789];
>>B=[468;556;322];
>>A*B
>>A.*B
四、【实验结果】
P16第4题
P27第2题
两种运算方式的区别：A*B是按矩阵的运算法则进行矩阵的乘法计算, A.*B是
两个矩阵中对应的元素相乘
五、【实验总结】
基本熟悉了matlab的基本命令和基本操作。
掌握了matlab命令文件和函数文件的建立、运行和保存。
学会了利用matlab编写简单的程序及使用。
北京科技大学matlab大作业
分析未来五年世界汽油价格问题机械1011实验目的通过本次实验，目的在于提高学生根据实际问题建立函数拟合与插值方面的问题的运用巩固我们的有关知识，在21世纪这个对能源问题而苦恼的世纪中，让我们认识到现实，为资源节约尽一份力。
实验问题根据以往油价，推测未来五年的油价走势，最终估计按此发展下去，在2012-2015年的原油价格会达到多少，我们应如何应对这些情况。
日期国际原油价(美元/桶)2003 40.332004 49.162005 59.842006 69.652007 93.512008 89.402009 68.58预计2012-2015年世界原油价格。
建立模型世界油价的走势一定可以用一个函数来拟合它，因为我们已经知道最近几年的油价走势，和2010,2011年的油价情况，所以，我们可以用拟合的函数图象与2010,2011年的原油价进行比对，找出最相近的拟合函数，对他进行分析找出其未来2012-2015年所对应的原油价格，并对他们的图形进行合理分析，看出偏差的原因。
模型求解首先我们要求出在用不同拟合曲线中最接近2010，与2011年的曲线，然后取出它单独作为分析对象。
（因世界油价取决因素过多，所以我们只能在短时间范围内进行拟合，所以我只取到2015年的世界原油价格）用matlab求数值拟合函数图象x=2003:1:2009;y=[40.33 ,49.16 , 59.84 , 69.65 ,93.51 ,89.40, 95.18 ];m2=polyfit(x,y,2);m3=polyfit(x,y,3);m5=polyfit(x,y,5);m9=polyfit(x,y,9);disp('二次拟合'),f2=poly2str(m2, 'x')disp('三次拟合'),f3=poly2str(m3,'x')disp('五次拟合'),f5=poly2str(m5, 'x')disp('九次拟合'),f9=poly2str(m9, 'x')x1=2003:1:2015;y2=polyval(m2,x1);y3=polyval(m3,x1);y5=polyval(m5,x1);y9=polyval(m9,x1);wch2=abs(y-polyval(m2,x))./ywch3=abs(y-polyval(m3,x))./ywch5=abs(y-polyval(m5,x))./ywch9=abs(y-polyval(m9,x))./ypjwch2=mean(wch2)pjwch3=mean(wch3)pjwch5=mean(wch5)pjwch9=mean(wch9)plot(x,y, 'rp',x1,y2, '--',x1,y3,'k-', x1,y5,x1,y9) legend('二次拟合'，'三次拟合'，'五次拟合'，'九次拟合')由图可以得出应选择二次拟合图象，再把二次拟合图象单独分析x=2003:1:2009;y=[40.33 ,49.16 , 59.84 , 69.65 ,93.51 ,89.40, 95.18 ];m=polyfit(x,y,3);disp('三次拟合'),f=poly2str(m,'x')wch3=abs(y-polyval(m,x))./y得到的函数m输出为m=-1.2683.*x.^3+7630.5197.*x.^2-153 02171.1633.*x+10228931022.1912由此可拟合估计2012-2015的值x=2003:2015;m=-1.2683.*x.^3+7630.5197.*x.^2-15302171.1633.*x+10228931022.1912 ;xi=2012:1:2015;mi1=interp1(x,m,xi,'*nearest');mi2=interp1(x,m,xi,'*linear');mi3=interp1(x,m,xi,'*spline');mi4=interp1(x,m,xi,'*cubic');plot(x,m,'ro',xi,mi1,'--',xi,mi2,'-',xi,mi3, 'k.-',xi,mi4,'m:')legend('原始数据', '最近点插值', '线性插值', '样条插值', '立方插值')油价已成为现今社会人们非常关注的一个话题，我们对油价的估计也许并非很准确，但也能说明一些问题，如下面专家所分析的那样。
matlab第二次训练题目
BLOCK1=idct2(BLOCK1);
BLOCK2=idct2(BLOCK2);
a=BLOCK2(1,1)/BLOCK1(1,1)-1;
ifa<0
W(p,q)=0;
else
W(p,q)=1;
end
end
end
%显示提取的水印
subplot(3,2,4);
imwrite(L,'attack.jpg','jpeg','Quality',85);
J=imread('attack.jpg','jpeg');
subplot(3,2,5);imshow(J);title('压缩后的图像');
I=imread('1.jpg');
%从压缩的图像中提取水印
forp=1:N
x1=x(1:10:256);
horz=1:10:256;
bar(horz,x1);
除了显示数据的直方图统计外，还可以对图像进行均衡处理，所用到的函数是histeq，这两个函数同样只使用于二维数据，所以也要对二维和三维数据分开处理。
直方图显示和均衡后的图像分别如下图所示
11二值图像
用j=im2bw(x);来对灰度图像到二值图像的转换。转换结果为：
case '中值滤波'
case '自适应滤波'
case '平滑滤波'
end
10直方图统计
用imhist函数对图像数据进行直方图统计，
x=imhist(handles.img(:,:,1));
北京科技大学matlab大作业
《数学实验》报告实验名称MATLAB在研究物体振动方面的应用学院专业班级姓名学号2015年 1月一、【实验目的】物体振动这样一个看似简单但又包含着很多复杂计算的运动中，在人为的计算时是很难精确的实现，而通过MATLAB可以处理诸多科学中的许多问题，利用它来研究物理学中的机械振动，不仅特别方便还非常有效。
二、【实验任务】本列举振动的一些实例，用matlab语言编制计算机程序进行仿真以达到研究简谐振动以及振动的合成，振动的计算以及受迫振动。
三、【实验程序】（一）简谐振动介绍最简单和最基本的振动是简谐振动．任何复杂的振动，都可以看成为许多简谐振动的合成．1．特点质点作简谐振动的条件是：在任何时候所受到的力与质点离开平衡位置的位移成正比，其指向与位移相反，始终指向平衡位置．所受的力与位移的关系表示为（1）式中为正的常数．对于弹簧振子，就是弹簧劲度系数2．运动的微分方程及其解根据牛顿第二定律，作简谐振动的质点的微分方程写成即（2）式中。
如下面的（3）和（4）所示，是简谐振动的圆频率。
微分方程（2）的解是（3）或（4）式（7.3）也可以表为复数形式（5）但要约定取其实数部分．利用三角公式，很容易导出A ，和B，C之间的关系即（6）3．速度和加速度作简谐振动的质点，它的速度和加速度很容易得到．只要将（7.3）对时间分别求导一次和求导两次即可，（7）（8）式（1）、（2）、（3）、（4）、（5）都是判别一个系统是否作简道振动的依椐．4．圆频率、周期和频率之间的关系，，（9），，三者不是独立的，只要知道其中一个，就可以由（7.9）求出其余两个。
它们是由振动系统的固有性质决定，常称为固有圆频率，固有周期和固有频率．5．振幅和初周相（3）中和是两个积分常数，可由初始条件决定．将初始条件：“，，”代入（3）和（7），得（10）解得（11）求解质点作简谐振动的具体运动情况，也就是要确定（7.3）中的，，三个值．其中和由初始条件决定，因此一般来说，首先必须确定初始值和，而根据（7.10）或（7.11）求出和值．至于（或或），它是由系统固有性质决定的，与初始情况无关．例如对于弹簧振子，，完全由弹簧劲度系数和物体质量所决定．弹簧的大（即所谓硬的弹簧），振动的圆频率也就大。
北科大matlab数学实验报告第二次
《数学实验》报告实验名称 MATLAB绘图学院材料科学与工程学院专业班级材料0801 姓名 xxxxxxx 学号 xxxxxxxxx2011年10月一、【实验目的】学会用MATLAB绘制二维曲线、三维曲线、三维曲面，掌握gtext, legend, title, xlabel, ylabel等指令用法。
二、【实验任务】第四章3、5、9题。
三、【实验程序】1.x1=-pi:pi/100:pi;y1=x1.*cos(x1);x2=pi:pi/100:4*pi;y2=x2.*tan(1./x2).*sin(x2.^3);x3=1:0.1:8;y3=exp(1./x3).*sin(x3);subplot(131),plot(x1,y1,'rp'),grid onaxis tighttitle('y1=x1cosx1')xlabel('x轴'),ylabel('y轴')gtext('y1=x1cosx1')legend('y1=x1cosx1')subplot(132),plot(x2,y2,'c-'),grid onaxis tighttitle('y2=x2tan(1/x2)sin(x2^3)')xlabel('x轴'),ylabel('y轴')gtext('y2=x2tan(1/x2)sin(x2^3)')legend('y2=x2tan(1/x2)sin(x2^3)')subplot(133);plot(x3,y3,'mx'),grid onaxis tighttitle('y3=e^(1/x3)sin(x3)')xlabel('x轴'),ylabel('y轴')gtext('y3=e^(1/x3)*sin(x3)')legend('y3=e^(1/x3)sin(x3)')2. t=0:pi/10:20*pi;x=t.*cos(pi/6.*t);uy=t.*sin(pi/6.*t);z=2*t;plot3(x,y,z),grid onxlabel('x轴 x=tcos(pi/6t)')ylabel('y轴 y=tsin(pi/6t)')zlabel('z轴 z=2t')title('圆锥螺线（0≤t≤20pi）')legend('圆锥螺线')3. %三维曲面t=-2:0.01:2;[x,y]=meshgrid(t);z1=5-x.^2-y.^2+eps;subplot(131),mesh(x,y,z1),title('三维曲面')%平面subplot(1,3,2)z2=3*ones(size(x))mesh(x,y,z2)title('平面')%交线r0=abs(z1-z2)<=1;zz=r0.*z2;yy=r0.*y;xx=r0.*x;subplot(1,3,3),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'m-')title('交线')四、【实验结果】1.2.3.五、【实验总结】这次作业明显比第一次轻松了许多，在上课认真听讲后，又研读了书上的讲解和例题，起初是模仿，后来自己就能完全理解并做出题目来。
北科数理统计与Matlab上机报告2
区间上限
6
【练习 2_03】
a1=[0.143,0.142,0.143,0.137]; a2=[0.140,0.142,0.136,0.138,0.140]; b1=mean(a1); b2=mean(a2); m1=var(a1); m2=var(a2); n1=length(a1); n2=length(a2);
切比雪夫不等式-----------------------------------------------------------------------样本数 样本均值 标准差 区间半径 区间下限 ------------------------------------------------------------------------------------10 1147.0000 87.0000 123.0366 1023.9634 1270.0366 ------------------------------------------------------------------------------------区间上限
1
2520 3540 2600 3320 3120 3400 2900 2420 3280 3100 2980 3160 3100 3460 2740 3060 3700 3460 3500 1600 3100 3700 3280 2880 3120 3800 3740 2940 3580 2980 3700 3460 2940 3300 2980 3480 3220 3060 3400 2680 3340 2500 2960 2900 4600 2780 3340 2500 3300 3640
已知方差的正态分布-------------------------------------------------------------------样本数 样本均值 标准差 区间半径 区间下限 ------------------------------------------------------------------------------------10 1147.0000 87.0000 53.9222 1093.0778 1200.9222 -------------------------------------------------------------------------------------
北京科技大学控制工程基础Matlab大作业
《控制工程基础》实验报告学院：工程学院专业：XXX姓名：XXX学号：XXX班级：2015年1月1日【实验一】1、系统传递函数G（s）,求出极点，画出零极点图，试判断其稳定性。
求解过程：【Matlab程序】num=[9,1.8,9];den=[1,3.2,11.4,18]；G=tf(num,den);pzmap(G),p=roots(den)图1-1 G（s）的零极点分布图【运行结果】p =-0.6000 + 2.9394i-0.6000 - 2.9394i-2.0000由计算结果可知，该系统的三个极点都位于S平面的左半平面，故系统稳定。
2、二阶系统（）,画出单位阶跃响应曲线，计算系统的闭环根，阻尼比，无阻尼震荡频率，记录实际测取的峰值时间tp上升时间tr调整时间ts 超调量σ%，并在所得图形上加网格线和标题。
求解过程：【Matlab程序】：G=zpk([],[-4+10.2i,-4-10.2i],[120]);C=dcgain(G);[y,t]=step(G);plot(t,y)grid on ;[Y,k]=max(y);timetopeak=t(k)percentovershoot=100*(Y-C)/Cn=1;while y(n)<Cn=n+1;endrisetime=t(n)i=length(t);while(y(i)>0.98*C)&(y(i)<1.02*C)i=i-1;endsetllingtime=t(i)图1-2 系统单位阶跃响应曲线【运行结果】timetopeak = 0.3016percentovershoot = 29.0989risetime =0.1967setllingtime =0.9836所以，峰值时间tp=0.3016,最大超调量σ%为29.0989%，上升时间t r为0.1967，调整时间t s为0.9836（2）计算系统的闭环根，阻尼比，无阻尼振荡频率num=[120];den=[1 8 120];G=tf(num ,den);[wn,z,p]=damp(G)【运行结果】wn =10.954510.9545z =0.36510.3651p =-4.0000 +10.1980i-4.0000 -10.1980i由上面的计算结果得，系统的闭环根为-4±10.1980i，阻尼比ζ=0.3651，无阻尼振荡频率为w=10.9545n【实验二】1、求解过程【Matlab程序】%画系统奈奎斯特曲线图sys=tf([0.316 31.6],[0.005 0.155 1.15 1 0]);figure(1);nyquist(sys);grid on;title('niquist plot of G(s)=31.6(0.1s+1)/s(s+1)(0.1s+1)(8s+1)')图2-1 系统奈奎斯特曲线%画系统伯德图%Matlab程序sys=tf([0.316 31.6],[0.005 0.155 1.15 1 0]);figure(1);bode(sys);grid on;title('Bode Digram of G(s)=31.6(0.1s+1)/s(s+1)(0.1s+1)(8s+1)')图2-2 系统伯德图%计算系统稳定裕度Lg和γcsys=tf([0.316 31.6],[0.005 0.155 1.15 1 0]);margin(sys);grid on图2-3 系统稳定裕度值由上图可知，系统稳定幅值裕度Lg为-12.2dB和γc为-27.92、【求解过程】【%Matlab程序】Num1=[1];den1=[1];Num2=[3.333 1];den2=[100 1];Num3=[100];den3=[0.005 0.15 1];Gc1=tf(num1, den1);Gc2=tf(num2,den2);G=tf(num3,den3);G11=series(Gc1,G);G22=series(Gc2,G);figure;bode(G,G11);grid on;title(‘伯德图曲线比较’);figure;bode(G,G22); grid on;title(‘伯德图曲线比较’);figure;bode(G,G11,G22); grid on;title(‘伯德图曲线比较’);figure;margin(G);grid on;figure;margin(G11);grid on;figure;margin(G22);grid on;图2-4 原系统与加入控制器Gc（s）=1系统的伯德图比较图2-5 原系统与加入控制器Gc（s）=(3.333s+1)/(100s+1)系统的伯德图比较图2-6 原系统分别与加入两个控制器系统的伯德图比较图2-7 原系统稳定裕度值图2-8 加入控制器Gc（s）=1系统稳定裕度值可见系统稳定性不变。
（完整版）MATLAB作业2参考答案
（完整版）MATLAB作业2参考答案MATLAB 作业⼆参考答案1、试求出如下极限。
【求解】极限问题可以由下⾯语句直接求解。
>> syms x; f=(x+2)A(x+2)*(x+3)A(x+3)/(x+5F(2*x+5); limit(f,x,i nf) ans = exp(-5) >> syms x yfa=(x A 2*y+x*y A 3)/(x+y)A 3; limit(limit(fa,x,-1),y,2) ans = -6>> fc=(1-cos(xA2+yA2))*exp(xA2+yA2)/(xA2+yA2); limit(limit(fc,x,O),y,O) ans = 02、试求出下⾯函数的导数。
(1) y(x) . xsinx 、1e x , (2) atan yIn(x 2 y 2)x【求解】由求导函数diff() 可以直接得出如下结果，其中 (2)为隐函数，故需要⽤隐函数求导公式得出导数。
>> syms x;f=sqrt(x*si n( x)*sqrt(1-exp(x))); simple(diff(f)) ans =1/2/(x*s in (x)*(1-exp(x))A(1/2))A(1/2)*(si n(x)*(1-exp(x))A(1/2)+ x*cos(x)*(1-exp(x))A(1/2)-1/2*x*sin (x)/(1-exp(x))A(1/2)*exp(x)) >> syms x,y; f=ata n(y /x)_log(xA2+yA2); f1=simple(-diff(f,x)/diff(f,y)) f1 =(y+2*x)/(x-2*y)2 23、假设U cos 1 J *，试验证⼀———。
\y x y y x【求解】证明⼆者相等亦可以由⼆者之差为零来证明，故由下⾯的语句直接证明。
>> syms x y; u=acos(x/y); diff(diff(u,x),y)-diff(diff(u,y),x) ans =(1) limx (x 2)x 2(x 3)x3(x 5)2x 52..x y lim y 21(x xy 3y)3H xy2Xos2Xe2o o2xyt 2x2f2f2f4、假设 f (x, y) e dt ，试求 22 2y x x yy【求解】由下⾯的命令可以得出所需结果。
matlab数学实验报告
MATLAB数学实验报告指导老师：班级：小组成员：时间：201_/_/_Matlab第二次实验报告小组成员：1题目：实验四;MATLAB选择结构与应用实验目的：掌握if选择结构与程序流程控制;重点掌握break;return;pause语句的应用..问题：问题1：验证“哥德巴赫猜想”;即：任何一个正偶数n>=6均可表示为两个质数的和..要求编制一个函数程序;输入一个正偶数;返回两个质数的和..问题分析：由用户输入一个大于6的偶数;由input语句实现..由if判断语句判断是否输入的数据符合条件..再引用质数判断函数来找出两个质数;再向屏幕输出两个质数即可..编程：function z1;z2=geden;n=input'please input n'if n<6disp'data error';returnendif modn;2==0for i=2:n/2k=0;for j=2:sqrtiif modi;j==0k=k+1;endendfor j=2:sqrtn-iif modn-i;j==0k=k+1;endendif k==0fprintf'two numbers are'fprintf'%.0f;%.0f';i;n-ibreakendendend结果分析：如上图;用户输入了大于6的偶数返回两个质数5和31;通过不断试验;即可验证哥德巴赫猜想..纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰;更快的解决问题..2题目：实验四;MATLAB选择结构与应用实验目的：用matlab联系生活实际;解决一些生活中常见的实际问题..问题：问题四：在一边长为1的四个顶点上各站有一个人;他们同时开始以等速顺时针沿跑道追逐下一人;在追击过程中;每个人时刻对准目标;试模拟追击路线;并讨论.. （1）四个人能否追到一起（2）若能追到一起;每个人跑过多少路程（3）追到一起所需要的时间设速率为1问题分析：由正方形的几何对称性和四个人运动的对称性可知;只需研究2个人的运动即可解决此问题..编程：hold onaxis0 1 0 1;a=0;0;b=0;1;k=0;dt=0.001;v=1;while k<10000d=norma-b;k=k+1;plota1;a2;'r.';'markersize';15;plotb1;b2;'b.';'markersize';15;fprintf'k=%.0f b%.3f;%.3f a%.3f;%.3f d=%.3f\n';k;b1;b2;a1;a2;da=a+b1-a1/d*dt;b2-a2/d*dt;b=b+b2-a2/d*dt;-b1-a1/d*dt;if d<=0.001breakendendfprintf'每个人所走的路程为：%.3f';k*v*dtfprintf'追到一起所需要的时间为%.3f';k*dt结果分析：上图为2人的模拟运动路线;有对称性可解决所提问题..-上图为运算过程和运算结果..四个人可以追到一起;走过的路程为1.003;时间也为1.003.纪录：此题利用正方形和运动的对称性可以简便运算..3题目：实验八;河流流量估计与数据插值目的：由一些测量数据经过计算处理;解决一些生活实际问题..问题：实验八上机练习题第三题：瑞士地图如图所示;为了算出他的国土面积;做以下测量;由西向东为x轴;由南向北为y轴;从西边界点到东边界点划分为若干区域;测出每个分点的南北边界点y1和y2;得到以下数据mm..已知比例尺1:2222;计算瑞士国土面积;精确值为41288平方公里..测量数据如下：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158 ;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;问题分析：先由题目给定的数据作出瑞士地图的草图;再根据梯形法;使用trapz语句;来估算瑞士国土的面积..编程：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;plotx;y1;'r.';'markersize';15;plotx;y2;'r.';'markersize';15;axis0 160 0 135grid;hold ont=7:158;u1=splinex;y1;t;u2=splinex;y2;t;plott;u1plott;u2s1=trapzt;u1;s2=trapzt;u2;s=s2-s1*2222*22222/10000000;fprintf'S=%.0f';s结果分析：上图为由所给数据绘制出的瑞士地图..上图为运算结果;计算出瑞士的国土面积为42472平方公里;与准确值41288较为接近..纪录：使用梯形分割的方法;trapz语句可以方便计算不规则图形面积;但存在一定误差..4题目：实验七：圆周率的计算与数值积分目的：将数值积分最基本的原理应用于matlab之中;解决一些与积分有关的问题..问题：实验七上机练习题第一题：排洪量某河床的横断面如图7.3所示;为了计算最大排洪量;需要计算其断面积;试根据所给数据m用梯形法计算其断面积..问题分析：河床断面可近似分割成若干曲边梯形;近似处理把它们当做梯形来计算面积可使问题得到简化..编程：clc;clear;x=0 4 10 12 15 22 28 34 40;y=0 1 3 6 8 9 5 3 0;y1=10-y;plotx;y1;'k.';'markersize';15;axis0 40 0 10;grid;hold ont=0:40;u=splinex;y1;t;plott;u;s=40*10-trapzt;u;fprintf's=%.2f\n';s结果分析：上图为河床的断面图..上图为计算结果面积约为180.70平方米..纪录：使用梯形法计算不规则图形面积十分简便易行..5题目：实验七：圆周率的计算与数值积分目的：使用matlab计算解决一些有关积分的问题..问题：实验七上机练习题第三题：从地面发射一枚火箭;在最初100秒内记录其加速度如下;试求火箭在100秒时的速度..Ts=0 10 20 30 40 50 60 70 80 90 100;Am/s*s=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.67 54.01 57.23;问题分析：加速度为速度的微分;已知微分求积分;类似于面积问题;可使用梯形法来计算..编程：clc;clear;x=0 10 20 30 40 50 60 70 80 90 100;y=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.6754.01 57.23;plotx;y;'k.';'markersize';15;axis0 100 20 60;grid;hold ons=0:10:100;z=splinex;y;s;plots;y;v=trapzx;y;fprintf'v=%.2f\n';v结果分析：上图为加速度变化图..上图为计算结果;求得火箭在100秒时速度约为4168.95m/s..纪录：梯形法可以推广解决许多已知微分求积分的其他问题..6题目：实验七：圆周率的计算与数值积分目的:计算曲线弧长闭曲线周长可使用微元法;ds=sqrtdx^2+dy^2;在转化微积分问题;累加即可得到结果..问题：实验七上机练习题第三题：计算椭圆想x^2/4+y^2=1的周长;使结果具有五位有效数字..问题分析：编程：s=0;dx=0.001;for x=0:0.001:1.999dy=1.-x+0.001.^2/4-1.-x.^2/4;ds=sqrtdx.^2+dy.^2;s=s+ds;ends=4*s;fprintf'the length is'fprintf'%.4f';s结果分析：上图为计算结果;给定椭圆的周长约为9.1823五位有效数字纪录：计算不规则曲线弧长;可使用微元法;划分为若干小的看做直角三角形;利用勾股定理解决..7题目：实验九人口预测与数据拟合目的：掌握一些曲线拟合的方法;了解曲线拟合常用函数..问题：用电压U=10v的电池给电容器充电;t时刻的电压Vt=U-U-V0exp-t/τ;其中V0是电容器的初始电压;τ是充电常数;由所给数据确定V0和τ..t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;问题分析：题中已给出函数关系式;为指数函数曲线拟合;将所给函数式整理可得标准的exp形函数曲线;从而便于解决..编程：t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;plott;V;'k.';'markersize';20;axis0 10 0 4;grid;hold onpause0.5n=8;a=sumt1:n;b=sumt1:n.*t1:n;c=sumlogV1:n;d=sumt1:n.*logV1:n;A=n a;a b;B=c;d;p=invA*Bx=0:10;y=expp1+p2*x;plotx;y;'r-';'linewidth';2结果分析：上图为电压与时间关系图..上图为计算结果;即U-V0=1.4766;所以V0=8.5234;-1/τ=-0.2835;所以τ=3.5273纪录：曲线拟合的一个重难点是选择合适的曲线函数;才能提高拟合度..8题目：实验七圆周率的计算与数值积分目的：拓展圆周率的各种计算方法;掌握其他数值的近似计算方法..问题：实验七练习2：计算ln2的近似值精确到10的-5次方（1）利用级数展开的方法来计算（2）利用梯形法计算（3）利用抛物线法问题分析：级数展开;梯形法;抛物线法是常见的近似运算方法..编程：1级数展开的方法clc;clear;n=0;r=1;p=0;k=-1;while r>=0.1e-5n=n+1;k=k*-1;p1=p+k/n;r=absp1-p;fprintf'n=%.0f;p=%.10f\n';n;p1;p=p1;end2梯形法clc;clear;f=inline'1./x';x=1:0.1:2;y=fx;p=trapzx;y;fprintf'p=%.6f\n';p3抛物线法clc;clear;f=inline'1./x';a=1;b=2;n=1;z=quadf;a;b;fprintf'z=%.10f\n';z结果分析：（1）级数展开的方法（2）梯形法3抛物线法纪录：级数展开法;梯形法;抛物线法;计算近似值时应合理利用..梯形法和抛物线法不易提高精确度;级数展开法可以提高精确度..9题目：实验八河流流量估计与数据插值目的：掌握求插值多项式的方法;并利用此计算近似值..问题：已知y=fx的函数表如下x=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382;求四次拉格朗日插值多项式;并由此求f0.596问题分析：利用所给函数表可计算拉格朗日插值多项式..编程：function p=lagrangex;yL=lengthx;a=onesL;for j=2:La:;j=a:;j-1.*x';endx=inva*y';for i=1:Lpi=xL-i+1;endx=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382; plotx;y;'k.';'markersize';15axis0 2 0 2grid;hold on;p=lagrangex;y;t=0:0.1:1.5;u=polyvalp;t;plott;u;'r-'a=polyvalp;0.596结果分析：上图为所求结果;估算值和插值多项式..纪录：插值多项式是一项十分实用的方法..10题目：求正整数n的阶乘：p=1*2*3*…*n=n;并求出n=20时的结果目的：练习使用循环变量解决数学问题问题：对程序：Clear;clc;n=20;p=1;for i=1:np=p*i;fprintf’i=%.0f;p=%.0f\n’;i;pend进行修改使它：利用input命令对n惊醒赋值问题分析：题中给出程序中“n=20”修改;使用input命令；讲题中的输出命令放出循环之外..编程：clear;clc;n=input'n=';p=1;for i=1:np=p*i;endfprintf'i=%.0f;p=%.0f\n';i;p结果：n=20i=20;p=2432902008176640000>>结果分析：使用input命令可以实现人机对话;使用户自由赋值；输出语句在程序中的位置对输出的结果有很大的影响；在循环内部可以在计算过城中不断输出结果；在循环之外则可以控制只输出最后结果..11题目：对于数列｛√2｝;n=1;2;…;求当其前n项和不超过1000时的n的值及合的大小..目的：运用条件循环解决文帝个项数的循环程序求解;问题：对程序：clear;clc;n=0;s=0;while s<=1000n=n+1;s=s+sqrtn;fprintf’n=%.0f;s=%.4f\n’;n;send问题分析：题中所给程序中的限制变量为上次循环之后的s;导致s超过上限后仍有一次的循环;若把循环变量改为这次的s;则可以避免这种情况的发生..编程：clear;clc;n=0;s=0;while s+sqrtn<=1000n=n+1;s=s+sqrtnfprintf'n=%.0f;s=%.4f\n';n;send结果：……s =970.8891n=128;s=970.8891s =982.2469n=129;s=982.2469s =993.6487n=130;s=993.6487>>结果分析：从结果中可以看出;最后一步为我们需要的答案;从这道题我们可以得出循环变量对一道编程的重要性..。
MATLAB-第二次实验课课堂作业
M A T L A B-第二次实验课课堂作业(4学时)(共4页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--MATLAB 第二次实验课课堂作业（4学时）注：1）此课堂作业作为本课程结业成绩的重要依据，请同学们认真、独立完成，不得抄袭。
2）请在授课教师规定的时间内完成；3）完成作业后，请以word格式保存，文件名为：学号+姓名1.请完成如下4小题：1) 在区间[0,1]s上生成均匀分布的长度为20的向量A。
2) 请生成均值为、方差为的5阶正态分布随机矩阵B。
2) 生成长度为100的单位阶跃序列，阶跃点n0＝25。
用离散信号图示。
3）幅度为1V，频率为5Hz的正弦信号和（0，）的白噪声信号的叠加。
信号长1秒。
设采样频率为100Hz。
将正弦信号和白噪声信号表示在一张图形框，不同坐标系内。
另起一个图形框对叠加信号作图。
要求: 1. 给出1和2 小题的矩阵A 和B2. 给出3和4小题中的2组图形3. 给出你的程序内容2. 已知y=conv(x,h)能够完成输入信号x（n）与系统单位样值响应h（n）的卷积操作，得到输出信号y（n）。
请通过调用子函数y=conv(x,h)，编制一个程序 [y ny ]=conv1(x,h,nx,nh)，其中ny是y （n）第一个样值的n值， nx是x（n）第一个样值的n值， nh是h（n）第一个样值的n值.完成: 1）若 x（n）=[3 1 21 5 2 ] ,nx =-1; h（n）=[1 1 3 ] ,nh=-2 ,则y（n）和ny为多少？请用你编制的函数conv1完成。
2）给出conv1的程序内容解1）程序clc;clear;closeA=linspace(0,1,20)A=Columns 1 through 11Columns 12 through 202）程序clc;clear;close allB=+sqrt*randn(5)B =c3)程序内容Clc;clear;close all N=100;n0=25;u=ones(1,N);u(1,1:n0)=0;stem(0:N-1,u) 4)00.51-1-0.8-0.6-0.4-0.20.20.40.60.8100.51-0.25-0.2-0.15-0.1-0.050.050.10.150.200.10.20.30.40.50.60.70.80.91-1.5-1-0.50.511.5程序 clc;clear;close allts=0;te=1;fs=100 ;f1=5;t=ts:1/fs:te;y1=sin(2*pi*f1*t);y2=sqrt*randn(size(t));subplot(1,2,1);plot(t,y1);subplot(1,2,2);plot(t,y2);figure(2)plot(t,y1+y2)2.1)y =3 4 31 29 70 17 6ny =-32) function [y,ny]=conv1(x,nx,h,nh)y =conv(x,h);ny=nx+nh;。
北京科技大学《matlab和数学实验》 第二次课
数学实验
例1 画出函数y = sin x2 在-5 x 5 的图形。
解: Matlab 命令: x=-5:0.1:5;↙ y=sin(x.^2); ↙ plot(x,y),grid on
X=a:h:b; 步长h可以任意选取， 步长越小，曲线越光滑。
grid on作用是：在图形中标出网格线
例 x=[ 3 4 7]; y=[4 5 6 5 4795 1 4 2 5]; plot (x,y)
6
数学实验
• 思考： 1、当x为n维向量，y为m x n 矩阵时， plot(x,y)怎么画？ 2、当x为n维向量，y为m x w 矩阵时， plot(x,y)怎么画？ 3、当x为n维向量，y为n x n 矩阵时， plot(x,y)怎么画？
数学实验
• Mesh命令介绍
– 利用函数meshgrid生成格点矩阵后，然后 求出各格点对应的函数值，就可以利用三维 网格命令mesh画出空间曲面。
– 函数mesh用来生成函数的网格曲面，即各 网格线段组成的曲面。
25
例:x=[1,2,3] y=[1,2,3,4] [X,Y]=meshgrid(x,y)
19
数学实验
• 例7 在同一坐标系中画出两个函数y = cos 2x ,y = x的图形, 自变量范围为:-2 x 2,函数y = cos 2x为红色实线, 函数y = x 为蓝色点划线。并加注标题，坐标轴。 解： Matlab命令为 clf x=-2:.1:2; y1=cos(2*x);y2=x; subplot(1,2,1) plot(x,y1,'r-',x,y2,'b-.'),grid on title('曲线y1=cos(2x)与y2=x') xlabel('x轴'),ylabel('y轴') gtext('y1=cos(2x)'),gtext('y2=x') legend('y1=cos(2x)','y2=x') subplot(1,2,2) plot (x,y1,x,y2)
matlab2022实验2参考答案
matlab2022实验2参考答案报告名称：MATLAB试验二符号计算姓名：学号：专业：班级：MATLAB实验二MATLAB符号计算试验报告说明：1做试验前请先预习，并独立完成试验和试验报告。
2报告解答方式：将MATLAB执行命令和最后运行结果从命令窗口拷贝到每题的题目下面，请将报告解答部分的底纹设置为灰色，以便于批阅。
3在页眉上写清报告名称，学生姓名，学号，专业以及班级。
3报告以Word文档书写。
一目的和要求1熟练掌握MATLAB符号表达式的创建2熟练掌握符号表达式的代数运算3掌握符号表达式的化简和替换4熟练掌握符号微积分5熟练掌握符号方程的求解二试验内容1多项式运算（必做）1.1解方程:f(某)=某^4-10某某^3+34某某^2-50某某+25=0%采用数值方法：>>f=[1-1034-5025];>>root(f)%采用符号计算方法：f1=ym('某^4-10某某^3+34某某^2-50某某+25')olve(f1)1.2求有理分式R=(3某^3+某)(某^3+2)/((某^2+2某-2)(5某^3+2某^2+1))的商多项式和余多项式.a1=[3010];a2=[1002];a=conv(a1,a2);b1=[12-2];b2=[5201];b=conv(b1,b2);[p,r]=deconv(a,b);%注意：ab秩序不可颠倒。
%reidue用于实现多项式的部分分式展开，此处用deconv函数报告名称：MATLAB试验二符号计算姓名：学号：专业：班级：%%此题，有同学程序如下：某1=[3010]，某2=[1002],某3=[12-2],某4=[5201]某5=conv(某1,某2)[y6,r]=deconv(某5,某3)R=deconv(y6,某4)%%这种方法较第一种解法缺点：在除法运算中，会产生误差，故此题应先将分母的多项式相乘后，再与分子部分的多项式进行运算。
最新Matlab数学实验第二版课后习题答案
Matlab数学实验第二版课后习题答案胡良剑第一章%Page20,ex1(5) 等于[exp(1),exp(2);exp(3),exp(4)](7) 3=1*3, 8=2*4(8) a为各列最小值，b为最小值所在的行号(10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture(11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10)(12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10)%Page20, ex2(1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b(2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码%Page20,ex3>> r=2;p=0.5;n=12;>> T=log(r)/n/log(1+0.01*p)T =11.5813%Page20,ex4>> x=-2:0.05:2;f=x.^4-2.^x;>> [fmin,min_index]=min(f)fmin =-1.3907 %最小值min_index =54 %最小值点编址>> x(min_index)ans =0.6500 %最小值点>> [f1,x1_index]=min(abs(f)) %求近似根--绝对值最小的点f1 =0.0328x1_index =24>> x(x1_index)ans =-0.8500>> x(x1_index)=[];f=x.^4-2.^x; %删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) %求另一近似根--函数绝对值次小的点f2 =0.0630x2_index =65>> x(x2_index)ans =1.2500%Page20,ex5>> z=magic(10)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6679 6 13 95 97 29 31 38 45 7210 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59>> sum(z)ans =505 505 505 505 505 505 505 505 505 505>> sum(diag(z))ans =505>> z(:,2)/sqrt(3)ans =57.157746.188046.765450.229553.693613.85642.88683.46416.928210.3923>> z(8,:)=z(8,:)+z(3,:)z =92 99 1 8 15 67 74 51 58 4098 80 7 14 16 73 55 57 64 414 81 88 20 22 54 56 63 70 4785 87 19 21 3 60 62 69 71 2886 93 25 2 9 61 68 75 52 3417 24 76 83 90 42 49 26 33 6523 5 82 89 91 48 30 32 39 6683 87 101 115 119 83 87 101 115 11910 12 94 96 78 35 37 44 46 5311 18 100 77 84 36 43 50 27 59 1小时单位负责人接到报告后，应当于1小时内向事故发生地县级以上人民政府安全生产监督管理部门和负有安全生产监督管理职责的有关部门报告。
数学实验第二次测验题及参考答案（09级）
数学实验第二次测验题及参考答案（09级）数学实验第二次测验题及参考答案一、写出下列MATLAB指令的运算结果.1. A=[1;2;3]; transpose(A)1 2 31 2 32. A=[1, 2, 3 ; 4, 5, 6 ]; B=A([1 2], [1,3]) , d=size(A)B =1 34 6d =2 33. a=1:3; b=linspace(1,3,3); x=sum(a.*b), y=cross(a, b)x = 14y = 0 0 04. A=[1,2, 3; 4,5,6; 7,8,9]; B=ones(3); C=A-BC =0 1 23 4 56 7 85. v=[1, 2, 3]; A=diag(v); E=eig(A), D=det(A)E=123D =66. x=[1,2,3,4,5]; [mean(x), median(x), range(x),sum(x), prod(x)]3 34 15 1207. x=[2,3,4]; a=cumsum(x) ,b=sort(x)a =2 5 9b =2 3 48．format rat; v=[1, 2, 3]; A=diag(v); inv(A)ans =1 0 00 1/2 00 0 1/39. [m,v]=normstat(1,4) % 求参数为1,4的正态分布的均值与方差m =1, v =16二、写出下列MATLAB指令的实验目的.1. dsolve('x*Dy+y-exp(-x)=0', 'y(1)=2*exp(1)', 'x')求微分方程0=-+'-x e y y x 在初始条件e y x 2|1==下的特解.2. u=[1,2,3],v=[0,3,2], w=[5, 2, 1]; dot(w, cross(u, v))计算向量u, v, w 的混合积.3. A=[1 2 3; 2 2 5; 3 5 1]; b=[1;2;3]; det(A); inv(A)*b利用逆矩阵解线性方程组=++=++=++3532522132321 321321x x x x x x x x x .4. A=[0 0 1; 0 1 1; 1 1 1; 1 0 0]; rref(A)求向量组)1,0,0(1=α,)1,1,0(2=α,)1,1,1(3=α,)0,0,1(4=α的秩.或对矩阵A 做行初等变换。
matlab第二次实验报告
《数学建模实验》上机实验报告班级：计算机35班学号：2130505099姓名：田博文【实验一】一上底面半径2米、下底面半径4米、高4米的圆台形水池内盛满了水，由池底一横截面积为0.001平方米的小孔放水。
求在任意时刻的水面高度和将水放空所需的时间。
一、问题分析：在dt很小一段时间内，水的流速v可以看成是不变的与水面高度h有关的函数，即v=gh2；在dt时间内流出的水量，可以近似的为一圆柱，对应的圆柱高为dh；由于容器下降的水量与流出的水量相等，令小孔面积为b，即dV=π2r dh=bvdt；可得dh与dt的关系为：dh=bvdt/(π2r)；根据相似三角形原理可得，r=4-h/2；初始条件有；h(0)=4；为表示方便，引入变量L=4-h故r=2+L/2二、代码：T=0;for L=0:0.001:4-0.001 %选取积分微元为0.001，即dhV=0.001*pi*(2+(1/2*L))^2; %水面处的dV，b=0.001t=V/(sqrt(2*(4-L)*9.8)*0.001) %dV与dt关系T=T+t;plot(T,(4-L),'red'); hold onend三、运行结果：四、结论：在任意时刻的水面高度和将水放空所需的时间如上图红线所示，横轴为时间，纵轴为水面高度。
【实验二】有A、B、C三个场地，每一个场地都出产一定数量的原料，同时也消耗一定数量的产品，具体数据如下表所示。
已知制成每吨产品需要消耗3吨原料，A、B两地，A、C两地和B、C两地之间的距离分别为150千米、100千米和200千米，假设每万吨原料运输1千米的运费为5000元，每万吨产品运输1千米的运费为6000元。
由于地区条件的差异，在不同地区设厂的费用不同，由于条件的限制，在B处建厂的规模不能超过5万吨，问：在这三地如何建厂、规模建多大才能使得总费用最小？地点年产原料（万吨）年销产品（万吨）生产费用（万元/万吨）A 20 7 150B 16 13 120C 24 0 100一、问题分析：设nij为i地运往j地的原料量，mij为i地运往j地的产品量，设A地为1地，B 为2地，C地为3地。
北科大Matlab_数学实验报告1~6次(全)
《数学实验》报告实验名称 Matlab 基础知识学院专业班级姓名学号2014年 6月一、【实验目的】1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。
2.掌握Matlab基本操作和常用命令。
3.了解Matlab常用函数，运算符和表达式。
4.掌握Matlab工作方式和M文件的相关知识。
5.学会Matlab中矩阵和数组的运算。
二、【实验任务】P16 第4题编写函数文件，计算1!nkk =∑，并求出当k=20时表达式的值。
P27第2题矩阵A=123456789⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，B=468556322⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦，计算A*B,A.*B,并比较两者的区别。
P27第3题已知矩阵A=5291⎡⎤⎢⎥⎣⎦，B=1292⎡⎤⎢⎥⎣⎦，做简单的关系运算A>B,A==B,A<B,并做逻辑运算(A==B)&(A<B),(A==B)&(A>B)。
P34 第1题用11114357π=-+-+……公式求π的近似值，直到某一项的绝对值小于-610为止。
三、【实验程序】P16 第4题function sum=jiecheng(n) sum=0;y=1;for k=1:nfor i=1:ky=y*i;endsum=sum+y;endsumP27第2题>>A=[1 2 3;4 5 6;7 8 9]>>B=[4 6 8;5 5 6;3 2 2]>>A*BP27第3题>> A=[5 2;9 1];B=[1 2;9 2];>>A>B>>A==B>>A<B>> (A==B)&(A<B)>> (A==B)&(A>B)P34 第1题t=1;pi=0;n=1;s=1;while abs(t)>=1e-6pi=pi+t;n=n+2;s=-s;t=s/n;endpi=4*pi;四、【实验结果】P16 第4题P27第2题两者的区别：A*B是按正规算法进行矩阵的计算, A.*B是对应元素相乘。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《数学实验》报告
实验名称 MATLAB绘图
学院
专业班级
姓名
学号
年月
一、 【实验目的】
了解并学习绘制MATLAB 二维曲线和三维曲线的图形。
二、 【实验任务】
1.绘制)π4,0(),3sin(3
∈=x x e y x 的图像，要求用蓝色的星号画图，并且画出其包络线3
x e y ±=的图像，用红色的点划线画图。
3.在同一图形窗口画三个子图，要求使用指令gtext ，axis ，legend ，title ，xlabel ，ylabel ： （1）)ππ,(,cos -∈=x x x y （2）)π,4π(,sin 1
tan
3∈=x x x
x y （3）]8,1[,sin 1∈=x x e y x
5.绘制圆锥螺线的图像并加各种备注，圆锥螺线的参数方程为：
π]20,0[,t
2z 6π
sin 6π
cos {
∈===t t t y t
t x
三、 【实验程序】
1.
x=0:pi/50:4*pi; y1=exp(x/3).*sin(3*x); y2=exp(x/3); y3=-exp(x/3); plot(x,y1,'b*') hold on
plot(x,y2,'r-.') hold on plot(x,y3,'r-.')
3.
x1=-pi:pi/50:pi; y1=x1.*cos(x1); x2=pi:pi/50:4*pi;
y2=x2.*tan(1./x2).*sin(x2.^3); x3=1:0.01:8;
y3=exp(1./x3).*sin(x3);
subplot(221),plot(x1,y1,'r-'),grid on axis tight
xlabel('x 轴'),ylabel('y 轴') title('y=xcosx')
gtext('y=xcosx')
legend('y=xcosx')
subplot(222),plot(x2,y2,'b--'),grid on
axis tight
xlabel('x轴'),ylabel('y轴')
title('y=xtan(1/x)sinx^3')
gtext('y=xtan(1/x)sinx^3')
legend('y=xtan(1/x)sinx^3')
subplot('position',[0.2,0.05,0.6,0.45]),plot(x2,y2,'m'),grid on axis tight
xlabel('x轴'),ylabel('y轴')
title('y=exp(1/x)sinx')
gtext('y=exp(1/x)sinx')
legend('y=exp(1/x)sinx')
5.
t=0:pi/10:20*pi;
x=t.*cos(pi*t/6);
y=t.*sin(pi*t/6);
z=2*t;
plot3(x,y,z)
四、【实验结果】
1.
3.
5．
五、【实验总结】
了解了MATLAB绘制二维和三维图形的方法，掌握了各种指令和控制方法，并且经过实践运用MATLAB进行了图形绘制。