matlab综合大作业(附详细答案)

《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.6724

1.2507 0.9247 -0.1766 1.1186

2.4286

1.5754 1.6546 5.3664 0.8087 4.2471

-1.2929 1.3493 0.7272 -0.6647 -0.3836

2)将矩阵A按列拉长得到矩阵B；

实验程序：B=A(:)

实验结果：

B =

0.1349

-2.3312

1.2507

1.5754

-1.2929

3.3818

3.3783

0.9247

1.6546

1.3493

0.6266

2.4516

-0.1766

5.3664

0.7272

1.2279

3.1335

1.1186

0.8087

-0.6647

1.5888

-1.6724

2.4286

4.2471

-0.3836

3)提取矩阵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.1335

0.9247 1.1186

4)寻找矩阵A中大于0的元素；]

实验程序：G=A(find(A>0))

实验结果：

G =

0.1349

1.2507

1.5754

3.3818

3.3783

0.9247

1.6546

1.3493

0.6266

2.4516

5.3664

0.7272

1.2279

3.1335

1.1186

0.8087

1.5888

2.4286

4.2471

5)求矩阵A的转置矩阵D；

实验程序：D=A'

实验结果：

D =

0.1349 -2.3312 1.2507 1.5754 -1.2929

3.3818 3.3783 0.9247 1.6546 1.3493

0.6266 2.4516 -0.1766 5.3664 0.7272

1.2279 3.1335 1.1186 0.8087 -0.6647

1.5888 -1.6724

2.4286 4.2471 -0.3836

6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；

实验程序：E=flipud(fliplr(A))

实验结果：

E =

-0.3836 -0.6647 0.7272 1.3493 -1.2929

4.2471 0.8087

5.3664 1.6546 1.5754

2.4286 1.1186 -0.1766 0.9247 1.2507

-1.6724 3.1335 2.4516 3.3783 -2.3312

1.5888 1.2279 0.6266 3.3818 0.1349

7)删除矩阵A的第2列和第4列得到矩阵F；

实验程序：F=A;

F(:,[2,4])=[]

实验结果：

F =

0.1349 0.6266 1.5888

-2.3312 2.4516 -1.6724

1.2507 -0.1766

2.4286

1.5754 5.3664 4.2471

-1.2929 0.7272 -0.3836

8)求矩阵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 0

0 -0.2877 + 3.4850i 0 0 0

0 0 -0.2877 - 3.4850i 0 0

0 0 0 0.5915 0

0 0 0 0 -2.3024

9)求矩阵A的每一列的和值；

实验程序：lieSUM=sum(A)

实验结果：

lieSUM =-0.6632 10.6888 8.9951 5.6240 6.2087

10)求矩阵A的每一列的平均值；

实验程序：average=mean(A)

实验结果：

average =-0.1326 2.1378 1.7990 1.1248 1.2417

2．符号计算（10分，每小题5分）：

1) 求方程组20,0uy vz w y z w ++=++=关于,y z 的解；

实验程序：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 dy

y 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) 二维图形绘制：绘制方程22

22

125x 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 4

y 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 滤波器恢复出原始信号。（20分）

实验程序：

t=0:0.001:0.2;

x1=sin(2*pi*80*t)+2*sin(2*pi*150*t);

x2=randn(size(t)); %添加随机噪声； x=x1+x2;

A=[1 -1.143 0.4128]; %滤波器分母系数；

B=[0.06745 0.1348 0.06745]; %滤波器分子系数；

y=filter(B, A, x); plot(t, x, ':' , t, y, '-');

legend('带有噪声的信号波形','FIR滤波后信号波形');

t=0:0.2:200;

x=sin(2*pi*80*t)+2*sin(2*pi*150*t)+randn(size(t));

n=3;

f=9;

[b,g]=sgolay(n,f);

halfwin=((f+1)/2)-1;

for n=(f+1)/2:996-(f+1)/2;

SG0(n)=dot(g(:,1),x(n-halfwin:n+halfwin));

end

subplot(2,1,1);plot(x(1:length(SG0)));

subplot(2,1,2);plot(SG0);

实验结果：

5. 人体心电图测量信号在测量的过程中经常会受到工业高频干扰，所以必须经过低通滤波处理后，才能判断心脏功能的有用信息。下面是一组实际心电图信号采样的样本x(n)，其中存在高频干扰。试在实验中，通过MATLAB程序，以x(n)作为输入序列，滤出其中的干

扰成分。x(n) = {-4，-2，0，-4，-6，-4，-2，-4，-6，-6，-4，-4，-6，-6，-2，6，12，8，0，-16，-38，-60，-84，-90，-66，-32，-4，-2，8，12，12，10，6，6，4，0，0，0，0，0，-2，-2，0，0，-2，-2，-2，-2，0}。（20分，每小题4分））

