matlab实验五

matlab实验五
matlab实验五

实验五 基于matlab 的数据拟合与回归分析

一、 实验类型

验证性实验

二、 实验学时

2学时

三、 实验目的

1、掌握利用MATLAB 中实现单因变量的多元线性回归分析(经典多元线性回归分析)的方法;

2、掌握利用MATLAB 中实现多因变量的多元线性回归分析(多对多线性回归分析)的方法。

3、掌握MATLAB 有关逐步回归的命令。

四、 实验内容及要求 实验内容:

1 学习MATLAB 中有关经典多元线性回归分析的命令;

(1) [b,bint,r,rint,stats] = regress(y,X,alpha)

建立回归分析模型

01122,i i i ip ip i y b b x b x b x ε=+++

++其中()()2E 0,D i i εεσ==,1,2,

,i n =.

注:在该命令中,设计矩阵()

(1)

ij n p X x ?+=(X 的第1列全为1 ) 或响应值向量()1i n y y ?=中

的无效值Nan 将被免忽略。

输入参数至少有2个,alpha 是检验的显著性水平,默认值为0.05。

输出参数至少有1个,b 是回归系数的估计值;bint 是各回归系数的置信度为1-alpha 的置信区间(第1列是它们的下界,第2列是它们的上界);r 是残差,rint 是残差的置信区间;stats 给出一个1*4的向量,依次是判定系数2R ,F 统计量的观测值及检验的p 值,以及误差的

方差的估计值2?σ

. 这里//(1)U p MMS

F Q n p MSE

==-- 是检验回归模型是否显著的检验统计

量,当012:0p H b b b ==

==成立时, ~(,1)F F p n p --.

可以通过分析残差的置信区间来观察原始数据是否存在异常点,若残差的置信区间不包括0,则可认为对应的观测值是异常值,将异常值全部去掉,重新建立的回归模型将更加精确。

可以调用命令rcoplot(r,rint)按观测顺序逐个画出各组观测值对应的残差和残差置信区间,从而观测异常值的位置。

(2) stats = regstats(y,X,model ,whichstats ) % 根据指定模型model 建立回归方程

输入参数:响应变量()1i n y y ?=,设计变量()

ij

n p

X x ?=,model 参数

表5.1 model 参数及说明

包括系数的估计值、因变量的预测值、残差、判定系数等23个选项供用户选择,通过这一界面,matlab将用户所选择的选项对应的回归分析结果显示在matlab工作空间中。

图5.1 regstats弹出的图形用户界面

当带有输出选项stats时,直接返回一个结构体stats,它有24字段,包含了图形界面中的所有选项,此时,不弹出用户界面。若输入中带有第4个参数whichstats,则stats仅包含whichstats 指定的统计量。

表5.2 stats字段名称及说明

(3) 乘方法。

输入参数至少两个,X(n p ?)是设计矩阵,y(1n ?)响应变量的观测值。wfun 给出权函数的名称,tune 给出权函数相应的调节系数,const 控制模型中是否包含常数项。若const=’on’或1,则模型中包含常数项,此时自动在X 的第1列加入1列1,若const=’off’或0,则不包含常数项。

表5.3 权函数与调节系数

其中

r =

,resid 是上一步迭代的残差,/0.6745s MAD =,MAD 是残差

绝对值的中位数,若X 有p 列,计算MAD 时应残差绝对值的前p 个值最小值去掉。h 是最小二乘拟合得到的中心化杠杆值。

输出:至少1个参数,b 是回归模型中的系数向量;stats 返回用于进行模型诊断的结构体变量。其字段名称和意义如下:

表5.4 stats 的字段名称及说明

2 [beta,Sigma,E,CovB,logL]= mvregress(X,Y ,Name,Value) 建立多对多的线性回归模型 输入:至少两个参数。d 维响应变量的观测值矩阵()ij n d Y y ?=,设计矩阵X 可以是矩阵或元胞数组,当1d =时,()

ij

n k

X x ?=,当模型包含常数项时,X 的第1列应全为1;当d>1,

并且Y 的每一个维度有相同的设计矩阵时,()

ij

n k

X x ?=或是一个只包含单个n k ?的设计矩

阵的元胞数组;当 d>1,并且Y 的每一个维度有不同的设计矩阵时,X 将是一个长度为d 的元胞数组,每个元胞都是一个n k ?的设计矩阵。Name 和Value 参数可指定算法与相应的一些选项。其值如表5.5所示。

表5.5 mvregress 函数的name 选项及其值(部分)

的协方差矩阵的估计;E 是残差矩阵;CovB 是beta 的协方差矩阵;,logL 是最后一次迭代对应的对数似然函数的值。

3 学习MATLAB 有关逐步回归的命令。

表5.6 逐步回归命令函数

4 (1) 设某种水泥在凝固时所释放的热量Y(卡/克)与水泥中下列4种化学成分有关

1x 为233CaO Al O ?的成分(%);2x 为23CaO SiO ?的成分(%); 3x 为23234CaO Al O Fe O ??的成分(%);4x 为22CaO SiO ?的成分(%)。

共观测了13组数据(见表5.7)。试求出Y 与1234,,,x x x x 的回归方程,并对该回归方程和各个回归系数进行检验。

表5.7 水泥数据

(2) 设发电量1,工业总产值2与钢材产量1,水泥产量2,机械工业总产值3,棉纱产量4x ,

机制纸产量5x 之间有线性相关关系。现收集了1949到1978年共30年的数据(见表4.2)。试建立12,Y Y 与12345,,,,x x x x x 的线性模型。

表5.8 发电量与经济发展数据

(1)实验过程要求用Matlab 软件完成;

(2)实验结果逐个导入Word 文档,并按问题作出解释; (3)实验报告按照既定格式书写。

五、 实验的重点和难点

