matlab作业

合集下载

matlab1-8章课后作业

MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件，是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境，已被广泛应用于不同领域。

Matlab的基本功能包括：数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。

Matlab第二章1．创建double的变量，并进行计算。

(1)a=87，b=190，计算a+b、a-b、a*b。

(2)创建uint8 类型的变量，数值与(1)中相同，进行相同的计算。

>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u)ans =15.4189>> sqrt(u-3*v)/(u*v)ans =0 + 0.4410i4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。

MATLAB作业

MATLAB作业⼀、必答题：1. MATLAB系统由那些部分组成？答：MATLAB系统主要由开发环境、MATLAB语⾔、MATLAB数学函数库、图形功能和应⽤程序接⼝五个部分组成。

2. 如何启动M⽂件编辑/调试器？答：在操作界⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M⽂件编辑/调试器将被启动。

在命令窗⼝中键⼊“edit”命令也可以启动M⽂件编辑/调试器。

3. 存储在⼯作空间中的数组能编辑吗？如何操作？答：存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。

4. 在MATLAB中有⼏种获得帮助的途径？答：在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器；（2）help命令：在命令窗⼝键⼊“help” 命令可以列出帮助主题，键⼊“help 函数名”可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗⼝键⼊“lookfor 关键词”可以搜索出⼀系列与给定关键词相关的命令和函数（4）模糊查询：输⼊命令的前⼏个字母，然后按Tab键，就可以列出所有以这⼏个字母开始的命令和函数。

5. 有⼏种建⽴矩阵的⽅法？各有什么优点？答：（1）以直接列出元素的形式输⼊；（2）通过语句和函数产⽣；（3）.在m⽂件中创建矩阵；（4）从外部的数据⽂件中装⼊。

6. 命令⽂件与函数⽂件的主要区别是什么？答：命令⽂件: M⽂件中最简单的⼀种，不需输出输⼊参数，⽤M ⽂件可以控制⼯作空间的所有数据。

运⾏过程中产⽣的变量都是全局变量。

运⾏⼀个命令⽂件等价于从命令窗⼝中顺序运⾏⽂件⾥的命令，程序不需要预先定义，只要依次将命令编辑在命令⽂件中即可。

函数⽂件：如果M⽂件的第⼀个可执⾏⾏以function开始，便是函数⽂件，每⼀个函数⽂件定义⼀个函数。

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“MATLAB”练习题要求:抄题、写出操作命令、运行结果，并根据要求，贴上运行图。