1)绘制原数据图形；

2)设计巴特沃斯低通滤波器并绘制出其幅频响应曲线；

3)用设计的滤波器对原数据进行滤波；绘制滤波后的数据图；

4)绘制原数据功率谱图；

5)绘制滤波后的数据功率谱图。

解：（1）实验程序：

xn=[-4 -2 0 -4 -6 -4 -2 -4 -6 -6 -4 -4 -6 -6 -2 6 12 8 0 -16 -38 -60 -84 -90 -66 -32 -4 -2 8 12 12 10 6 6 4 0 0 0 0 0 -2 -2 0 0 -2 -2 -2 -2 0];

N=49;

n=0:N-1;

stem(n,xn); %绘制原数据图形

title('原数据图形');

实验结果：

原数据图形

（2）实验程序：

N_filter=10;fs_filter=5000; %20阶滤波器，采样频率5000Hz

wn_filter=1000/(fs_filter/2); % 截止频率1000Hz

[b,a]=butter(N_filter,wn_filter);

freqz(b,a); %画出频率特性

title('低通滤波器频率特性');

实验结果：

Normalized Frequency (?π rad/sample)

P h a s e (d e g r e e s )

Normalized Frequency (?π rad/sample)

M a g n i t u d e (d B )

低通滤波器频率特性

（3）实验程序：

yn=filter(b,a,xn); %滤波 stem(n,yn);

title('滤波后的数据图');

实验结果：

滤波后的数据图

（4）实验程序：

Nfft=1024;

periodogram(xn,window,nfft,fs_filter); %绘制原数据功率谱图

实验结果：

Frequency (kHz)

P o w e r /f r e q u e n c y (d B /H z )

Periodogram Power Spectral Density Estimate

（5）实验程序：periodogram(yn,window,nfft,fs_filter);

%绘制滤波后数据功率谱图

实验结果：

Frequency (kHz)

P o w e r /f r e q u e n c y (d B /H z )

Periodogram Power Spectral Density Estimate

6.已知x(n)=[2 1 0 1]，计算如下表达式：（10分）

1) 计算()x n 的6点DFT 结果()X k ;

实验程序：xn=[2 1 0 1];

Xk=fft(xn,6)

实验结果：

Xk=4.0000 1.5000-0.8660i 2.5000-0.8660i 0 2.5000+0.8660i 1.5000+0.8660i

2) 已知()()36k Y k W X k =，求()()y n IDFT x n =????；

实验程序： 实验结果：

3) 已知()()Re W k X k =????，求()()w n IDFT X k =????；

实验程序：

实验结果：

4) 已知()()2Q k X k =，求()q n ；

实验程序： 实验结果：

5) 已知()z k =()2X k ，求()z n ；

实验程序： 实验结果：

MATLAB基础训练作业(含答案)