实验重点: MA TLAB 中关于多元线性回归模型的实现方法。 实验难点: 无

六、 实验准备

1、实验室电脑要求安装有Matlab软件;

2、上实验课程的学生要对涉及到的统计概念有所了解;

3、要求学生事先查阅并熟悉Matlab的相关命令。

七、实验程序及结果分析

1 学习MATLAB中有关经典多元线性回归分析的命令

(1)

程序如下:

y=[7613.51 7850.91 8381.86 9142.81 10813.6 8631.43 8124.94 9429.79 10230.81 10163.61 9737.56 8561.06 7781.82 7110.97]';

x1=[7666 7704 8148 8571 8679 7704 6471 5870 5289 3815 3335 2927 2758 2591]'; x2=[16.22 16.85 17.93 17.28 17.23 17 19 18.22 16.3 13.37 11.62 10.36 9.83 9.25]'; X=[ones(size(y)) x1.^2 x2.^2 x1 x2 x1.*x2];

[b,bint,r,rint,stats] = regress(y,X)

结果:

结果说明:b是回归系数的估计值;bint是各回归系数的置信度为1-alpha的置信区间(第1列是它们的下界,第2列是它们的上界);r是残差,rint是残差的置信区间;stats给出一个1*4的向量,依R,F统计量的观测值及检验的p值,以及误差的方差的估计值2? .

次是判定系数2

(2)

程序如下:

x=[ones(10,1) (1:10)'];

y=x*[10;1]+normrnd(0,0.1,10,1);

[b,bint]=regress(y,x,0.05)

结果:

(3)

程序如下:

y=[3.73 6.27 5.93 5.77 5.72 5.80 5.87 5.78 5.96 ...

6.03 6.50

7.00 6.80 6.68 7.03 7.67 7.59 6.96 ...

7.17 6.99 7.17 6.99 6.64 6.71 7.01 7.40 7.49 ...

7.75 8.17 8.09 8.11 8.48 8.99];

x=1:33;

scatter(x,y);%画图,离散点

hold on;

p=regress(y',[ones(33,1) x']);%线性拟合

r=robustfit(x,y);%稳健拟合

plot(x,p(1)+p(2)*x,':',x,r(1)+r(2)*x,'r');%画图

legend('原始数据','线性拟合','稳健拟合')%设置图例

结果:

2. 学习MATLAB有关多对多的多元线性回归分析的命令

程序如下:

x=[143 144 145 147 148 150 153 154 155 156 157 158 159 160 161 162]'; X=[ones(16,1),x];

Y=[87 85 88 91 92 90 93 95 98 98 97 95 97 99 100 102]';

[beta,Sigma,E,CovB,logL]= mvregress(X,Y,'covtype','full')

%建立多对多的线性回归模型

结果:

结果说明:beta是回归方程的系数估计矩阵;Sigma是对响应向量的协方差矩阵的估计;E是残差矩阵;CovB是beta的协方差矩阵;,logL是最后一次迭代对应的对数似然函数的值。

3.学习MATLAB有关逐步回归的命令。

程序如下:

y=[11.2 13.4 40.7 5.3 24.8 12.7 20.9 35.7 8.7 ...

9.6 14.5 26.9 15.7 36.2 18.1 28.9 14.9 25.8 21.7 25.7];

%犯罪率y,低收入家庭百分比x1,失业率x2,人口总数x3

x1=[16.5 20.5 26.3 16.5 19.2 16.5 20.2 21.3 17.2 14.3 ...

18.1 23.1 19.1 24.7 18.6 24.9 17.9 22.4 20.2 16.9];

x2=[6.2 6.4 9.3 5.3 7.3 5.9 6.4 7.6 4.9 6.4 6.0 7.4 ...

5.8 8.6

6.5 8.3 6.7 8.6 8.4 6.7];

x3=[587 643 635 692 1248 643 1964 1531 713 749 7895 ...

762 2793 741 625 854 716 921 595 3353];

figure(1),plot(x1,y,'*');

figure(2),plot(x2,y,'*');

figure(3),plot(x3,y,'*');