1、求230x e x -=的所有根.（先画图后求解）（要求贴图）>> solve('exp(x)—3*x^2',0)ans =—2*lambertw （—1/6＊3^（1/2））-2＊lambertw(—1,—1/6＊3^(1/2）)—2*lambertw （1/6*3^(1/2))3、求解下列各题：1）30sin lim x x x x ->->> sym x ；〉> limit((x-sin （x))/x^3)ans =1/62） (10)cos ,x y e x y =求>> sym x;>> diff （exp(x ）＊cos(x)，10）ans =（-32）＊exp(x)*sin （x)3）21/20(17x e dx ⎰精确到位有效数字）〉〉 sym x;〉〉 vpa((int(exp(x^2），x,0,1/2)），17)ans =0.544987104183622224）42254x dx x+⎰〉> sym x ；>〉 int （x^4/(25+x^2）,x)ans =125*atan （x/5) - 25＊x + x^3/35)求由参数方程arctan x y t⎧⎪=⎨=⎪⎩dy dx 与二阶导数22d y dx 。

〉> sym t;>> x=log(sqrt （1+t^2)）;y=atan(t)；〉> diff （y ，t ）/diff （x ，t)ans =1/t6)设函数y =f (x ）由方程xy +e y = e 所确定,求y ′(x ).>> syms x y ；f=x ＊y+exp(y ）—exp （1);〉> -diff(f,x ）/diff （f,y)ans =-y/（x + exp （y))7)0sin 2x e xdx +∞-⎰>〉 syms x ；>〉 y=exp(-x)＊sin(2*x ）；〉> int(y ，0,inf ）ans =2/58） 08x =展开（最高次幂为）〉> syms xf=sqrt(1+x);taylor(f,0，9)ans =— (429＊x^8）/32768 + (33*x^7）/2048 — (21＊x^6)/1024 +（7*x^5）/256 - (5＊x^4)/128 + x^3/16 - x^2/8 + x/2 + 19) 1sin (3)(2)x y e y =求〉> syms x y ；>〉 y=exp(sin （1/x)）;>〉 dy=subs （diff(y,3）,x ，2）dy =—0.582610)求变上限函数2x x ⎰对变量x 的导数.>> syms a t ；>〉 diff （int(sqrt(a+t),t,x ，x^2）)Warning: Explicit integral could not be found 。

MATLB实验作业

实验一MATLAB运算基础1、用逻辑表达式求下列分段函数的值。

t=0:0.5:2.5y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*(( t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2 )&(t<3))2、求[100，999]之间能被21整除的数的个数。

p=rem([100:999],21)==0;sum(p)3、建立一个字符串向量，删除其中的大写字母。

ch='KdDdfKaWdsfCI',k=find(ch>=' A'&ch<='Z'),ch(k)=[]4、输入矩阵，并找出A中大于或等于5的元素。

A=[1 2 3;4 5 6;7 8 9],[m,n]=find(A>=5),for j=1:length(m)x(j)=A(m(j),n(j))xend5、求矩阵的行列式值、逆和特征根。

a11=input('a11='),a12=input('a12 ='),a21=input('a21='),a22=input('a22 ='),A=[a11,a12;a21,a22],DA=det(A),IA=inv(A),EA=eig(A) 6、不采用循环的形式求出和式的数值解。

sum(2.^[0:63])实验二1、1行100列的Fibonacc数组a，a(1)=a(2)=1，a(i)=a(i-1)+a(i-2)，用for循环指令来寻求该数组中第一个大于10000的元素，并指出其位置i。

n=100;a=ones(1,n);for i=3:na(i)=a(i-1)+a(i-2);if a(i)>10000a(i),break;end;end,i2、编写M脚本文件，定义下列分段函数，并分别求出当（）、（）和（）时的函数值。

MATLAB与控制系统仿真大作业

>> G1=tf(num,den);
>> G1=zpk(G1)
G1 =
5 (s+0.6) ----------------(s+3) (s+2) (s+1)
>> num=[2 1]; den=[1 2.9 1]; G2=tf(num,den); >> G2=zpk(G2)
G2 =
2 (s+0.5)
>> xlabel('x') ylabel('exp(-2.*x)')
三、simulink 建模（每题 10 分，共 20 分）
1、已知单位负反馈的开环传递函数为 G(s) =
2 s2 + 4s
，试利用
simulink
建立系在单位阶跃输入作用下的模型。
要求答案包括：（1）simulink 建模结构图；
（2）在同一个坐标中的阶跃信号和响应曲线图。
2、已知系统的开环传递函数为
G1(s)
=
2s2 + 5s + 6 s2 + 2s + 3
、G2 (s)
=
s2
s+6 + 7s +1
，H
(s)
=
5(s + 2) s +10
求：建立 Simulink 仿真模型，并求出其系统在单位阶跃响应；
要求答案包括：（1）simulink 建模结构图；
x=
-2.9709
0.5491
3.6000
0.0509
2. 已知下列矩阵
2 3 1
−1 3 5
A
=

MATLAB平时作业(图文版)

MATLAB 平时作业第一章习题16. 以下两种说法对吗？（1）“MATLAB 的数值表达精度与其指令窗中的数据显示精度相同。

”答：此种说法错误。

MATLAB 提供了控制数据显示格式的控制指format ，该指令并不改变MATLAB 内存中变量的精度，只是改变其显示精度。

（2）“MATLAB 指令窗中显示的数值有效位数不超过7位。

”答：此种说法错误。

当变量小于1000时，使用format 或format short 后，或者默认情况下，变量的显示精度最多不超过7位，但显示精度不等于变量的精度。

7. 想要在MATLAB 中产生二维数组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321S ，下面哪些指令能实现目的？ S=[1,2,3;4,5,6;7,8;9]S=[1 2 3;4 5 6;7 8 9]S=[1，2，3；4，5，6；7，8，9] %整个指令在中文状态下输入答：操作如图：第1、2条指令可以实现，第3条指令不可实现。

第三章习题31．在MATLAB 中，先运行指令A=magic(3), B=[1,2,1;3,4,3;5,6,7], C=reshape(1:6,3,2)生成阵列33⨯A ，23⨯B ，23⨯C ，然后根据运行结果回答以下问题：运行结果如图：（1）计算A*B, B*A ，这两个乘积相同吗？计算结果如图：答：不同。

（2）计算A\B, B/A ，左除、右除结果相同吗？计算结果如图：答：不同。

（3）计算B( : ,[1,2]).*C和C.*B( : , [1,2])，这两个乘积相同吗？计算结果如图答：相同。

（4）计算A\A和A.\A，这两个计算结果相同吗？计算结果如图：答：相同。

（5）计算A\eye(3)和inv(A)，这两个计算结果相同吗？计算结果如图：答：不同。

（提示：根据对计算结果的目测回答问题）2．在MATLAB中，先运行A=[1, 2; 3, 4],b=0.5,C=[4, 2; 1, 0.5], 然后根据计算结果回答以下问题：创建数据步骤略（1）计算A^b和A.^b, 这两个计算结果相同吗？答：不同。

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习题及答案

现代计算方法Matlab 作业答案1.绘出函数f(x)=sin x x ，在[0,4]上的图形解：在M 文件输入：x=0:pi/100:4;y=x.*sin(x);plot(y)运行2. 求3x +2x +5 = 0的根解：在命令窗口输入：>> solve('x^3+2*x+5=0')ans =((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3) - 2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3))1/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) - ((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)/2 -(3^(1/2)*i*(2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) + ((108^(1/2)*707^(1/2))/108 -5/2)^(1/3)))/21/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) - ((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)/2 +(3^(1/2)*i*(2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) + ((108^(1/2)*707^(1/2))/108 -5/2)^(1/3)))/23.321436min x x x z ++=120..321=++x x x t s301≥x5002≤≤x203≥x解：运用单纯形法计算此题，首先把约束条件化成标准形式：,,,,,205030120654321635241321≥=-=+=-=++x x x x x x x x x x x x x x x（1）在M 文件输入SimpleMthd 函数：function [x,minf] = SimpleMthd(A,c,b,baseVector)sz = size(A);nVia = sz(2);n = sz(1);xx = 1:nVia;nobase = zeros(1,1);m = 1;for i=1:nViaif (isempty(find(baseVector == xx(i),1)))nobase(m) = i;m = m + 1;else;endendbCon = 1;M = 0;while bConnB = A(:,nobase);ncb = c(nobase);B = A(:,baseVector);cb = c(baseVector);xb = inv(B)*b;f = cb*xb;w = cb*inv(B);for i=1:length(nobase)sigma(i) = w*nB(:,i)-ncb(i);end[maxs,ind] = max(sigma);if maxs <= 0minf = cb*xb;vr = find(c~=0 ,1,'last');for l=1:vrele = find(baseVector == l,1);if (isempty(ele))x(l) = 0;elsex(l)=xb(ele);endendbCon = 0;elsey = inv(B)*A(:,nobase(ind));if y <= 0disp('不存在最优解!');x = NaN;minf = NaN;return;elseminb = inf;chagB = 0;for j=1:length(y)if y(j)>0bz = xb(j)/y(j);if bz<minbminb = bz;chagB = j;endendendtmp = baseVector(chagB);baseVector(chagB) = nobase(ind);nobase(ind) = tmp;endendM = M + 1;if (M == 1000000)disp('找不到最优解!');x = NaN;minf = NaN;return;endend（2）在命令窗口输入：clear allA=[1 1 1 0 0 0;1 0 0 -1 0 0;0 1 0 0 1 0;0 0 1 0 0 -1];c=[6 3 4 0 0 0];b=[120;30;50;20];[xm,mf]=SimpleMthd(A,c,b,[3 4 5 6])xm =0 50 70mf =4304.计算下面函数在区间(0，1)内的最小值。

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综合大作业(附详细答案)

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编程作业