实验一 MATLAB 工作环境熟悉及简单命令的执行 一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。 二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单 （Menu ）、工具条 （Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、启动MATLAB ，熟悉MATLAB 的桌面。 2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。 （1）（365-52?2-70）÷3 （2）>>area=pi*2.5^2 （3）已知x=3，y=4，在MATLAB 中求z ： ()2 3 2 y x y x z -= （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。 m1=????? ???? ???115 14 4 12679810115 133216 执行以下命令 >>m1( 2 , 3 ) >>m1( 11 ) >>m1( : , 3 ) >>m1( 2 : 3 , 1 : 3 ) >>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) （5）执行命令>>help abs 查看函数abs 的用法及用途，计算abs( 3 + 4i ) （6）执行命令 >>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y) （6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。 四、思考题 1、以下变量名是否合法？为什么？ （1）x2 （2）3col （3）_row （4）for 2、求以下变量的值，并在MATLAB 中验证。 （1）a = 1 : 2 : 5 ; （2）b = [ a' a' a' ] ; （3）c = a + b ( 2 , : )

MATLAB大作业

选 题 说 明 本人选做第2、4、5、9、11、12、13、14、16、19、24 题。 作业内容题目2：问题描述：在[0 ， 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先写出x的取值范围，再用plot函数画出y的图像。 （2）软件说明及源代码 >> x = 0:pi/100:2.*pi; y=cos(5*x).*sin(2*x); >> plot(x,y) （3）实验结果 题目4：问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c

(2)求f=0的解 （1）问题分析 这是符号计算问题，首先要确定符号变量，然后创建符号函数，最后利用subs函数求解特值。 （2）软件说明及源代码 >> syms a b c x f; f=a*x^2+b*x+c; subs(f,0) （3）实验结果 ans = c 题目5：问题描述：求积分 （1）问题分析 这是符号计算的积分求解问题，首先需要确定符号变量，然后利用int函数计算积分。 （2）软件说明及源代码 >> syms x y; y=sqrt(1-2*sin(2*x)); >> int(y,x,0,pi/2) （3）实验结果 ans = ellipticE(-pi/4, 4)*1i - ellipticE(pi/4, 4)*1i - ellipticE(-pi/6, 4)*2i + ellipticE(pi/6, 4)*2i 题目9：问题描述：按水平和竖直方向分别合并下述两个矩阵：

（1）问题分析 这是考查矩阵的基本操作，首先定义矩阵，然后合并矩阵。 （2）软件说明及源代码 >> A=[1,0,0;1,1,0;0,0,1]; B=[2,3,4;5,6,7;8,9,10]; >> a=[A,B],b=[A;B] （3）实验结果 a = 1 0 0 2 3 4 1 1 0 5 6 7 0 0 1 8 9 10 b = 1 0 0 1 1 0 0 0 1 2 3 4 5 6 7 8 9 10 题目11：问题描述：计算z＝yx2+3y2x+2y3的和： （1）问题分析 这是符号计算问题，首先确定符号变量，然后构造函数，最后利用diff函数进行求导。 （2）软件说明及源代码 >> syms x y z; >> z=y*x^2+3*y^2*x+2*y^3; >> diff(z,y,1),diff(diff(z,y,1),x,1) （3）实验结果 ans = x^2 + 6*x*y + 6*y^2

Matlab课后习题解答

P16 Q2: 计算表达式()2 tan arccos x x -在0.25x =和0.78x π =时的函数值。 function y=jie(x) y=tan(-x.^2)*acos(x); >> jie(0.25) ans = -0.0825 >> jie(0.78*pi) ans = 0 + 0.4418i Q3：编写M 命令文件，求5010 2 1 1 1k k k k ==+ ∑∑ 的值。 a=0;b=0; for i=1:50 a=a+i*i; end for j=1:10 b=b+1/j; end c=a+b; >> c c = 4.2928e+004 P27 Q2：矩阵1234567 8 9A ????=??????,4 685563 2 2B ?? ?? =? ????? ,计算A B *,.A B *,并比较两者的区别。 >> 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 20 25 36 21 16 18 A*B表示A与B两矩阵相乘。 A.*B表示A与B对应元素相乘。 P34 Q2：编写一个转换成绩等级的程序，其中成绩等级转换标准为：考试分数在[] 90,100显示为优秀；分数在[) 0,60的 60,80的显示为及格；分数在[) 80,90的显示为良好；分数在[) 显示为不及格。 if x>=90 disp('优秀'); elseif x>=80 disp('良好'); elseif x>=60 disp('及格'); else disp('不及格'); end >> x=85 x = 85 良好

MATLAB第一章作业答案

第一章 M A T L A B 概况与基本操作 1.选择题（每题2分，共20分）： (1)最初的MATLAB 核心程序是采用D 语言编写的。 (2)即将于2011年9月发布的MATLAB 新版本的编号为C 。 2011Ra 2011Rb R2011b R2011a (3)在默认设置中，MATLAB 中的注释语句显示的颜色是B 。 A.黑色 B.绿色 C.红色 D.蓝色 (4)如果要以科学计数法显示15位有效数字，使用的命令是A 。 long e long long g long d (5)在命令窗口新建变量a 、b ，如果只查看变量a 的详细信息，使用的命令为A 。 a a (6)如果要清除工作空间的所有变量，使用的命令为C 。 all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上B 。 A.冒号 B.分号 C.空格 D.逗号 (8)如果要重新执行以前输入的命令，可以使用D 键。 A.下箭头↓ B.右箭头→ C.左箭头← D.上箭头↑ (9)如果要查询函数det 的功能和用法，并显示在命令窗口，应使用命令C 。 D.三者均可 (10)如果要启动Notebook 文档，下列D 操作是可行的。 A.在命令窗口输入notebook 命令 B.在命令窗口输入notebook filename 命令 C.在Word 中启动M-book 文档 D.三者均可 2.填空题（每空1分，共20分）： (1)MATLAB 是matrix 和laboratory 两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler 和Jack Little 。 (2)在MATLAB 的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为褐红色，注释显示的字体为绿色，错误信息显示的字体为红色。 (3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format loose ，各行之间不添加空行的命令为format compact 。 (4)在MATLAB 中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组每行各个元素，逗号的作用是分隔数组每行各个元素或函数的各个输入参数，分号的作用是作为不显示命令结果的命令行的结尾或分隔数组各列，冒号的作用是生成一维数组或表示数组全部元素，百分号的作用是引导一行注释，…的作用是连接相邻两行，感叹号的作用是调用操作系统命令。 3.程序设计题（每题10分，共40分） (1)以25m/s 的初速度向正上方投球（g=s 2 ），计算到达最高点的时间tp 以及球从出发点到最高点的距离hp 。 解：根据物理学知识，物体上抛运动的速度与经过的时间之间的关系为0p p v v gt =-，因此所需要的时间为0p p v v t g -= 。而到达最高点时的速度0p v =，因此可根据此公式求出tp ： v0=25;g=;vp=0; tp=(v0-vp)/g tp =

MATLAB大作业

MATLAB大作业 作业要求： （1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。 （2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。打印稿不要求MATLAB程序，但电子稿要包含MATLAB 程序。 （3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。 （4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。 （5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。 第一类：绘制图形。（B级） 问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据 斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。 问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢 尔宾斯基在1915年提出,它是一种典型的自相似集。其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。

问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮 亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。这方面的资料很多（如https://www.wendangku.net/doc/f512352144.html,/content/16/0103/14/5315_525141100.shtml），请分析构图原理并用MATLAB实现。 问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。 点数 1 2 3 4 5 6 出现次数166 150 164 162 184 174 问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。 提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。 第二类：插值与拟合。（B级） 问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未 （1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。 （2）绘制插值图形并标注样本点。 问题二：估算矩形平板各个位置的温度。已知平板长为5m，宽为3m，平板上3×5栅格 点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。 （1）分别使用最近点插值、线性插值和三次样条插值进行计算。 （2）用杆图标注样本点。 （3）绘制平板温度分布图。 问题三：自行车道的设计。对9条道路上的自行车道宽度以及自行车与过往机动车之间 （1）对数据进行线性拟合。 （2）绘制拟合曲线和样本点。 （3）如果自行车与过往机动车之间安全距离的最小距离是1.8m，试计算相应的自行车道宽度的最小值。 问题四：在水资源工程学中，水库的大小与为了蓄水而拦截的河道中的水流速度密切相 关。对于某些河流来说，这种长时间的历史水流记录很难获得。然而通常容易得到过去若干年间关于降水量的气象资料。鉴于此，推导出流速与降水量之间的关系式往往特别有用。只

matlab课后习题及答案详解

第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？

Matlab大作业

Matlab 大作业 （组内成员：彭超杰、南彦东、江明伟） 一、研究模型 （电车）通过控制油门（保持一定角度）来调节电动机能输出稳定的转速，从而控制车速稳定。 数学依据说明如下： 由图可知存在以下关系：a d a a u w k R i dt di L =++ (w k e d d =) L M M dt dw J -= a m i k M = L a m M i k dt dw J -=

k为反电势常数，m k为电动机电磁力矩常数，这里忽略阻尼力矩。d

二、数学模型 再看整个研究对象，示意图以课本为依据，不同点是这里将数控的进给运动，转换为汽车行驶所需要的扭矩。（这里不说明扭矩的具体产生过程，仅仅说明输出车轮旋转的角速度w ） 对照课本不同，() s θ变为()s N ，1 221z z w w =，1w 为电动机的转速，2w 为轮胎的转速，1z 为电动机的光轴齿轮的齿数，2z 为与轮胎相连光轴的 齿轮齿数。 )(*10110w x w k x ==，1 21z z k = ()c a m m d b a m x K K K k s k k JRs JLs K K K k s G i 1231+++= () c a m m d M K K K k s k k JRs JLs R Ls K s G L 1231)(++++-= 同理，忽略电枢绕组的电感L ，简化系统传递函数方框图如下

()JR K K K k JR s k k s JR K K K k s G c a m m d b a m x i 121++= ()JR K K K k JR s k k s K K K K k s k k Rs R K s G c a m m d c a m m d M L 121121++-=++-=

matlab课后习题解答第二章doc

第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象？ 3/7+0.1; sym(3/7+0.1); sym('3/7+0.1'); vpa(sym(3/7+0.1)) 〖目的〗 ●不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+0.1 c2=sym(3/7+0.1) c3=sym('3/7+0.1') c4=vpa(sym(3/7+0.1)) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = 0.5286 c2 = 37/70 c3 = 0.52857142857142857142857142857143 c4 = 0.52857142857142857142857142857143 Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 ●理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1) ans = a

symvar(sym('z*exp(j*th)'),1) ans = z 3 求以下两个方程的解 （1）试写出求三阶方程05.443 =-x 正实根的程序。注意：只要正实根，不要出现其他根。 （2）试求二阶方程022=+-a ax x 在0>a 时的根。 〖目的〗 ● 体验变量限定假设的影响 〖解答〗 （1）求三阶方程05.443 =-x 正实根 reset(symengine) %确保下面操作不受前面指令运作的影响 syms x positive solve(x^3-44.5) ans = (2^(2/3)*89^(1/3))/2 （2）求五阶方程02 2 =+-a ax x 的实根 syms a positive %注意：关于x 的假设没有去除 solve(x^2-a*x+a^2) Warning: Explicit solution could not be found. > In solve at 83 ans = [ empty sym ] syms x clear syms a positive solve(x^2-a*x+a^2) ans = a/2 + (3^(1/2)*a*i)/2 a/2 - (3^(1/2)*a*i)/2 4 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。 a =@, b = sym( @ ), c = sym( @ ,' d ' ), d = sym( '@ ' ) 在此，@ 分别代表具体数值 7/3 , pi/3 , pi*3^(1/3) ；而异同通过vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d))等来观察。 〖目的〗 ● 理解准确符号数值的创建法。 ● 高精度误差的观察。 〖解答〗 （1）x=7/3 x=7/3;a=x,b=sym(x),c=sym(x,'d'),d=sym('7/3'), a =