X1=[x1',x2',x3'];

stepwise(X1,y)

结果:

4.

(1)

程序如下:

a=[7 26 6 60;1 29 15 52;11 56 8 20;

11 31 8 47;7 52 6 33;11 55 9 22;

1 31 2

2 44;1 31 22 44;2 54 18 22;

21 47 4 26;1 40 23 34;11 66 9 12;

10 68 8 12];

x1=a(:,1);x2=a(:,2);

x3=a(:,3);x4=a(:,4);

y=[78.5;74.3;104.3;87.6;95.9;109.2;

72.5;72.5;93.1;115.9;83.8;113.3;109.4]; X=[ones(size(x1)),x1,x2,x3,x4];

[b,bint,r,rint,stats] = regress(y,X)

结果如下:

b =

69.7661

1.4118

0.4621

-0.0399

-0.2226

bint =

-97.7881 237.3202

-0.3563 3.1798

-1.2765 2.2008

-1.8203 1.7406

-1.9175 1.4722

r =

0.4328

1.8950

-2.1031

-1.2393 -0.1932 3.7441 -2.3311 -2.3311 1.1711 0.7144 2.6235 0.5346 -2.9178 rint =

-3.7577 4.6233 -3.1648 6.9548 -5.8987 1.6924 -6.4122 3.9337 -5.3272 4.9408 -1.2481 8.7364 -7.3498 2.6876 -7.3498 2.6876 -3.9131 6.2554 -2.6165 4.0453 -2.0795 7.3264 -4.7525 5.8216 -7.2782 1.4427 stats =

0.9833 117.4117 0.0000 6.5281

说明:b 是线性方程的系数估计值,并且第一值表示常数系数。bint 是系数估计值的置信度为95%的置信区间,r 表示残差,rint 表示各残差的置信区间,stats 给出一个1*4的向量,依次是判定

系数2R ,F 统计量的观测值,检验的p 值以及误差的方差的估计值2?σ

。当P<α时 拒绝H0,回归模型成立。相关系数2R 越大,说明回归方程越显著。 结论:

Y 与1234,,,x x x x 的回归方程为:

123469.7661+1.4118 0.46210.03990.2226y x x x x =+--

(2) 程序如下:

b=[0.9 0.8 0.14 6.63 0.24 1.47 7.31; 1.0 2.1 0.15 7.07 0.46 1.25 7.42; 2.9 6.3 0.33 7.60 1.02 2.05 11.13; 5.0 4.4 0.78 12.88 1.61 2.49 16.08;

8.2 13.3 1.18 15.86 1.63 3.16 22.86; 13.1 16.8 1.56 18.79 1.93 3.87 29.52; 23.8 17.8 2.11 14.63 2.31 4.50 34.54;

34.8 27.8 0.09 19.79 3.32 6.09 41.22;

35.4 22.1 3.58 16.50 4.44 6.78 47.54; 47.0 32.2 7.31 26.22 7.18 10.73 60.00; 62.6 33.2 9.61 28.00 8.77 17.65 78.00; 68.0 55.6 12.85 27.56 9.89 26.84 96.20;

35.3 24.4 6.76 10.95 5.58 24.20 52.37; 31.3 17.9 5.08 10.15 6.03 20.08 37.77; 35.2 24.8 5.54 14.23 7.18 19.28 40.07; 45.3 37.8 7.14 20.38 8.80 22.89 50.36; 49.5 78.8 11.20 26.56 10.45 28.94 65.33;

59.7 101.6 15.89 33.18 12.51 39.05 83.64;

47.8 74.9 10.86 23.90 11.42 39.09 68.16;

17.7 40.2 5.10 17.56 9.03 26.81 41.64; 36.0 73.3 13.14 27.20 8.05 37.19 67.30;

62.0 138.6 25.54 36.28 10.30 54.09 103.57;

97.0 247.0 31.31 41.53 14.18 77.39 135.80;

95.2 270.0 28.79 40.24 15.19 84.02 118.10; 118.4 233.5 28.03 38.20 15.77 88.39 119.62;

99.9 205.0 26.50 31.54 12.29 86.32 112.39; 151.0 288.0 38.61 46.87 17.36 107.94 144.41; 108.0 262.2 31.46 38.62 15.10 102.76 130.66; 162.5 358.6 46.21 52.48 20.48 118.84 175.10; 238.2 454.8 55.86 55.96 26.40 139.30 214.44] x1=b(:,1);x2=b(:,2);x3=b(:,3);x4=b(:,4);x5=b(:,5);

y1=b(:,6);y2=b(:,7);

X=[ones(size(x1)),x1,x2,x3,x4,x5];

[b1,bint1,r1,rint1,stats1] = regress(y1,X)

[b2,bint2,r2,rint2,stats2] = regress(y2,X)

部分结果如下:

b1 =

8.4549

-0.1439

0.1881

1.4499

-0.6775

1.8788

bint1 =

0.1575 16.7523

-0.3315 0.0436

0.0725 0.3036

0.2561 2.6438 -1.2848 -0.0702 0.4690 3.2887 stats1 =

0.9798 232.6942 0.0000 40.3087 b2 = 3.6644 0.3059 -0.1148 1.9173 1.0957 0.9190 bint2 =

-5.0812 12.4101 0.1082 0.5036 -0.2366 0.0070 0.6590 3.1757 0.4556 1.7358 -0.5671 2.4050 stats2 =

0.9862 343.1469 0.0000 44.7810 结论:

1y 与12345,,,,x x x x x 的回归方程为:

1123458.45490.1439 0.1881 1.44990.6775 1.8788y x x x x x =-++-+

由于检验的p 值=0.0000<0.05,因此回归方程是显著的。

2y 与12345,,,,x x x x x 的回归方程为:

2123453.66440.30590.1148 1.9173 1.09570.919y x x x x x =--+++ 由于检验的p 值=0.0000<0.05,因此回归方程是显著的。

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:

matlab实验报告

数学实验报告 班级: 学号: 姓名: 实验序号:1 日期:年 月 日 实验名称:特殊函数与图形 ◆ 问题背景描述:绘图是数学中的一种重要手段,借助图形,可以使抽象的对象得到 明白直观的体现,如函数的性质等。同时,借助直观的图形,使初学者更容易接受新知识,激发学习兴趣。 ◆ 实验目的:本实验通过绘制一些特殊函数的图形,一方面展示这些函数的特点属性, 另一方面,就 Matlab 强大的作图功能作一个简单介绍。 实验原理与数学模型: 1、 球2222x y z R ++= ,x=Rsin φcos θ, y= Rsin φsin θ, z= cos φ, 0≤θ≤2π , 0≤φ≤π 环面 222222222()4(),(cos )cos ,x y z a r a x y x a r φθ+++-=+=- (cos )sin ,sin ,02,02y a r z r φθφφπθπ=-=≤≤≤≤ 2、 平面摆线:2 22 31150,(sin ),(1cos ),0233 x y x a t t y a t t π+-==-=-≤≤ 3、 空间螺线:(圆柱螺线)x=acost , y=asint , z=bt ;(圆锥螺线)22 cos ,sin ,x t t y t t z t === 4、 椭球面sin cos ,sin sin ,cos ,02,0x a y b z c φθφθφθπφπ===≤<≤≤ 双叶双曲面3 tan cos ,tan sin ,sec ,02,22 x a y b z c π φθφθφθπφπ===≤<- << 双曲抛物面2 sec ,tan 2 u x au y bu z θθ=== 实验所用软件及版本:mathematica(3.0) 主要内容(要点): 1、 作出下列三维图形(球、环面) 2、 作出下列的墨西哥帽子 3、 作出球面、椭球面、双叶双曲面,单叶双曲面的图形 4、 试画出田螺上的一根螺线 5、 作出如图的马鞍面

MATLAB第一次实验

一、【实验目的】 了解matlab基础知识和操作 一、【实验任务】 P16 4. 编写函数文件,计算∑k!,并求出当n=20时表达式的值。 P27 2. 矩阵A=[1,2,3;4,5,6;7,8,9],B=[4,6,8;5,5,6;3,2,2],计算A*B,A.*B,并比较两者的区 别。 P27 3. 已知矩阵A=[5,2;9,1],B=[1,2;9,2],做简单的关系运算A>B,A= =B,AB)。 P34 1. 用π/4=1-/3+1/5-1/7+…公式求π的近似值,直到某一项的绝对值小于10-6为止。 二、【实验程序】 P16 第4题 function s=s(n) y=1; sum=0; for i=1:n; y=y*i; sum=sum+y; end sum >> s(20) sum = 2.5613e+18 P27 第2题 A=[1 2 3;4 5 6;7 8 9]; B=[4 6 8;5 5 6;3 2 2]; >>A*B ans = 23 22 26 59 61 74 95 100 122 >> A.*B