《Matlab编程训练》作业专业学生姓名班级学号指导教师完成日期实训一 MATLAB 语言介绍和数值计算1．先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存变量。

122sin 851z e =+.2. 已知 1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦，求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） A+6*B:A-B+I:(2)A*B和A.*BA*B程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]c=A*B结果：A.*B程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]D=A.*B结果：(3)A^3和A.^3A^3程序：A=[12 34 -4;34 7 87;3 65 7]E=A^3结果：A．^3程序：A=[12 34 -4;34 7 87;3 65 7]C=A.^3(4)A/B及B\AA/B程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]C=A/B结果：B\A程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]D=B\A结果：(5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并（3）保存变量（mat文件）程序：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7];C=B*inv(A);D=C(2:3,2:3)结果：3. 求得矩阵⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=34157864653434533145A 的每行最大元素所在的位置？(至少两种方法) 第一种：A=[5 14 33;45 43 3;65 4 6;78 15 34][RowMax Order]=max(A')结果：第二种：A=[5 14 33;45 43 3;65 4 6;78 15 34][Max_num,index]=max(A,[],2)结果：实训二 MATLAB 编程基础1. 求[25,1258]之间能被15整除的数的个数。

MATLAB第一章作业答案

第一章M A T L A B概况与基本操作1.选择题：(1)最初的MATLAB核心程序是采用A语言编写的。

(2)即将于2011年9月发布的MATLAB新版本的编号为D。

2011Ra 2011Rb R2011a R2011b(3)在默认设置中，MATLAB中的注释语句显示的颜色是D。

A.黑色B.蓝色C.红色D.绿色(4)如果要以科学计数法显示15位有效数字，使用的命令是B。

long long e long g long d(5)在命令窗口新建变量a、b，如果只查看变量a的详细信息，使用的命令为B。

a a(6)如果要清除工作空间的所有变量，使用的命令为 C 。

all C.两者都可 D.两者都不可(7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上D。

A.冒号B.逗号C.空格D.分号(8)如果要重新执行以前输入的命令，可以使用B键。

A.下箭头↓B.上箭头↑C.左箭头←D.右箭头→(9)如果要查询函数inv的相关信息，并显示在命令窗口，应使用命令A。

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。

备注：本题布置给大家时有一点小错误，现在予以更正。

（完整版）武汉大学matlab期末课程作业

（完整版）武汉⼤学matlab期末课程作业“MATLAB及其应⽤”课程作业院系：姓名：学号：联系⽅式：1. 请指出如下5个变量名中，哪些是合法的？ abcd-2 xyz_3 3chan a 变量 ABCDefgh答：xyz_3，ABCDefgh 是合法的。

2. 在命令窗中，运⾏命令a=sqrt(2)。

然后请回答以下问题：计算结果a吗？该计算结果只是5近似吗？请在命令窗中，显⽰出具有最多位有MATLAB 数值结果显⽰的默认设置。

该计算结果只是5近似。

3. 命令clear, clf, clc 各有什么⽤处？答：clear 可以清除matlab ⼯作空间中保持的变量。

clf 可以清除图形窗。

clc 清除命令窗中显⽰内容。

4. 想要在MATLAB 中产⽣⼆维数组=987654321S ，下⾯哪些命令能实现⽬的？ S=[1,2,3;4,5,6;7,8;9]S=[1 2 3;4 5 6;7 8 9]S=[1，2，3；4，5，6；7，8，9] %整个命令在中⽂状态下输⼊答：第⼆条S=[1 2 3;4 5 6;7 8 9]能实现⽬的。

中⽂状态下逗号不是有效字符。

1．说出以下四条命令产⽣的结果各属于哪种数据类型，是“双精度”对象，还是“符号”对象？3/7+0.1, sym(3/7+0.1), vpa(sym(3/7+0.1),4), vpa(sym(3/7+0.1))答：3/7+0.1结果是双精度。

sym(3/7+0.1)结果是符号。

vpa(sym(3/7+0.1),4)结果是符号。

vpa(sym(3/7+0.1))结果是符号。

过程如图：2．已知a1=sin(sym(pi/4)+exp(sym(0.7)+sym(pi/3)))产⽣精准符号数字，请回答：以下产⽣的各种符号数哪些是精准的？若不精准，误差⼜是多少？能说出产⽣误差的原因吗？a2=sin(sym(pi/4)+exp(sym(0.7))*exp(sym(pi/3)))a3=sin(sym('pi/4')+exp(sym('0.7'))*exp(sym('pi/3')))a4=sin(sym('pi/4')+exp(sym('0.7+pi/3')))a5=sin(sym(pi/4)+exp(sym(0.7+pi/3)))a6=sin(sym(pi/4)+sym(exp(0.7+pi/3)))a7=sin(sym(pi/4+exp(0.7+pi/3)))a8=sym(sin(pi/4+exp(0.7+pi/3)))（提⽰：可⽤vpa观察误差；注意数位的设置）。

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

实验一 MATLAB 工作环境熟悉及简单命令的执行一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。

二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条（Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。

三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。

2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。

（1）（365-522-70）3（2）>>area=pi*^2（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。

m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>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::6*pi;>>y=5*sin(x);>>plot(x,y)（6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。

四、思考题1、以下变量名是否合法为什么（1）x2（2）3col（3）_row（4）for2、求以下变量的值，并在MATLAB 中验证。

离散卷积MATLAB编程作业

一、离散卷积MATLAB编程作业（1）请编程实现一维卷积功能，即自己编写MATLAB中的CONV函数，自己写出一个CONV_NEW函数。

要求，请自己学习MATLAB向量运算操作，基于学习到的内容，用尽量少的代码实现CONV_NEW。

(a)请用例子检查你写的CONV_NEW是否与CONV的输出是否一致。

(b)请用两个10000维数组检测对比你写的CONV_NEW和CONV的时间复杂度。

请查资料大致回答，为什么MATLAB自带的CONV更快一些？（2）请查询自学二维卷积的知识(a)请写一个函数CONV2_NEW实现MATLAB中的二维卷积函数CONV2功能，请用例子检查你写的CONV2_NEW与CONV2的输出是否一致(b)请用此函数，做以下两个卷积核与附件图片lenna.jpg的卷积，将卷积结果打印出来。

并想一下为什么是这样。

（请查阅imread和imwrite函数）两点提示：由于imread和imwrite函数读入和写出的数值范围都是0到255，因此需要将结果归一化。

请参考一下函数程序读入写出例程。

A = imread(‘lenna.jpg’);B = double(A);Imwrite(uint8(B),’test.bmp’);（3）请编写MATLAB程序H = DECONV(X, Y)要求输入序列X和Y，输出序列H，使Y=X*H，若没有H满足要求，则输出错误信息”No sequence satisfies this condition.”二、傅里叶变换编程题对于卷积公式：x(t)*h(t) = y(t)，则根据傅里叶变换和卷积性质，在已知x(t)和y(t)的情况下，h(t)可以写成：ℎ(t)=12π∫Y(jω)X(jω)e jωt dω+∞−∞（1）若x(t)和y(t)分别是如下函数，求h(t)。

（2）请用MATLAB画出h(t)在t=[-50,50]的函数图像。

（提示：用分成小段求矩形面积的形式近似求积分，注意规避分母等于0时候的情况）。

MATLAB大作业

M A T L A B大作业(总15页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--MATLAB大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。

（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。

打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。

（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。

不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。

（4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。

（5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。

第一类：绘制图形。

（B级）问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。

斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。

问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。

其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。

接下来对其余三个小三角形重复上述操作，如图所示。

问题三：其他分形曲线或图形。

分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。

matlab作业练习及答案(有程序和截图)

1. 给出如下方程组：16323313915111081697612534141543s u p w s u p w s u p w s u p w +++=+++=+++=+++=求s ，u ，p ，w 的值，并求出系数行列式。

[答案：s ＝-0.1258，u ＝-8.7133，p ＝11.2875，w ＝-0.0500。

行列式=7680。

]>> A=[16 32 33 13;5 11 10 8;9 7 6 12;34 14 15 1];>> B=[91 16 5 43]';>> C=A\B2. 求矩阵H1()H X X X X -''=其中1731565419289121110X ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦〔答案：H ＝〔0.7294，0.9041，0.4477，0.9188〕'〕>> X=[17 31 5;6 5 4;19 28 9;12 11 10]; >> H=X/(X'*X)*X'>> X=[17 31 5;6 5 4;19 28 9;12 11 10]; >> H=X*inv((X'*X))*X'3．用MA TLAB语言实现下面的分段函数,()/,||,h x Dy f x h Dx x Dh x D>⎧⎪==≤⎨⎪-<-⎩。

y=h*(x>D)+h/D*x.x>=-D&x<=D)-h*(x<-D)4．636263 0S=2124822 ii==++++++∑①试不采用循环的形式，用数值方法求出上式的解。

②由于数值方法采用double形式进行计算，难以保证有效数字，试采用符号运算的方法求该式的精确解。

并给出保留16位有效数字的结果。

>> sum(2.^[1:63])>> sum(sym(2).^[1:63])5．编写一个矩阵相加函数mat_add（），使其具体的调用格式为，要求该函数能接受任意多个矩阵进行加法运算。

4个Matlab实验的作业

现给出bpsk、qpsk及“书上习题”的调制解调程序，理解各程序，完成以下习题。

将程序运行结果及各题目的解答写入word中：1.用matlab运行书上习题中的“bpskqpsk125.m”（a）说明bpsk、qpsk解调判决方法（b）误比特率为1e-2、1e-3及1e-4时的Eb/N0分别是多少？（c）从物理意义上说明为什么bpsk、qpsk误比特率曲线是重叠的（d）当samples减少为100000，10000，1000时观察误比特率曲线的变化，你得出什么结论。

2.用matlab运行“bpsk.m”、“qpsk.m”（a）在各程序中标注“注释”处加上注释（英文或中文）（b）说明加性高斯白噪声的产生方法，请再给出一种加性高斯白噪声的产生方法，并验证其正确性。

（c）参考“bpskqpsk125.m”的画图功能，给出Eb/N0---误比特率曲线和高斯信道下的理论误比特率曲线。

（d）观察nd及nloop参数变化时，曲线的现象，并说明原因。

（e）画出不同信噪比条件下的的星座图，解释其对误码率的影响。

（f）通过程序画出QPSK和BPSK的Eb/N0---误比特率曲线，观察曲线的现象，能得出什么结论。

3.若信源是你的学号，结合程序说明其在qpsk (调制mod)和（解调demod）子程序中的具体实现过程。

4.针对题目2中的BPSK、QPSK，(a)若信道使信号幅度呈瑞利衰落，画出Eb/N0---误比特率曲线和瑞利衰落下的理论误比特率曲线，说明与题目2观察结果的异同，并说明原因。

(b)若信道使信号幅度呈莱斯衰落，更改K值的大小，画出Eb/N0---误比特率曲线、瑞利衰落下的理论误比特率曲线和高斯信道下的理论误比特率曲线，观察曲线的现象，能得出什么结论。

现给出循环码及卷积码的编解码程序，理解各程序，完成以下习题。

将程序运行结果及各题目的解答写入word中：1.用matlab运行书上习题中的“clockcode.m”（a）说明（7，4）码的纠错检错方法（b）在程序中标注“注释”处加上注释（英文或中文）（c）对于编码和未编码的情况，误比特率为1e-2、1e-3及1e-4时的Eb/N0分别是多少？（d）从物理意义上说明编码增益问题（e）当采用（15，11）码时，观察与（7，4）码相比编码增益的变化。