MATLAB作业4参考答案

MATLAB 作业四参考答案 1、 用2sin(103)y t =+在（0，3）区间内生成一组较稀疏的数据，并用一维数据插值的方法对给出的数据进行曲线拟合，并将结果与理论曲线相比较。 【求解】类似于上面的例子，可以用几乎一致的语句得出样本数据和插值效果。 >> t=0:0.2:3; y=sin(10*t.^2+3); plot(t,y,'o') ezplot('sin(10*t^2+3)',[0,3]); hold on x1=0:0.001:3; y1=interp1(t,y,x1,'spline'); plot(x1,y1) 由于曲线本身变换太大，所以在目前选定的样本点下是不可能得出理想插值效果的，因为样 本数据提供的信息量不够。为了得到好的插值效果，必须增大样本数据的信息量，对本例来 说，必须在快变化区域减小样本点的步长。 >> hold off t=[0:0.1:1,1.1:0.04:3]; y=sin(10*t.^2+3); plot(t,y,'o') ezplot('sin(10*t^2+3)',[0,3]); hold on x1=0:0.001:3; y1=interp1(t,y,x1,'spline'); plot(x1,y1) 2、 用2422 3 1(,)sin()3x y f x y e xy x y x y --= ++原型函数生成一组网络数据或随机数据，分别拟合出曲面，并和原曲面进行比较。 【求解】由下面的语句可以直接生成一组网格数据，用下面语句还可以还绘制出给定样本点是三维表面图。 >> [x,y]=meshgrid(0.2:0.2:2); z=exp(-x.^2-y.^4).*sin(x.*y.^2+x.^2.*y)./(3*x.^3+y); surf(x,y,z) 选择新的密集网格，则可以通过二元插值得出插值曲面。对比插值结果和新网格下的函数值精确解，则可以绘制出绝对插值误差曲面。由插值结果可见精度是令人满意的。 >> [x1,y1]=meshgrid(0.2:0.02:2); z1=interp2(x,y,z,x1,y1,'spline'); surf(x1,y1,z1) >> z0=exp(-x1.^2-y1.^4).*sin(x1.*y1.^2+x1.^2.*y1)./(3*x1.^3+y1); surf(x1,y1,abs(z1-z0)) 现在假设已知的样本点不是网格形式分布的，而是随机分布的，则可以用下面语句生成样本点，得出分布的二维、三维示意图。 >> x=0.2+1.8*rand(400,1); y=0.2+1.8*rand(400,1); % 仍生成(0.2,2) 区间的均匀分布随机数 z=exp(-x.^2-y.^4).*sin(x.*y.^2+x.^2.*y)./(3*x.^3+y); plot(x,y,'x') figure, plot3(x,y,z,'x')