ans = 4 12 24 2 25 36 21 16 18 P27 第3题 >> A=[5 2;9 1]; >> B=[1 2;9 2]; >> A>B ans = 1 0 0 0 >> A==B ans = 0 1 1 0 >> A> (A==B)&(A> (A==B)&(A>B) ans = 0 0 0 0

Matlab实验第一次实验答案

实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法; 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容: 1、帮助命令 使用help命令,查找sqrt(开方)函数的使用方法; 解:sqrt Square root Syntax B = sqrt(X) Description B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i

1.0000 1.4142 2、矩阵运算 (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 (2)矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B 解:A=[1 2 3;4 5 6;7 8 9 ]; B=[1 0 0;0 2 0;0 0 3 ]; A\B,A/B (3)矩阵的转置及共轭转置

MATLAB实验报告50059

实验一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-4

3、利用帮助功能查询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命令包含的字符很多,需要分成多行输入,该如何处理?

MATLAB实验报告(5-9)

第五次MATLAB实验报告 一.实验报告 1.学会运用MATLAB分析连续系统的频率特性; 2.学会运用MATLAB进行连续系统的频域分析; 二.实验任务 熟练运用MATLAB进行连续时间LTI系统的频率特性及频率分析,三.实验内容 (1).MATLAB信号处理工具箱提供的freqs函数可直接计算系统的频率响应 的数值解,其语句格式为 H=freqs(b,a,w)。 b,a分别表示H(w)的分子和分母多项式的系统向量; w为系统频率响应的频率范围,其一般形式为w1:p:w2, w1为频率起始值,w2位频率终止值,p为频率取样间隔, H返回的样值可能包括实部和虚部的复数。 因此,如果想要得到系统的幅频特性或相频特性,还需利用abs和angle 函数来分别求得。 (2).连续LTI系统的频域分析法也称傅里叶变换分析法。 四.实验小结 该实验让我们学会了在MATLAB软件上用傅里叶分析法对连续系统进行频率特性和频域的分析,并输出相应的图形。 五.实验作业

(1).实例8-1 实例8-4

六.

第六次MATLAB实验报告 一、实验目的 1.学会运用MATLAB求拉普拉斯变换。 2.学会运用MATLAB求拉普拉斯反变换。 3.学会运用MATLAB求解系统函数的零极点。 4.学会运用MATLAB分析系统函数的零极点分布与其时域特性的关系。 5.学会运用MATLAB分析系统函数的极点分布与系统稳定性的分布。 6.学会运用MATLAB绘制波特图。 二、实验任务 熟练运用MATLAB进行一系列的编程输出结果和输出图形。 三、实验内容 1.拉普拉斯变换: 利用MATLAB的符号数学工具箱中laplace函数来实现拉布拉斯变换,语句格式:L=laplace(f); 2.拉普拉斯反变换: (1).语句格式:f=ilaplace(L); (2).用MATLAB函数residue可以得到复杂有理分式F(s)的部分分式展 开式,语句格式:[r,p,k]=residue(B,A); 例题: format rat;

matlab实验五答案1

实验五 1、编写程序,该程序在同一窗口中绘制函数在 []0,2π之间的正弦曲线和余弦曲线,步长 为200/π,线宽为2个象素,正弦曲线设置为蓝色实线,余弦曲线颜色设置为红色虚线,两条曲线交点处,用红色星号标记并标注sin(x)=cos(x);通过函数方式在生成的图形中添加注释,至少应包括:标题,文本注释,图例和坐标轴标注。%x=linspace(0,2*pi,1000);x=[0:pi/200:2*pi]sinx =sin(x);cosx =cos(x); k=find(abs(sinx-cosx)<1e-2);x1=x(k); plot(x,sinx,'LineWidth',2) hold on ,plot(x,cosx,'r:','LineWidth',2)hold on ,plot(x1,sin(x1),'r*')xlabel('x:(0-2\pi)'); ylabel('y:sin(x)/cos(x)'); title('正弦-余弦曲线'); text(x1+0.1,sin(x1),'sin(x)=cos(x)');legend('sin(x)','cos(x)'); 1 2 34 5 6 7 x:(0-2π) y :s i n (x )/c o s (x ) 正弦-余弦曲线 2、绘制图像:双曲抛物面:22 x y z =- ,1616x -<<,44y -<<,并对绘制的双曲抛物面尝试进行视点控制。 [X,Y]=meshgrid(-16:0.4:16,-4:0.1:4);Z =X.^2/16-Y.^2/4; subplot(1,3,1),plot3(X,Y,Z),view(0,180),title('azimuth =0,elevation =180'); subplot(1,3,2),plot3(X,Y,Z),view(-37.5,-30),title('azimuth =

实验二 MATLAB程序设计 含实验报告

实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台(带有MATLAB7.0以上的软件环境)。 M 文件的编写: 启动MATLAB 后,点击File|New|M-File ,启动MATLAB 的程序编辑及调试器(Editor/Debugger ),编辑以下程序,点击File|Save 保存程序,注意文件名最好用英文字符。点击Debug|Run 运行程序,在命令窗口查看运行结果,程序如有错误则改正 三、 实验内容 1.编写求解方程02=++c bx ax 的根的函数(这个方程不一定为一元二次方程,因 c b a 、、的不同取值而定) ,这里应根据c b a 、、的不同取值分别处理,有输入参数提示,当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 (提示:提示输入使用input 函数) 2.输入一个百分制成绩,要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+,90分~99分为A ,80分~89分为B ,70分~79分为C ,60分~69分为D ,60分以下为E 。 要求:(1)用switch 语句实现。 (2)输入百分制成绩后要判断该成绩的合理性,对不合理的成绩应输出出错信息。 (提示:注意单元矩阵的用法) 3.数论中一个有趣的题目:任意一个正整数,若为偶数,则用2除之,若为奇数,则与3相乘再加上1。重复此过程,最终得到的结果为1。如: 2?1 3?10?5?16?8?4?2?1 6?3?10?5?16?8?4?2?1 运行下面的程序,按程序提示输入n=1,2,3,5,7等数来验证这一结论。 请为关键的Matlab 语句填写上相关注释,说明其含义或功能。 4. 的值,调用该函数后,

matlab第一次实验报告

Matlab第一次实验报告 2012029010010 尹康 1. 编程实现以下数列的图像,用户能输入不同的初始值以及系数。并以x,y为坐标显示图像 x(n+1) = a*x(n)-b*(y(n)-x(n)^2); y(n+1) = b*x(n)+a*(y(n)-x(n)^2) 程序代码: n=input('input the number of pionts:'); a=input('input a:'); b=input('input b:'); x=[]; y=[]; x(1)=input('input x1:'); y(1)=input('input y1:'); %输入点数、初始值以及系数for i=2:n x(i)=a*x(i-1)-b*(y(i-1)-x(i-1)^2); y(i)=a*x(i-1)+b*(y(i-1)-x(i-1)^2); %根据已输入的数据进行迭代end figure;plot(x,y,'linewidth',2) axis equal %横纵坐标等比例 text(x(1),y(1),'1st point') %标记初始点 运行结果:

心得体会及改进:在输入某些数据时,所绘曲线可能是一条折线(如:n=5,a=b=x1=1,y1=2)甚至只有一个点(如:n=5,a=b=x1=y1=1),此时可能出现曲线与坐标轴重合或无法看到点的情况,为了更清晰地展现曲线,可以使线宽适当加宽并标记初始点。 2.编程实现奥运5环图,允许用户输入环的直径。 程序代码: 函数circle: %在指定的圆心坐标处,用指定颜色、宽度的线条绘出指定半径、圆心角的弧 function f=circle(r,x,y,color,linw,alp1,alp2) alp=linspace(alp1,alp2); X=r*cos(alp)+x; Y=r*sin(alp)+y; plot(X,Y,color,'linewidth',linw) end 主程序代码: r=input('input r:');

实验5 Matlab程序设计1

实验5 Matlab 程序设计1 实验目的: 1、 掌握建立和执行M 文件的方法; 2、 掌握实现选择结构的方法; 3、 掌握实现循环结构的方法。 实验内容: 1. 从键盘输入一个4位整数,按如下规则加密后输出。加密规则:每位数字都加上7,然 后用和除以10的余数取代该数字;再把第一位与第三位交换,第二位与第四位交换。 2. 求分段函数的值。 ,x x x x y x x x x x x x ?+- <≠-?=-+ ≤<≠≠??-- ?2226035605231且且及其他 用if 语句实现,分别输出x=-5,-3,0,1,2,2.5,3,5时的y 值。请输入x 的值 x=input('请输入x 的值'); if x<0&x~=-3 y=x*x+x-6; elseif 0<=x&x<5&x~=2&x~=3 y=x*x-5*x+6; else y=x*x-x-1; end y

second 请输入x的值-5 y = 14 second 请输入x的值-3 y = 11 >> second 请输入x的值0 y = 6 >> second 请输入x的值1

y = 2 >> second 请输入x的值2 y = 1 >> second 请输入x的值2.5 y = -0.2500 >> second 请输入x的值3

5 >> second 请输入x的值5 y = 19 >> 3.输入一个百分制成绩,要求输出成绩等级A、B、C、D、E,其中90~100分为A,80~89 分为B,70~79分为C,60~69分为D,60分以下为E。 要求: (1)分别用if语句和swich语句实现。 (2)grade=input('input grade:'); (3)if grade<=100&grade>=90 (4)level='A'; (5)elseif grade<=89&grade>80 (6)level='B'; (7)elseif 70<=grade&grade<=79

昆明理工大学MATLAB实验指导书(第二次实验)

************************ MATLAB上机指导书 ************************ 昆明理工大学机电学院 彭用新 2015年3月

实验三符号计算 一、操作部分:在命令窗口执行命令完成以下运算,记录运算结果。 1.findsym:帮助我们获取系统定义的自变量 f= sym('sin(a*x+b*y)'); findsym(f) 2.numden(获取分子分母), sym2poly,(获取多项式时系数)poly2sym(根据多项式系 数获得符号表达式) [n,d]=numden(sym('x*x+y')+sym('y^2')) p=sym('2*x^3+3*x^2+4'); sym2poly(p) x=[2,3,0,4]; poly2sym(x) 3. collect :合并同类项;expand:展开多项式;horner: 分解成嵌套形式;factor:因式 分解;simplify: 对表达式化简 syms x y; collect(x^2*y+y*x-x^2-2*x) collect((x+y)*(x^2+y^2+1), y) syms x y; expand((x-2)*(x-4)) syms x;horner(x^3-6*x^2+11*x-6) syms x;factor(x^3-6*x^2+11*x-6) syms x;simplify((x^2+5*x+6)/(x+2)) 4. finverse :求得符号函数的反函数。 syms x y; finverse(1/tan(x)) f= x^2+y; finverse(f,y) finverse(f) https://www.360docs.net/doc/2c13030274.html,pose 求符号函数的复合函数 syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g) 6. subs :表达式替换。 syms a b;subs(a+b,a,4)

matlab实验报告

实验一小球做自由落体运动内容:一小球竖直方向做自由落体,并无损做往返运动。程序: theta=0:0.01:2*pi x=cos(theta) y=sin(theta) l=1 v=1 while l<10 for t=1:10 y=y+(-1)^l*v*t plot(x,y,[-1,1],[-56,2],'.') axis equal pause(0.1) end l=l+1 end 结果:

-50 -40 -30 -20 -10 收获:通过运用小球自由落体规律,及(-1)^n 来实现无损往 返运动! 实验二 旋转五角星 内容:一个五角星在圆内匀速旋转 程序:x=[2 2 2 2 2 2] y=[0 4/5*pi 8/5*pi 2/5*pi 6/5*pi 0] y1=2*sin(y) x1=2*cos(y) theta=0:4/5*pi:4*pi

x2=2*cos(theta) y2=2*sin(theta) plot(x,y,x1,y1,x2,y2) axis equal theta1=theta+pi/10 x2=2*cos(theta1) y2=2*sin(theta1) plot(x2,y2) axis equal theta=0:4/5*pi:4*pi for rot=pi/10:pi/10:2*pi x=2*cos(theta+rot) y=2*sin(theta+rot) plot(x,y) pause(0.1) end 结果:

-2 -1.5-1-0.500.51 1.52 -2-1.5-1-0.500.511.5 2 收获:通过theta1=theta+pi/10,我们可以实现五角星在圆内匀速 旋转! 实验三 转动的自行车 内容:一辆自行车在圆内匀速转动 程序:x=-4:0.08:4; y=sqrt(16-x.^2); theta1=-pi/2:0.01*pi:3*pi/2; x3=0.5*cos(theta1); y3=0.5*sin(theta1); theta=-pi/2+0.02*pi for k=1:100

MATLAB实验五 函数文件

MATLAB实验报告 学院:光电学院 班级:073-1 姓名:刘颖 学号:200713503117

实验五 函数文件 1.定义一个函数文件,求给定复数的指数、对数、正弦和余弦,并在命令文件中调用该函数文件。 程序设计: function [e ln s c]=num(x) e=exp(x) ln=log(x) s=sin(x) c=cos(x) end 运行结果: >> num(5i) e = 0.2837 - 0.9589i ln = 1.6094 + 1.5708i s = 0 +74.2032i c = 74.2099 ans = 0.2837 - 0.9589i 2.一物理系统可用下列方程组来表示: ??? ? ??? ???????= ?????? ??? ??? ???????????? ??----g g m m N N a a m m m m 2121212 111001cos 0 0sin 00cos 0 sin 0sin cos θ θθ θθθ 从键盘输入 m 1 、 m 2 和θ的值,求 N a a 121、、和 N 2 的值。其中g 取9.8,输入θ时以角度为单位。 程序设计: 函数文件in.m: function [a1,a2,N1,N2]=in(m1,m2,t) g=9.8; A=[m1*cos(t) -m1 -sin(t) 0;m1*sin(t) 0 cos(t) 0;0 m2 -sin(t) 0;0 0 -cos(t) 1]; C=[0;m1*g;0;m2*g]; B=inv(A)*C; a1=B(1); a2=B(2); N1=B(3); N2=B(4); end 调用in.m 的命令文件: >> m1=1;m2=2;t=30*pi/180; >> [a1,a2,N1,N2]=in(m1,m2,t) 运行结果: a1 = 6.5333 a2 = 1.8860 N1 = 7.5440 N2 = 26.1333 4.设 f(x)= 01 .01 1 .01 ) 3() 2(4 2 +++--x x , 编写一个MATLAB 函数文件fx.m ,使得调用f(x)时,x 可用矩阵代入,得出的f(x)为同阶矩阵。 程序设计: 函数文件fx.m: function A=fx(x) A=1./((x-2).^2+0.1)+1./(((x-3).^4)+0.01) end 调用fx.m 的命令文件: >> A=fx([1 2;2 3;4 3]) 运行结果: A = 0.9716 10.9901 10.9901 100.9091 1.2340 100.9091 5.已知y= ) 20()30() 40(f f f + (1)当f(n)=n+10ln(n 2+5)时,求y 的值。

MATLAB第二次上机实验报告

电子科技大学电子工程学院标准实验报告(实验)课程名称MATLAB与数值分析 学生姓名: 学号: 指导教师:

一、实验名称 实验二 线性方程组求解和函数的数值逼近 二、实验目的 通过上机实验,使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。 实验涉及的核心知识点:病态方程求解、矩阵分解和方程组求解、Lagrange 插值。 实验重点与难点:算法设计和MATLAB 编程 三、实验内容 1. 对高阶多项式 ()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 2. 对2 20n =,生成对应的Hilbert 矩阵,计算矩阵的条件数;通过先确定解获得常向量 b 的方法,确定方程组 ()n H x b = 最后,用矩阵分解方法求解方程组,并分析计算结果。 3. 对函数 ()2 1 125f x x = + []1,1x ∈- 的Chebyshev 点 ()()21cos 21k k x n π ?? -= ? ?+? ? ,1,2,,1k n =+ 编程进行Lagrange 插值,并分析插值结果。 四、实验数据及结果分析 1. 对高阶多项式

()()() ()()20 1 1220k p x x x x x k ==---=-∏ 编程求下面方程的解 ()190p x x ε+= 并绘图演示方程的解与扰动量ε的关系。 p=[1,-1]; for i=2:20 n=[1,-i]; p=conv(p,n); % 求多项式乘积 end m=zeros(1,21); % m 的最高次幂为20,有21项 hold on x=1:20; d=[-1,0,0.1,0.5,1]; for i=1:5 delt=d(i); m(2)=delt; y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); end title('方程p(x)=0的解与扰动量delt 的关系') legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1') 2468101214161820 010 20 30 40 50 60 方程p(x)=0的解与扰动量delt 的关系 delt=-1delt=0delt=0.1delt=0.5delt=1

参考答案Matlab实验报告

实验一 Matlab基础知识 一、实验目的: 1.熟悉启动和退出Matlab的方法。 2.熟悉Matlab命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握Matlab各种表达式的书写规则以及常用函数的使 用。 二、实验内容: 1.求[100,999]之间能被21整除的数的个数。(rem) 2.建立一个字符串向量,删除其中的大写字母。(find) 3.输入矩阵,并找出其中大于或等于5的元素。(find) 4.不采用循环的形式求出和式 63 1 2i i= ∑ 的数值解。(sum) 三、实验步骤: ●求[100,199]之间能被21整除的数的个数。(rem) 1.开始→程序→Matlab 2.输入命令: ?m=100:999; ?p=rem(m,21); ?q=sum(p==0) ans=43 ●建立一个字符串向量,删除其中的大写字母。(find) 1.输入命令:

?k=input('’,’s’); Eie48458DHUEI4778 ?f=find(k>=’A’&k<=’Z’); f=9 10 11 12 13 ?k(f)=[ ] K=eie484584778 ●输入矩阵,并找出其中大于或等于5的元素。(find) 1.输入命令: ?h=[4 8 10;3 6 9; 5 7 3]; ?[i,j]=find(h>=5) i=3 j=1 1 2 2 2 3 2 1 3 2 3 ●不采用循环的形式求出和式的数值解。(sum) 1.输入命令: ?w=1:63; ?q=sum(2.^w) q=1.8447e+019

实验二 Matlab 基本程序 一、 实验目的: 1. 熟悉Matlab 的环境与工作空间。 2. 熟悉M 文件与M 函数的编写与应用。 3. 熟悉Matlab 的控制语句。 4. 掌握if,switch,for 等语句的使用。 二、 实验内容: 1. 根据y=1+1/3+1/5+……+1/(2n-1),编程求:y<5时最大n 值以及对应的y 值。 2. 编程完成,对输入的函数的百分制成绩进行等绩转换,90~100为优,80~89为良,70~79为中,60~69为及格。 3. 编写M 函数文件表示函数 ,并分别求x=12和56时的函数值。 4. 编程求分段函数 2226;03 56;0532 1;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它,并求输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。 三、 实验步骤: 根据y=1+1/3+1/5+……+1/(2n-1),编程求:y<5时最大n 值以及对应的y 值。 1. 打开Matlab ,新建M 文件 2. 输入命令: 51022-+x

matlab第五章思考与实验

第五章思考与实验 一思考题 1在matlab数据处理中,数据序列是如何表示的? 答:可以让矩阵的每列或每行代表不同的被测变量,相应的行或列的元素代表被测向量的观测值,这样就很容易通过对矩阵元素的访问进行数据的处理和分析。 2函数sum与cumsum有何区别? 答:sum(A):如果A是一个向量,则返回向量各元素的和。如果A是一个矩阵,则返回一个行向量,其第i个元素是A的第i列的元素和。 sum(A,dim):当dim为1时,该函数等同于sum(A);当dim为2时,返回一个列向量,其第i 个元素是A的第i行的个元素之和。 cumsum(A):如果A是一个向量,则返回向量的累加和。如果A是一个矩阵,则返回一个矩阵,其第i列是A的第i列的累加和向量。 cumsum(A,dim):当dim为1时,该函数等同于cumsum(A);当dim为2时,返回一个矩阵,其第行是的第i行的累加和向量。 3利用randn函数生成符合正态分布的10*5随机矩阵A,写出完成下列操作的命令。 (1)A各列元素的均值和标准方差。 答:>> a=randn(10,5) a = 1.3005 -1.3165 -0.5492 0.0601 -1.5809 0.2691 -0.6103 1.0014 0.0510 -1.0978 -0.1551 0.4468 0.7482 0.7101 0.2103 0.0342 -2.4119 0.7242 -0.0953 -0.4456 0.9913 1.6895 0.7019 -1.2037 -0.2694 -1.3618 -0.9683 0.8053 0.2633 -1.1721 0.9792 1.4889 -0.1651 -0.1608 0.1277 -0.8863 -1.0166 -1.0022 -0.1653 1.0411 -0.3562 -2.1758 0.5817 1.1027 1.2827 -0.1428 -0.5737 2.1113 0.3762 0.3490 >> mean(a) ans = 0.0672 -0.5448 0.4957 0.0938 -0.1555 >> std(a) ans =

(完整版)Matlab实验5选择结构程序结构

实验五、选择与循环结构 一、实验目的: 1、 掌握建立和执行M 文件的方法。 2、 掌握利用if 语句实现选择结构的方法。 3、 掌握利用switch 语句实现多分支选择结构的方法。 4、 掌握try 语句的使用。 5、 掌握利用for 语句实现循环结构的方法。 6、 掌握利用while 语句实现循环结构的方法。 7、 熟悉利用向量运算来代替循环的操作方法。 二、实验内容: 1、 列分段函数的值。 ?? ???--≠≠<≤+--≠<-+=其他且且,632,100,6530,6222x x x x x x x x x x x y 要求: (1) 用if 语句实现,分别输出x =-0.5,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。 提示:x 的值从键盘输入,可以是向量。 %homework_5_1_1.m x=input('请输入x 的值:x='); if (x<0 & x~=-3) y= x.*x + x - 6 elseif (x>=0 & x<10 & x~=2 & x~=3) y=x.*x-5.*x+6 else y=x.*x-x-6 end >> homework_5_1 请输入x 的值:x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0] y = -5.2500 6.0000 -6.0000 -4.0000 -2.2500 0 14.0000 (2) 用逻辑表达式实现上述函数。 %homework_5_1_2.m x=input('请输入x 的值:x=') y=(x<0 & x~=-3).*(x.*x+x-6)... +(x>=0 & x<10 &x~=2 &x~=3).*(x.*x-5.*x+6)... +(x>=10 | x==-3 | x==3 | x==2).*(x.*x-x-6) >> homework_5_1_2 请输入x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0] x = -0.5000 -3.0000 1.0000 2.0000 2.5000 3.0000 5.0000 y = -6.2500 6.0000 2.0000 -4.0000 -0.2500 0 6.0000

第二次数学实验报告Matlab 二维曲线绘图

《数学实验》报告实验名称 Matlab 二维曲线绘图 2011年 5月

一、【实验目的】 学习Matlab 绘图的运用,学会制作二维曲线,三维图形的绘画。 二、【实验任务】 P79 第3,5,9题。 1,在同一图形窗口画三个子图…… 2,绘制圆锥螺线的图像并加各种标注…… 3,画三维曲面z=5-x^2-y^2与平面z=3的交线。 三、【实验程序】 1. >> clear >> x=-pi:pi/50:4*pi; y1=x.*cos(x); y2=x.*tan(1./x).*sin(x.^3); y3=exp(1./x).*sin(x); subplot(3,1,1) plot(x,y1,'r*'),grid on title('y1=xcosx') xlabel('x轴'),ylabel('y轴') axis([-pi pi -pi pi]) gtext('y1=xcosx'),legend('y1=xcosx') subplot(3,1,2),plot(x,y2,'b'),grid on title('y=xtan(1/x)sin(x^3)') gtext('y=xtan(1/x)sin(x^3)') legend('y=xtan(1/x)sin(x^3)') axis([pi 4*pi -2 2]) subplot(3,1,3),plot(x,y3,'y'),grid on title('y=exp(1/x)sinx') xlabel('x轴'),ylabel('y轴') gtext('y=exp(1/x)sinx') legend('y=exp(1/x)sinx') axis([1 8 -3 3]) 2. >> clear >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t); y=t.*sin(pi/6.*t); z=2.*t; plot3(x,y,z) title('圆锥螺线') xlabel('x轴'),ylabel('y轴'),zlabel('z轴') >> t=0:pi/50:20*pi; x=t.*cos(pi/6.*t);

