北京科技大学MATLAB作业3
matlab作业3
信息科学与工程学院Matlab作业**: ***专业: 电子信息工程2班学号: ********** 日期 2015 年 6 月 28 日1.线性方程求解：使用“\”，求解下面的线性系统，计算并显示误差向量。
364152773a b c a b b c ++=+=+= sy3_1.mA=[3,6,4;1,5,0;0,7,7]; B=[1;2;3]; X=A\B w=A*X-B 运行结果：X = -0.5824 0.5165 -0.0879 w =1.0e-15 *-0.1110 0 02.数值积分：使用trapz 或者quad 计算积分5/3x xedx -⎰，我们知道根据分部积分法可以计算这个积分得到/35/355/3003|9|249x x xe ee -----=-+，那么请比较实用trapz 计算得到的结果和直接使用积分后数值计算的结果差异，并显示在命令行窗口中。
Sy3_2.mx=0:0.0001:5; y=x.*exp(-x/3); z0=trapz(x,y) z1=-24*exp(-5/3)+9 w=z0-z1运行结果：z0 = 4.4670z1 = 4.4670 w =-9.3826e-10 3.矩阵求逆：计算矩阵1234⎡⎤⎢⎥⎣⎦的逆，使用inv函数。
然后将结果乘以原始矩阵，检查结果是否为单位阵。
Sy3_3.ma=[1,2;3,4];b=inv(a)c=b*a运行结果：b =-2.0000 1.00001.5000 -0.5000c =1.0000 00.0000 1.00004.多项式拟合：使用xlsread函数分别载入数据randomDatax.xls和randomDatay.xls，使用ployfit函数对该两对变量分别做1，2，3，4，5阶拟合，为了得到比较好的拟合效果，使用polyfit的中心与标尺模式，具体使用方法为[p,S,mu] = polyfit(x,y,n)，它会返回三个参数，接下来使用对应的[y,delta] = polyval(p,x,S,mu)来根据x对y进行拟合。
MATLAB习题及参考答案经典.doc
习题：1， 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
2， 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3， 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。
(应用sin,cos,tan.cot)5， 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574（2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []2965318772546， 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
(应用poly,polyvalm)7， 求解多项式x 3-7x 2+2x +40的根。
(应用roots)8， 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。
(应用poly,polyvalm)9， 计算多项式9514124234++--x x x x 的微分和积分。
(应用polyder,polyint ，poly2sym)10， 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。
(应用x=a\b)11， 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。
(应用pinv) 12， 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ，计算a 的行列式和逆矩阵。
(应用det,inv)13， y =sin(x )，x 从0到2π，∆x =0.02π，求y 的最大值、最小值、均值和标准差。
MATLAB作业3参考答案
MA TLAB 作业三参考答案1、 请将下面给出的矩阵A 和B 输入到MA TLAB 环境中，并将它们转换成符号矩阵。
若某一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵还是数值矩阵。
576516535501232310014325462564206441211346,39636623515212107600774101201724407734737812486721711076815A B ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦【求解】矩阵的输入与转换是很直接的。
>> A=[5,7,6,5,1,6,5; 2,3,1,0,0,1,4; 6,4,2,0,6,4,4; 3,9,6,3,6,6,2; 10,7,6,0,0,7,7; 7,2,4,4,0,7,7; 4,8,6,7,2,1,7]; A=sym(A) A =[ 5, 7, 6, 5, 1, 6, 5] [ 2, 3, 1, 0, 0, 1, 4] [ 6, 4, 2, 0, 6, 4, 4] [ 3, 9, 6, 3, 6, 6, 2] [ 10, 7, 6, 0, 0, 7, 7] [ 7, 2, 4, 4, 0, 7, 7] [ 4, 8, 6, 7, 2, 1, 7]>> B=[3,5,5,0,1,2,3; 3,2,5,4,6,2,5; 1,2,1,1,3,4,6; 3,5,1,5,2,1,2; 4,1,0,1,2,0,1; -3,-4,-7,3,7,8,12; 1,-10,7,-6,8,1,5]; B=sym(B) B =[ 3, 5, 5, 0, 1, 2, 3] [ 3, 2, 5, 4, 6, 2, 5] [ 1, 2, 1, 1, 3, 4, 6] [ 3, 5, 1, 5, 2, 1, 2] [ 4, 1, 0, 1, 2, 0, 1][ -3, -4, -7, 3, 7, 8, 12] [ 1, -10, 7, -6, 8, 1, 5]2、 利用MA TLAB 语言提供的现成函数对习题1中给出的两个矩阵进行分析，判定它们是否为奇异矩阵，得出矩阵的秩、行列式、迹和逆矩阵，检验得出的逆矩阵是否正确。
Matlab作业习题与答案详解（附程序）
Matlab作业习题与答案详解（附程序）clear all;clc;close all;x=-10:0.01:20;y=4*sin(x)./x;ymin=min(y)二、蒙特卡罗算法的数值计算当前的油位高度是2.3米，见图1。
模拟油流进储油罐的过程（图维数任选），请计算罐内油量。
三维的效果图参见图2。
储油罐由两部分组成，中间是圆柱体，两端是球罐体。
（本题简化自2011年UCMCM A题《储油罐的变位识别与罐容表标定》，细节参见原题原题附件2 cumcm2010A.doc。
）图1图2主程序：clc;clear all;close all;center1=[-3.375,0,1.5]; %左球罐中心center2=[3.375,0,1.5]; %右球罐中心n=10000; %每次的撒点数delta=0.02; %层高h=3;en=h/delta;Show; %画出油罐for i=0:en-1x=(rand(1,n)-0.5)*10; %随机生成n个点y=(rand(1,n)-0.5)*h;z=(rand(1,n)*delta+i*delta);Z=[x;y;z];[dis1dis2]=juli(center1,center2,Z); %算出各点对应的距离index=find(((x>-4&x<4)&dis2<1.5)|(x<-4|x>4)&dis1<1.625); %找出在罐内的点plot3(x(index),y(index),z(index),'.k'); %画出在罐内的点drawnowend子程序1：function [dis1 dis2]=juli(a,b,q)d11=q(1,:)-a(1);d12=q(2,:)-a(2);d13=q(3,:)-a(3);d1=sqrt(d11.^2+d12.^2+d13.^2);d21=q(1,:)-b(1);d22=q(2,:)-b(2);d23=q(3,:)-b(3);d2=sqrt(d21.^2+d22.^2+d23.^2);d1(d1>d2)=d2(d1>d2);dis1=d1;dis2=sqrt(d12.^2+d13.^2);子程序2：function tu=Show%===圆柱部分==figure('color','w')h=3;y='3/2*cos(s)';z='3/2*sin(s)+1.5';x='t';ezmesh(x,y,z,[0 pi*2 -4 4]);axis equalhidden offhold on%===左罐部分==z='cos(s)*(2.6406-t.^2).^(1/2)+1.5'; y='sin(s)*(2.6406-t.^2).^(1/2)';x='t-3.375';ezmesh(x,y,z,[0 pi*2 -1.625 -0.625]);axis equalhidden offhold on%===右罐部分===z='cos(s)*(2.6406-t.^2).^(1/2)+1.5';y='sin(s)*(2.6406-t.^2).^(1/2)';x='t+3.375';ezmesh(x,y,z,[0 pi*2 0.625 1.625]);colormap(gray)axis equalhold onaxis offhidden off三、元胞自动机的简单应用1.简单交通流模拟请模拟一个4车道的交通流，车辆密度为0.3，其余规则自行定义。
北航MATLAB第三次作业
tri=0; s=0; for k=1:m for l=1:n if k==l; tri=tri+A(k,l); end s=s+A(k,l); end end tri s
7、试用 for 语句和 if 语句创建下列矩阵
5 1 A 0 0 0
function f7 A=[]; for i=1:5 for j=1:5 if i==j A(i,j)=5; else if abs(i-j)==1 A(i,j)=1; else A(i,j)=0; end end end end A
8. ®•õ‘ªP1 (x) = 3x + 2, • P (x) = P1 P2 P3 • P (x) = 0
P2 (x) = 5x2 − x + 2,
P2 (x) = x2 − 0.5.
ê¶
ÜŠ"2源自作业题答案 第 1 题： >> x=-1.42; >> y=0.52; >> sin(abs(x)+abs(y))/sqrt(abs(x-y)) ans = 0.6696 第 2 题： >> a=1; >> b=2; >> c=3; >> p=[a b c]; >> roots(p) ans = -1.0000 + 1.4142i -1.0000 - 1.4142i 第 3 题： >> x=linspace(0,20,50) x= Columns 1 through 7 0 0.4082 0.8163 1.2245 1.6327 2.0408 2.4490
>> g=[3 2]; >> g1=[0 0 0 3 2]; >> f+g1 ans =
北京科技大学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实验三答案
运算符
（）
^
~
*
> < >= == <=
|&
矩阵进行关系运算时，应该保证两个矩阵的维数一致或其中一 个矩阵为标量。比较关系运算是对两个矩阵的对应运算进行比较。
1.2 租用一个交通工具，前 100 公里需要 0.5 美元一公里，在 下面的 200 公里中 0.23 美元每公里，超过 300 公里的部分一律按 0.2 美元每公里算，编写一个程序，自定义总公里数 x，求总花销及平 均每公里的花销。（分别定义 x<=100,100<x<=300,x>300 三个总公里 数验证结果）。
>> d|b>a
ans =
1
>> a*b^2>a*c
ans =
(2)
a
1
1 2
0 1
>> a=[1 0;-2 1];b=[0 2;-2 -1];
>> a>=b
ans =
>> a|b
ans =
(3)
a>c&b>c, c<=d
10
11
11
11
>> format compact
a
，
2,
M 文件： x=input('请输入总公里数 x='); if x<=100
w=0.5*x; disp(['总花销 w=',num2str(w),' 美元']) elseif 100<x<=300 w=0.5*100+(x-100)*0.23; disp(['总花销 w=',num2str(w),' 美元']) else
北京科技大学控制工程基础matlab大作业
实验一MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性；二、要点内容1、系统的典型响应有哪些？2、如何判断系统稳定性？3、系统的动态性能指标有哪些？三、实验方法（一）四种典型响应1、阶跃响应：阶跃响应常用格式：1、step(sys)；其中sys可以为连续系统，也可为离散系统。
2、step(sys,Tn)；表示时间范围 0---Tn。
3、step(sys,T)；表示时间范围向量 T 指定。
4、Y=step(sys,T)；可详细了解某段时间的输入、输出情况。
2、脉冲响应：∞脉冲函数在数学上的精确定义：∫f(x)dx=1f (x)=0, t 0其拉氏变换为：f(s)=1Y (s)= G(s) f (s)= G(s)所以脉冲响应即为传函的反拉氏变换。
脉冲响应函数常用格式：①impulse(sys)；②impulse(sys,Tn);impulse(sys,T );③ Y = impulse(sys,T )（二）分析系统稳定性有以下三种方法：1、利用 pzmap 绘制连续系统的零极点图；2、利用 tf2zp 求出系统零极点；3、利用 roots 求分母多项式的根来确定系统的极点（三）系统的动态特性分析Matlab 提供了求取连续系统的单位阶跃响应函数step、单位脉冲响应函数 impulse、零输入响应函数initial以及任意输入下的仿真函数lsim.四、实验内容(一) 稳定性%MATLAB计算程序：num=[9 1.8 9];den=[1 3.2 2 18];G=tf(num,den);pzmap(G);p=roots(den) p =-3.8801 + 0.0000i0.3401 + 2.1268i0.3401 - 2.1268i由计算结果可知，该系统的两个极点具有正实部，故系统不稳定。
键入程序，观察并记录单位阶跃响应曲线：num=[120];den=[1 8 120];step(num,den);title('Step Response ofG(s)=120/s^2+8s+120)');>> 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所以系统的闭环根为s=-4+10.1980i， s=-4-10.1980i；阻尼比为0.3651；无阻尼振荡频率为10.9545.由上图实测峰值时间tp=0.292s；最大超调量Mp=1.29-1=0.29，即Mp=29%。
matlab作业3参考答案
matlab作业3参考答案Matlab作业3参考答案Matlab作业3是一个综合性的编程任务，要求学生运用Matlab的各种功能和工具来解决实际问题。
本文将提供Matlab作业3的参考答案，并对其中的关键步骤和思路进行详细解释。
一、问题描述在本次作业中，学生需要解决一个关于图像处理的问题。
具体来说，给定一张彩色图像，学生需要编写Matlab代码来实现以下功能：1. 将彩色图像转换为灰度图像；2. 对灰度图像进行高斯滤波；3. 对滤波后的图像进行边缘检测；4. 对边缘图像进行二值化处理。
二、解决方案1. 将彩色图像转换为灰度图像首先，我们需要读取彩色图像。
可以使用Matlab的imread函数来实现。
然后，使用rgb2gray函数将彩色图像转换为灰度图像。
代码如下：```matlabrgbImage = imread('image.jpg');grayImage = rgb2gray(rgbImage);```2. 对灰度图像进行高斯滤波接下来，我们需要对灰度图像进行高斯滤波。
高斯滤波是一种常用的图像平滑方法，可以有效地去除图像中的噪声。
Matlab提供了fspecial函数来生成高斯滤波器。
代码如下：```matlabh = fspecial('gaussian', [3 3], 1);filteredImage = imfilter(grayImage, h);```3. 对滤波后的图像进行边缘检测在这一步中，我们需要对滤波后的图像进行边缘检测。
边缘检测可以帮助我们找到图像中的边缘和轮廓。
Matlab提供了多种边缘检测算法，如Sobel算子和Canny算子。
代码如下：```matlabedgeImage = edge(filteredImage, 'canny');```4. 对边缘图像进行二值化处理最后，我们需要对边缘图像进行二值化处理，将图像中的边缘转换为黑白两种颜色。
matlab综合大作业(附详细答案)
m a t l a b综合大作业(附详细答案)-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）：1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5)实验结果：A =0.1349 3.3818 0.6266 1.2279 1.5888-2.3312 3.3783 2.4516 3.1335 -1.67241.2507 0.9247 -0.1766 1.11862.42861.5754 1.6546 5.3664 0.8087 4.2471-1.2929 1.3493 0.7272 -0.6647 -0.38362)将矩阵A按列拉长得到矩阵B；实验程序：B=A(:)实验结果：B =0.1349-2.33121.25071.5754-1.29293.38183.37830.92471.65461.34930.62662.4516-0.17665.36640.72721.22793.13351.11860.8087-0.66471.5888-1.67242.42864.2471-0.38363)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)]实验结果：C =3.3783 3.13350.9247 1.11864)寻找矩阵A中大于0的元素；]实验程序：G=A(find(A>0))实验结果：G =0.13491.25071.57543.38183.37830.92471.65461.34930.62662.45165.36640.72721.22793.13351.11860.80871.58882.42864.24715)求矩阵A的转置矩阵D；实验程序：D=A'实验结果：D =0.1349 -2.3312 1.2507 1.5754 -1.29293.3818 3.3783 0.9247 1.6546 1.34930.6266 2.4516 -0.1766 5.3664 0.72721.2279 3.1335 1.1186 0.8087 -0.66471.5888 -1.67242.4286 4.2471 -0.38366)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；实验程序：E=flipud(fliplr(A))实验结果：E =-0.3836 -0.6647 0.7272 1.3493 -1.29294.2471 0.80875.3664 1.6546 1.57542.4286 1.1186 -0.1766 0.9247 1.2507-1.6724 3.1335 2.4516 3.3783 -2.33121.5888 1.2279 0.6266 3.3818 0.13497)删除矩阵A的第2列和第4列得到矩阵F；实验程序：F=A;F(:,[2,4])=[]实验结果：F =0.1349 0.6266 1.5888-2.3312 2.4516 -1.67241.2507 -0.17662.42861.5754 5.3664 4.2471-1.2929 0.7272 -0.38368)求矩阵A的特征值和特征向量；实验程序：[Av,Ad]=eig(A)实验结果：特征向量Av =-0.4777 0.1090 + 0.3829i 0.1090 - 0.3829i -0.7900 -0.2579 -0.5651 -0.5944 -0.5944 -0.3439 -0.1272-0.2862 0.2779 + 0.0196i 0.2779 - 0.0196i -0.0612 -0.5682 -0.6087 0.5042 - 0.2283i 0.5042 + 0.2283i 0.0343 0.6786 0.0080 -0.1028 + 0.3059i -0.1028 - 0.3059i 0.5026 0.3660 特征值Ad =6.0481 0 0 0 00 -0.2877 + 3.4850i 0 0 00 0 -0.2877 - 3.4850i 0 00 0 0 0.5915 00 0 0 0 -2.30249)求矩阵A的每一列的和值；实验程序：lieSUM=sum(A)实验结果：lieSUM =-0.6632 10.6888 8.9951 5.6240 6.208710)求矩阵A的每一列的平均值；实验程序：average=mean(A)实验结果：average =-0.1326 2.1378 1.7990 1.1248 1.24172．符号计算（10分，每小题5分）：1)求方程组20,0++=++=关于,y z的解；uy vz w y z w实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z');y= S. y, z=S. z实验结果：y =[ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] z =[ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))]2)利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解； 实验程序：[x,y]=dsolve('Dx=y','Dy=-x')实验结果：x =-C1*cos(t)+C2*sin(t)y = C1*sin(t)+C2*cos(t)3．数据和函数的可视化（20分，每小题5分）：1)二维图形绘制：绘制方程2222125x y a a +=-表示的一组椭圆，其中0.5:0.5:4.5a =；实验程序：t=0:0.01*pi:2*pi; for a=0.5:0.5:4.5; x=a*cos(t); y=sqrt(25-a^2)*sin(t); plot(x,y) hold on end实验结果：2) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线；实验程序：x=0:0.1:4; y1=sin(x); y2=10.^x;[ax,h1,h2]=plotyy(x,y1,x,y2); set(h1,'LineStyle','.','color','r'); set(h2,'LineStyle','-','color','g'); legend([h1,h2],{'y=sinx';'y=10^x'});实验结果：3)用曲面图表示函数22z x y =+；实验程序：x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; surf(X,Y,Z)实验结果：4)用stem 函数绘制对函数cos 4y t π=的采样序列；实验程序：t=-8:0.1:8;y=cos(pi.*t/4); stem(y)实验结果：4. 设采样频率为Fs = 1000 Hz ，已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=，由于某一原因，原始信号被白噪声污染，实际获得的信号为))((ˆt size randn x x+=，要求设计出一个FIR 滤波器恢复出原始信号。
北科数理统计与Matlab上机报告材料3
效应A303.60475.9015.180.0003
误差50.00105.00**
总和353.6014临界值=3.4780(0.05),5.9943(0.01)
>> lx4_01_01
序号个数求和平方和方差均值标准差
13423842049121143.985770.117633.8229
fprintf('%d \t%d\t %d\t %d\t %.4f\t %.4f\t %.4f\n',oo,ni(oo),S1(oo),S2(oo),S3(oo),S4(oo),S5(oo));
end
fprintf('------------------------------------------------------------------------------------------------------------\n')
la1=finv(1-alpha1,a-1,n-a); %由F分布的累积概率，求临界值，P{F<la}=1-alpha
alpha2=0.01; %显著性水平
la2=finv(1-alpha2,a-1,n-a); %由F分布的累积概率，求临界值，P{F<la}=1-alpha
p=1-fcdf(F,a-1,n-a); %计算F比值做为临界点的右侧概率p=1-P{X<F}
if F>la2
xzx='**';
elseif F>la1
xzx='*';
else
xzx='-';
end
fprintf('序号\t个数\t求和\t平方和\t\t方差\t \t均值\t\t标准差\n');
北京科技大学matlab和数学实验绘图
怎么画？
北京科技大学数学实验
问题：plot(x,y)命令可以用来画通常的函数 f(x)（其中a<x<b)的图像？
可以。 x=a:h:b 函数f(x)在绘图区间[a,b]上的自变量点向
量数据 y=f(x) 对应的函数值向量
步长h可以任意选取，步长越小，曲线越光滑。
北京科技大学数学实验
例1 在区间[- , ]上，绘制函数y = sin(x)图形。
， 例：绘制三维曲线的图像：
解：matlab命令为：
t=0:pi/10:20*pi; x=t.* sin(t); y=t.*cos(t) ; z=t ; plot3(x,y,z)
x t sin t
y
t cos t
z t
(0t20)
北京科技大学数学实验
四、Matlab空间曲面绘图 Meshgrid命令：
思考：
北京科技大学数学实验
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)怎么画？
4、当x，y为n x n 矩阵时，plot(x,y)
例：绘制火柴棍图 x=0:.1:2; y= exp(-x.^2); stem(y)
北京科技大学数学实验
例：绘制阶梯图 x=0:pi/20:2*pi; y=sin(2*x); stairs(x,y)
北京科技大学数学实验
例：极坐标图 polar(t,r) 使用极角t和极径r绘制极坐标图形；
t=0:pi/50:2*pi; polar(t,abs(sin(4*t)),'r')
《Matlab语言及其应用》形考作业1-3答案
《Matlab语言及其应用》形考作业1-3答案形考作业1一、判断题（8小题，每小题2分，共16分）1.Matlab提供了两种创建符号变量和表达式的函数：sym和syms。
（A）判断题(2 分) 2分A.对B.错2.xy12_是合法的Matlab变量名。
（A）判断题(2 分) 2分A.对B.错3.Matlab中的变量必须遵循先定义后使用的原则。
（B）判断题(2 分) 2分A.对B.错4.Matlab的含义为矩阵实验室（Matrix Laboratory）。
（A）判断题(2 分) 2分A.对B.错5.Matlab中注释行的开头是#。
（B）判断题(2 分) 2分A.对B.错6.plot( )是Matlab最基本的绘图函数。
（A）判断题(2 分) 2分A.对B.错7.双纵坐标图形的绘制函数是plotyy（）。
（A）判断题(2 分) 2分B.错8.在Matlab中,二维饼图由函数pie( )来实现。
（A）判断题(2 分) 2分A.对B.错二、单选题（9小题，每小题3分，共27分）9.清空Matlab工作空间内所有变量的指令是（C）。
单选题(3 分) 3分A.clcB.clsC.clearD.clf10.在循环结构中跳出循环，执行循环后面代码的命令为（B）。
单选题(3 分) 3分A.returnB.breakC.continueD.keyboard11.图形窗口的分割函数subplot(2,2,[1 3])可以把一个窗口分割为（B）个子窗口。
单选题(3 分) 3分A.2B.3C.4D.无法确定12.一般情况下，绘图命令每执行一次就刷新当前图形窗口，图形窗口原有图形将不复存在。
若希望在已绘图形上叠加新的图形，可以使用（A）命令。
单选题(3 分) 3分A.hold onB.hold offD.plot13.Matlab中求极限的函数为（B）。
单选题(3 分) 3分A.diff( )B.limit( )C.int( )D.solve( )14.Matlab中求导数的函数为（A）。
北科大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是对应元素相乘。
matlab第三次作业
i 3 10000 的最小 m 值，由命令窗 i 1
m
口(command window)验证运行结果。
值， ， 函数输入值为 n， 返回值为 m 和 s， 由命令窗口(command window)输入[m， sum]=myfun7(20000)验证。
function [m,sum]=myfun7(n) m=0; sum=0;
>> ans = 30.6998 >> 9．使用 switch+case 语句编写如下函数 jijie=myfun9(yue)，当时输入 12，1，2 时输出 winter， 当时输入 3， 4， 5 时输出 spring， 当时输出 autumn。 结果由命令窗口(command window)输入 myfun9(7)和 s=myfun9(11)验证。 function jijie=myfun9(yue) switch(yue) case {3 4 5} jijie='spring'; case {6 7 8} jijie='summer'; case {9 10 11} jijie='autumn'; case {12 1 2} jijie='winter'; end > myfun9(7) ans = summer >> myfun9(11) ans =
贴到 word 中）
12．编写一个 M 脚本文件 ti12.,m，使用 plot 命令绘制 y=sinx+cosx 的结果图, 其中 x=0:0.1:2π ，计算结果，要求曲线为红色线，菱形符数据点形，图名： “cosx+sinx” ； x 坐标轴名： “x 值” ；y 坐标名： “cosx+sinx 值” ；显示分格线。 绘图结果粘贴到下面：
北京科技大学概率论与数理统计上机报告材料3
精彩文档概率论与数理统计第三次上机报告专业：信息与计算科学班级：信计1502（35组）学生姓名：吕瑞杰陈炎睿何芝芝指导教师：张志刚完成时间：2019年5月7日Matlab 概率论与数理统计上机练习（3）五、假设检验【例】（离散型分布检验）某工厂近五年发生了63起事故，按星期几可以分为[9 10 11 8 13 12]，问该厂发生的事故数是有与星期几有关？-------------------------------------------------------------------------------------63 6 0 5 1.6667 11.0705 0.8931 --------------------------------------------------------------------------------------α=），对数学分析I【练习3.1】（基本计算，两个正态总体的假设检验，检验水平0.05（1）求课程中“专业（数学、信计）””的考试人数、平均分、最小值、最大值、极差、标准差、及格人数、及格率、优良人数（大于等于80）、优良率；写出标准差的计算公式。
（2）对“专业（数学、信计）”，检验方差、平均分是否相等。
（3）对“专业（数学、信计）”，检验及格率、优秀率是否相等。
（4）对“全体成绩”的分布进行检验，首先估计期望和方差，画出正态分布的密度函数曲线以及样本密度散点，对假设的正态分布进行检验。
Matlab程序实现：sy=[60 60 63 63 40 69 65 60 72 67 62 78 82 90 69 60 72 76 78 93 69 68 9571 83 60 73 73 60 74 77 71 85 70 89 60 61 77 62 68 60 70 66 84 74 6961 60 86 73 69 74 71 74];se=[50 81 67 65 77 71 76 62 89 65 65 62 62 60 78 81 66 70 80 53 69 66 6148 66 69 61 60 60 85 52 68 60 74 60 62 43 61 60 60 64 70 74 65 73 7960 43 76 66 63 60 60 68 60 60 60 67 74 64];alpha=0.05; %取显著水平为0.05sy1=length(sy);se1=length(se);%人数sy2max=max(sy);sy2min=min(sy);se2max=max(se);se2min=min(se);%最大值，最小值sy3=range(sy);se3=range(se);%极差sy4=mean(sy);se4=mean(se);%平均分sy5=sqrt(sum((sy-sy4).^2)/(sy1));se5=sqrt(sum((se-se4).^2)/(se1));%标准差%sy5=std(sy);se5=std(se);或sy6=length((find(sy>=60)));se6=length((find(se>=60)));%及格人数sy7=length((find(sy>=80)));se7=length((find(se>=80)));%优秀人数sy8=sy6/sy1;se8=se6/se1;%及格率sy9=sy7/sy1;se9=se7/se1;%优秀率fprintf('\t人数\t 平均分\t 最小值\t 最大值 \t极差\t\t标准差\t\t及格人数及格率\t 优秀人数优秀率\n');fprintf('--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n');fprintf('数学 %4d\t%10.4f\t%4d\t\t%4d\t\t%4d\t\t%10.4f\t%4d\t %10.4f\t%4d\t%10.4f\n',sy1,sy4,sy2min,sy2ma x,sy3,sy5,sy6,sy8,sy7,sy9)fprintf('信计 %4d\t%10.4f\t%4d\t\t%4d\t\t%4d\t\t%10.4f\t%4d\t %10.4f\t%4d\t%10.4f\n',se1,se4,se2min,se2ma x,se3,se5,sy6,sy8,se7,se9)fprintf('\n');%方法一fprintf('检验数学和信计的方差是否相等\n');[h1,p1,varci1,stats1]=vartest2(sy,se,alpha,'both');if(h1==0)disp('结果：方差相等');elsedisp('结果：方差不相等');endfprintf('\n');% %方法二% F=sy5^2/se5^2;%统计量F，满足F分布% alpha=0.05; %取显著水平为0.05% Fla1=finv(alpha/2,sy1-1,se1-1);Fla2=finv(1-alpha/2,sy1-1,se1-1);%求F的临界值% if (F>Fla1 && F<Fla2)% MM='数学分析1和数学分析2的方差无显著差异';% else% MM='数学分析1和数学分析2的方差有显著差异';% end% fprintf('检验数学分析1和数学分析2的方差是否相等\n');% fprintf('统计量F的值\t\t\t显著性水平\t\t临界值\t\t\t\t\t检验结果\n');% fprintf(' %.4f\t\t\t\t%.4f\t\t\t%.4f\t\t\t%15s\n',F,alpha,Fla1,MM);% fprintf('\n\n');%方法一fprintf('检验数学和信计的平均分是否相等\n');[h2,p2,muci2,stats2]=ttest2(sy,se,alpha,'both');if(h2==0)disp('结果：平均分相等');elsedisp('结果：平均分不相等');endfprintf('\n');%%方法二% %%%%方法三% sw=((sy1-1)*sy5^2+(se1-1)*se5^2)/(sy1+se1-2);% T=(sy4-se4)/sw/sqrt(1/sy1+1/se1);%统计量T，满T分布% Tla1=tinv(alpha/2,sy1+se1-2);Tla2=tinv(1-alpha/2,sy1+se1-2);%求出T的临界值% if (abs(T)<tinv(1-alpha/2,sy1+se1-2))% XX='数学分析1和数学分析2的平均分无显著差异';% else% XX='数学分析1和数学分析2的平均分有显著差异';% end% fprintf('检验数学分析1和数学分析2的平均分是否相等\n');% fprintf('统计量T的值\t\t\t显著性水平\t\t临界值\t\t\t\t\t检验结果\n ');% fprintf(' %.4f\t\t\t\t%.4f\t\t\t%.4f??%.4f\t\t%15s\n',T,alpha,Tla1,Tla2,XX); % fprintf('\n\n');p=(sy7+se7)/(sy1+se1);U=(sy9-se9)/sqrt((sy9+se9)*p*(1-p));Ua=norminv(1-alpha/2);if(abs(U)>Ua)disp('优秀率无显著差异');elsedisp('优秀率有显著差异');endp=(sy6+se6)/(sy1+se1);U=(sy8-se8)/sqrt((sy8+se8)*p*(1-p));Ua=norminv(1-alpha/2);if(abs(U)>Ua)disp('及格率无显著差异');elsedisp('及格率有显著差异');end[h3,p3,kstat3,critval3]=lillietest(sy,alpha);if(h3==1)disp('数学不是正态分布')elsedisp('数学是正态分布')end[h4,p4,kstat4,critval4]=lillietest(se,alpha);if(h4==1)disp('信计不是正态分布')elsedisp('信计是正态分布')end%%hist(sy)%直方图%[h5,p5,stats5]=chi2gof(sy)%可以检验分布%cdf=[sy,normcdf(sy,sy4,sy5)]%[h5,p5,ksstat,cv5]=kstest(sy,cdf)% a=0:1:100;% a=a';% CDF=[a,cdf(a,sy4,sy5)];% h = kstest(sy,CDF,0.05);S=[65 65 68 81 74 76 68 69 82 77 74 66 73 72 77 60 62 81 66 68 76 60 7480 90 69 60 63 68 67 69 62 60 60 60 67 60 77 67 60 60 60 71 72 60 6661 86 64 60 60 89 73 74 43 40 61 95 69 70 62 66 63 62 78 74 60 50 7662 65 84 70 69 83 73 43 71 70 73 71 69 74 60 61 60 70 74 78 48 93 6461 79 71 53 60 60 52 63 60 61 65 62 78 60 65 60 85 85 89 69 66 60]; [h,p,jbstat,critval]=jbtest(S,alpha);if(h==0)disp('服从正态分布');elsedisp('不服从正态分布');endsavg=mean(S);svar=var(S);x=20:130;y=normpdf(x,savg,sqrt(svar));d=5;a=20:d:130;pdf=hist(S,a)./length(S)./d;plot(x,y,'r');hold onscatter(a,pdf,'filled');hold off输出：>> lx3_1_lrj_41521335人数平均分最小值最大值极差标准差及格人数及格率优秀人数优秀率--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------数学 54 70.6667 40 95 55 10.0885 53 0.9815 9 0.1667信计 60 65.5167 43 89 46 9.2313 53 0.9815 5 0.0833检验数学和信计的方差是否相等结果：方差相等检验数学和信计的平均分是否相等结果：平均分不相等优秀率有显著差异及格率有显著差异数学不是正态分布Warning: P is less than the smallest tabulated value, returning 0.001.> In lillietest (line 206)In lx3_1_lrj_41521335 (line 99)信计不是正态分布服从正态分布六、方差分析【例1】（单因素方差分析）考虑温度对某化工产品得率的影响，选择五种不同温度进行试验，每一温度各做三次试验。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《数学实验》报告
实验名称 MATLAB绘图
学院
专业班级
姓名
学号
2014年 5月
一、【实验目的】
学会用MA TLAB绘制二维、三维图形，并为其标注、添色等。
二、【实验任务】
1.用mesh与surf命令绘制三维曲面z=x^2+3y^2的图像，并使用不同的着色效果及光照效果
2.绘制由函数(x^2)/9+(y^2)/16+(z^2)/4=1形成的立体图，并通过改变观测点获得该图形在各个坐标平
面上的头影
3.画三维曲面z=5-x^2-y^2(-2<=x,y<=2)与平面z=3的交线
三、【实验程序】
1. t=-1:0.1:1;
[x,y]=meshgrid(t);
z=x^2+3*y^2;
subplot(1,2,1),mesh(x,y,z),colormap(bone),light('position',[20,20,5])
subplot(1,2,2),surf(x,y,z),colormap(cool)
2. [xx,yy,zz]=sphere(40);
x=xx*2;y=yy*3;z=zz*4;
subplot(2,2,1),surf(x,y,z);
subplot(2,2,2),surf(x,y,z);view(0,90)
subplot(2,2,3),surf(x,y,z);view(90,0)
subplot(2,2,4),surf(x,y,z);view(0,0)
3. t=-2:0.1:2;[x,y]=meshgrid(t);z1=5-x.^2-y.^2;
subplot(1,3,1),mesh(x,y,z1),title('曲面z1=5-x.^2-y.^2');
z2=3*ones(size(x));
subplot(1,3,2),mesh(x,y,z2),title('平面z=3');
r0=abs(z1-z2)<=1;
zz=r0.*z2;yy=r0.*y;xx=r0.*x;subplot(1,3,3);
subplot(1,3,3),plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'.'),title('交线')
四、【实验结果】
1.
2.
3.
五、【实验总结】
进一步掌握了MA TLAB绘图功能，并学习了对三维图形的添色及改变光源、视角变幻的处理。