matlab与数学实验大作业

《数学实验与MATLAB》 ——综合实验报告 实验名称：不同温度下PDLC薄膜的通透性 与驱动电压的具体关系式的研究学院：计算机与通信工程学院 专业班级： 姓名： 学号： 同组同学： 2014年 6月10日

一、问题引入 聚合物分散液晶(PDLC)是将低分子液晶与预聚物Kuer UV65胶相混合，在一定条件下经聚合反应，形成微米级的液晶微滴均匀地分散在高分子网络中，再利用液晶分子的介电各向异性获得具有电光响应特性的材料，它主要工作在散射态和透明态之间并具有一定的灰度。聚合物分散液晶膜是将液晶和聚合物结合得到的一种综合性能优异的膜材料。该膜材料能够通过驱动电压来控制其通透性，可以用来制作PDLC型液晶显示器等，具有较大的应用范围。已知PDLC薄膜在相同光强度及驱动电压下，不用的温度对应于不同的通透性，不同温度下的阀值电压也不相同。为了尽量得到不同通透性的PDLC薄膜，有必要进行温度对PDLC薄膜的特性的影响的研究。现有不同温度下PDLC 薄膜透过率与驱动电压的一系列数据，试得出不同温度下PDLC薄膜通透性与驱动电压的具体关系式，使得可以迅速得出在不同温度下一定通透性对应的驱动电压。 二、问题分析 想要得到不同温度下PDLC薄膜通透性与驱动电压的具体关系式可以运用MATLAB多项式农合找出最佳函数式，而运用MATLAB多项式插值可以得出在不同温度下一定通透性所对应的驱动电压。 三、实验数据 选择10、20、30摄氏度三个不同温度，其他条件一致。