matlab实验报告

Matlab实验报告 实验二图像处理 一、实验目的 (1)通过应用MA TLAB语言编程实现对图像的处理,进一步熟悉MATLAB软件的编程及应用; (2)通过实验进一步掌握图像处理的基本技术和方法。 二、实验内容及代码 ㈠.应用MA TLAB语言编写显示一幅灰度图像、二值图像、索引图像及彩色图像的程序,并进行相互之间的转换 首先,在matlab页面中的current directory下打开存放图像的文件夹。 1.显示各种图像 ⑴显示彩色图像: ①代码:>> mousetif=imread('tif.TIF'); >> image(mousetif) 显示截图: ②代码:>> mousetif=imread('tif.TIF'); >> imshow(mousetif) 显示截图:

③代码:mousetif=imread('tif.TIF'); subimage(mousetif) 显示截图: 显示截图:

⑵显示二值图像 ①代码:>> I=imread('单色bmp.bmp'); >> imagesc(I,[0 2]) 显示截图: ②代码:>> I=imread('单色bmp.bmp');

>> imshow(I,2) 显示截图: ③代码:>> I=imread('单色bmp.bmp'); >> subimage(I) 显示截图:

⑶显示灰度图像 ①代码:>> I1=imread('256bmp.bmp'); >> imagesc(I1,[0,256]) 显示截图: 代码:>> I1=imread('256bmp.bmp'); >> colormap(gray); >> subplot(1,2,1); >> imagesc(I1,[0,256]); >> title('灰度级为[0 256]的mouse.bmp图'); >> subplot(1,2,2); >> imagesc(I1,[0,64]); >> colormap(gray); >> title('灰度级为[0 64]的mouse.bmp图'); 显示截图:

matlab实验报告

MATLAB 数学实验报告 指导老师: 班级: 小组成员: 时间:201_/_/_

Matlab 第二次实验报告 小组成员: 1 题目:实验四,MATLAB 选择结构与应用实验 目的:掌握if 选择结构与程序流程控制,重点掌握break,return , pause语句的应用。 问题:问题1:验证“哥德巴赫猜想” ,即:任何一个正偶数(n>=6)均可表示为两个质数的和。要求编制一个函数程序,输入一个正偶数,返回两个质数的和。 问题分析:由用户输入一个大于6 的偶数,由input 语句实现。由if 判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数,再向屏幕输出两个质数即可。 编程:function [z1,z2]=gede(n); n=input('please input n')

if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n- i) break end

end end 结果分析 如上图,用户输入了大于6的偶数返回两个质数5和31,通过 不断试验,即可验证哥德巴赫猜想。 纪录:if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰,更快的解决问题。 2题目:实验四,MATLAB选择结构与应用实验 目的:用matlab联系生活实际,解决一些生活中常见的实际问 题。