（1）、10摄氏度 实验程序： x=2:2:40; y=[5.2,5.4,5.8,6.4,7.2,8.2,9.4,10.8,12.2,14.0,16.6,22.0, 30.4,39.8,51.3,55.0,57.5,58.8,59.6,60.2]; p3=polyfit(x,y,3); p5=polyfit(x,y,5); p7=polyfit(x,y,7); disp('三次拟合函数'),f3=poly2str(p3,'x') disp('五次拟合函数'),f5=poly2str(p5,'x') disp('七次拟合函数'),f7=poly2str(p7,'x') x1=0:1:40; y3=polyval(p3,x1); y5=polyval(p5,x1); y7=polyval(p7,x1); plot(x,y,'rp',x1,y3,'--',x1,y5,'k-.',x1,y7); legend('拟合点','三次拟合','五次拟合','七次拟合') 实验结果：

(完整版)matlab习题及答案(2)

第1章MATLAB概论 1.1与其他计算机语言相比较，MA TLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MA TLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MA TLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MA TLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MA TLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MA TLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。

matlab课后习题答案

第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff 或数值梯度gradient指令计算)(t y'曲线绘制 y'，然后把)(t y和)(t 在同一张图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1); yc=gradient(y)/dt;

plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0 sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间内的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间内的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =

MATLAB作业答案

题4.1 控制系统结构如图4.A 所示： 图4.A 习题 （1） 利用MA TLAB 对以上单位负反馈控制系统建立传递函数模型； （2） 将第一问中求得的传递函数模型转化为零极点增益形式和状态空间形式。 解：（1）MA TLAB 程序代码如下： z=[-1];p=[-1,-3];k=2 [num,den]=zp2tf(z,p,k) g_tf=tf(num,den) 程序运行结果如下： Transfer function: 2 s + 2 ------------- s^2 + 4 s + 3 由计算结果可知，系统的传递函数模型为： G （s ）= 3 4222 +++s s s (2) MA TLAB 程序代码如下： num=[2,2];den=[1,4,3] [z,p,k]=tf2zp(num,den) [A,B,C,D]=tf2ss(num,den) g_zpk=zpk(z,p,k) g_ss=ss(A,B,C,D) 程序运行结果如下： Zero/pole/gain: 2 (s+1) ----------- (s+3) (s+1) a = x1 x2 x1 -4 -3 x2 1 0 b = u1 x1 1 x2 0

c = x1 x2 y1 2 2 d = u1 y1 0 由计算结果可知，系统的零极点增益模和状态空间模型分别为： G （s ）=)1)(3()1(2+++s s s 和[]? ????=??????+??????-=x y u x x 2 2010 13- 4. 题5.1 设单位负反馈控制系统的开环传递函数为G （s ）= ) 177(2 ++s s s K 。 （1） 试绘制K=10、100时闭环系统的阶跃响应曲线，并计算稳态误差、上升时间、超调 量和过渡过程时间； （2） 绘制K=1000时闭环系统的阶跃响应曲线，与K=10、100所得结果相比较，分析增 益系数与系统稳定性的关系； （3） 利用roots 命令，确定使系统稳定时K 的取值范围。 解：（1） 题 6.1 已知单位负反馈控制系统的前向传递函数分别为G （s ）= ) 4)(2()1(2 +++s s s s K 、 ) 164)(1() 1(2 ++-+s s s s s K 和 ) 15)(7)(5)(3() 8(2 +++++s s s s s s K ，试利用MA TLAB 分别绘制各系 统的根轨迹图。 解：（1）G （s ）= ) 4)(2()1(2 +++s s s s K MA TLAB 程序代码如下： num=[1,1];den=conv([1,0],conv([1,0],[1,6,8])); sys=tf(num,den)

MATLAB期末大作业模版

《MATLAB》期末大作业 学院土木工程与建筑学院 专业 班级 姓名 指导教师李琳 2018 年 5 月16 日

明 作业内容题目2:问题描述：在[0 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先划定x的范围与间距，再列出y的表达式，利用plot函数绘制二维曲线。 （2）软件说明及源代码 >> x = 0:pi/10:2*pi; >>y = cos(5*x).*sin(x); >>plot(x,y) （3）实验结果 题目4:问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c (2)求f=0的解 （1）问题分析 这是一个符号函数显示以及符号函数的求解问题，第一问先定义常量与变量，在写出f表达式，利用pretty函数显示f。第二问利用solve函数求解f=0时的解。 （2）软件说明及源代码

第一问 >> syms a b c x; >> f=a*x^2+b*x+c; >> pretty(f) 第二问 >>syms a b c x; >>f=a*x^2+b*x+c; >> solve(f) （3）实验结果 1、 2、 题目5:问题描述：求积分 （1）问题分析 这是一个利用符号函数求积分的问题，先定义变量x，再列出I1表达式，利用int函数求在范围0到Pi/2上的积分。 （2）软件说明及源代码 >> syms x; >> I1=(1-2*sin(2*x))^0.5; >> int(I1,0,0.5*pi) （3）实验结果 题目6:问题描述：分别随机产生一个6×6的整数矩阵（元素可在［－20,20］之间），求该随机阵的秩，特征值和特征向量。 （1）问题分析 这是一个矩阵运算问题，先利用rand函数产生一个6*6的元素在-20到20

汽车理论课后作业答案MATLAB

汽车理论作业MA TLAB过程 汽车驱动力与阻力平衡图 加速度倒数-速度曲线图 0 10 20 30 40 50 60 70 80 90 100 u

汽车功率平衡图 u/(km/h)最高档等速百公里油耗曲线 Ua/(km/h)

燃油积极性-加速时间曲线 源程序： 《第一章》 m=3880; g=9.8; r=0.367; x=0.85; f=0.013; io=5.83; CdA=2.77; lf=0.218; Iw1=1.798; Iw2=3.598; Iw=lw1+lw2; ig=[6.09 3.09 1.71 1.00]; %变速器传动比 L=3.2; a=1.947; hg=0.9; n=600:1:4000;

T=-19.313+295.27* n/1000-165.44*(门/1000)人2+40.874*(门/1000)人3-3.8445*( n/IOOO).%; Ft1=T*ig(1)*io*x/r; %计算各档对应转速下的驱动力 Ft2=T*ig(2)*io*x/r; Ft3=T*ig(3)*io*x/r; Ft4=T*ig(4)*io*x/r; u1=0.377*r*n/(io*ig(1)); u2=0.377*r*n/(io*ig(2)); u3=0.377*r*n/(io*ig(3)); u4=0.377*r*n/(io*ig(4)); u=0:130/3400:130; F仁m*g*f+CdA*u”2/21.15;%计算各档对应转速下的驱动阻力 F2=m*g*f+CdA*u2.A2/21.15; F3=m*g*f+CdA*u3.A2/21.15; F4=m*g*f+CdA*u4.A2/21.15; figure(1); plot(u1,Ft1, '-r' ,u2,Ft2, '-m' ,u3,Ft3, '-k' ,u4,Ft4, '-b' ,u1,F1, '-r' ,u2,F2, '-m' ,u3,F3, ' k' ,u4,F4, '-b' , 'LineWidth' ,2) title( ' 汽车驱动力与阻力平衡图' ); xlabel( 'u_{a}/km.hA{-1}' ) ylabel( 'F/N' ) gtext( 'F_{t1}' ) gtext( 'F_{t2}' ) gtext( 'F_{t3}' ) gtext( 'F_{t4}' ) gtext( 'F_{f}+F_{w}' ) %由汽车驱动力与阻力平衡图知，他们无交点，u4在最大转速时达到最大 umax=u4(3401) Ft1max=max(Ft1); imax=(Ft1max-m*g*f)/(m*g) disp( ' 假设是后轮驱动' ); C=imax/(a/L+hg*imax/L) % 附着率 delta1=1+(Iw1+Iw2)/(m*rA2)+If*ig(1)*rA2*ioA2*x/(m*rA2); delta2=1+(Iw1+Iw2)/(m*rA2)+If*ig(2)*rA2*ioA2*x/(m*rA2); delta3=1+(Iw1+Iw2)/(m*rA2)+If*ig(3)*rA2*ioA2*x/(m*rA2); delta4=1+(Iw1+Iw2)/(m*rA2)+If*ig(4)*rA2*ioA2*x/(m*rA2); a1=(Ft1-F1)/(delta1*m); %加速度 a2=(Ft2-F2)/(delta2*m); a3=(Ft3-F3)/(delta3*m); a4=(Ft4-F4)/(delta4*m); h1=1./a1; %加速度倒数 h2=1./a2; h3=1./a3; h4=1./a4; figure(2);

MATLAB第一章作业答案

第一章M A T L A B概况与基本操作 1.选择题： (1)最初的MATLAB核心程序是采用A语言编写的。 A.FORTRAN B.C C.BASIC D.PASCAL (2)即将于2011年9月发布的MATLAB新版本的编号为D。 A.MATLAB 2011Ra B.MATLAB 2011Rb C.MATLAB R2011a D.MATLAB R2011b (3)在默认设置中，MATLAB中的注释语句显示的颜色是D。 A.黑色 B.蓝色 C.红色 D.绿色 (4)如果要以科学计数法显示15位有效数字，使用的命令是B。 A.format long B.format long e C.format long g D.format long d (5)在命令窗口新建变量a、b，如果只查看变量a的详细信息，使用的命令为B。 A.who a B.whos a C.who D.whos (6)如果要清除工作空间的所有变量，使用的命令为 C 。 A.clear B.clear all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上D。 A.冒号 B.逗号 C.空格 D.分号 (8)如果要重新执行以前输入的命令，可以使用B键。 A.下箭头↓ B.上箭头↑ C.左箭头← D.右箭头→ (9)如果要查询函数inv的相关信息，并显示在命令窗口，应使用命令A。 A.help B.lookfor C.doc D.三者均可 (10)如果要启动Notebook文档，下列D操作是可行的。 A.在命令窗口输入notebook命令 B.在命令窗口输入notebook filename命令 C.在Word中启动M-book文档 D.三者均可 2.填空题： (1)MATLAB是MATrix和LABoratory两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler和Jack Little。 (2)在MATLAB的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为紫色，注释显示的字体为绿色，错误信息显示的字体为红色。 (3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format compact，各行之间不添加空行的命令为format compact。 备注：本题布置给大家时有一点小错误，现在予以更正。 (4)在MATLAB中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组中每一行的各个元素，逗号的作用是分隔数组中每一行的各个元素或不同的命令，分号的作用是分隔数组中的各行或控制命令执行结果是否在命令窗口显示，冒号的作用是生成一维数组或显示全部元素，百分号的作用是注释行的开头，…的作用是把相邻两行的语句连接为一个命令，感叹号的作用是执行操作系统命令。 3.先建立自己的工作目录，再将自己的工作目录设置到MATLAB搜索路径下。请写出操作步骤或用Matlab命令实现。用help命令能查询到自己的工作目录吗？ 解：操作步骤： (1)在Windows环境中建立一个工作目录，如：c:\mywork； (2)启动MATLAB，在File菜单中选择Set Path…命令，显示出如下图的对话框：

matlab综合大作业(附详细答案)

《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.6724 1.2507 0.9247 -0.1766 1.1186 2.4286 1.5754 1.6546 5.3664 0.8087 4.2471 -1.2929 1.3493 0.7272 -0.6647 -0.3836 2)将矩阵A按列拉长得到矩阵B； 实验程序：B=A(:) 实验结果： B = 0.1349 -2.3312 1.2507 1.5754 -1.2929 3.3818 3.3783 0.9247 1.6546 1.3493 0.6266 2.4516 -0.1766 5.3664 0.7272 1.2279 3.1335 1.1186 0.8087 -0.6647 1.5888 -1.6724 2.4286 4.2471

-0.3836 3)提取矩阵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.1335 0.9247 1.1186 4)寻找矩阵A中大于0的元素；] 实验程序：G=A(find(A>0)) 实验结果： G = 0.1349 1.2507 1.5754 3.3818 3.3783 0.9247 1.6546 1.3493 0.6266 2.4516 5.3664 0.7272 1.2279 3.1335 1.1186 0.8087 1.5888 2.4286 4.2471 5)求矩阵A的转置矩阵D； 实验程序：D=A' 实验结果： D = 0.1349 -2.3312 1.2507 1.5754 -1.2929 3.3818 3.3783 0.9247 1.6546 1.3493 0.6266 2.4516 -0.1766 5.3664 0.7272 1.2279 3.1335 1.1186 0.8087 -0.6647 1.5888 -1.6724 2.4286 4.2471 -0.3836 6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E； 实验程序：E=flipud(fliplr(A)) 实验结果：