Matlab语言及应用大作业02答案

合集下载

matlab习题二带答案

matlab习题二带答案Matlab习题二带答案Matlab是一种强大的数值计算和科学编程软件，广泛应用于各个领域的科学研究和工程实践中。

掌握Matlab的使用对于科学家和工程师来说至关重要。

为了帮助大家更好地掌握Matlab，下面将介绍一些常见的Matlab习题，并提供相应的答案。

1. 习题一：计算平均值编写一个Matlab函数，输入一个向量，输出该向量的平均值。

答案：```matlabfunction avg = computeAverage(vector)avg = sum(vector) / length(vector);end```2. 习题二：矩阵操作编写一个Matlab函数，输入两个矩阵A和B，输出它们的乘积C。

答案：```matlabfunction C = matrixMultiplication(A, B)[m, n] = size(A);[~, p] = size(B);C = zeros(m, p);for i = 1:mfor j = 1:pfor k = 1:nC(i, j) = C(i, j) + A(i, k) * B(k, j);endendendend```3. 习题三：图像处理编写一个Matlab函数，输入一张彩色图像，输出该图像的灰度图像。

答案：```matlabfunction grayImage = convertToGray(image)grayImage = rgb2gray(image);end```4. 习题四：数据拟合给定一组数据点(x, y)，编写一个Matlab函数，拟合这些数据点为一条直线，并返回拟合直线的斜率和截距。

答案：```matlabfunction [slope, intercept] = fitLine(x, y)n = length(x);sx = sum(x);sy = sum(y);sxy = sum(x .* y);sxx = sum(x .* x);slope = (n * sxy - sx * sy) / (n * sxx - sx^2);intercept = (sy - slope * sx) / n;end```5. 习题五：数值积分编写一个Matlab函数，输入一个函数f(x)和积分区间[a, b]，输出该函数在该区间上的数值积分结果。

MATLAB课后习题集附标准答案

第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息.（3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数.注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？答：在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法？各有什么优点？答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷；（2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改；（3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵；（4）通过数据文件建立，该方法可以调用由其他软件产生数据.3.3在进行算术运算时，数组运算和矩阵运算各有什么要求？答：进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。

内蒙古国开电大一网一《Matlab语言及其应用》形考作业二

试卷总分:100 得分:97.33
一、判断题（6小题，每小题5分，共30分）
1.标量可以和数组进行乘除运算。

答案:正确
2.假设二维数组A，则A( {图} 代表数组A的全部元素。

答案:正确
3.Matlab中对虚数虚单元的描述使用 i或j 。

答案:正确
4.符号表达式的四则运算和数值运算一样，用+、－、*、/、∧运算符实现。

答案:正确
5.使用函数zeros(5)生成的是一个具有5个元素的向量。

答案:正确
6.在sin(x)运算中，x是角度。

答案:错误
二、单选题（6小题，每小题5分，共30分）
7.ones(n,m)函数是用来产生特殊矩阵的，由它形成的矩阵称为（）。

A.单位矩阵
B.全零矩阵
C.对角矩阵
D.空矩阵
答案:A
8.角度x=[20 30 60]，计算其余弦函数的运算为 {图} 。

A.COS（deg2rad(x)）
B.COS(x)
C.cos(x)
D.cos(deg2rad(x))
答案:D
9.如果x=1:-2:-8,则x(1)和x(5)分别是 {图} 。

A.1，-8
B.-1,-7
C.-1,-8
D.1,-7
答案:D
10.如果想将矩阵A第1行第1列的元素改为50，可以通过（）语句来完成。

matelab作业2参考答案

matelab作业2参考答案Matlab作业2参考答案Matlab作业2是一项综合性的任务，要求学生运用Matlab编程语言解决一系列数学问题。

本文将为大家提供一份参考答案，帮助学生更好地理解和完成这项作业。

首先，我们将讨论作业的第一个问题，即给定一个矩阵A，求解其特征值和特征向量。

在Matlab中，可以使用eig函数来实现这一功能。

例如，假设我们有一个3×3的矩阵A，可以按照以下方式计算其特征值和特征向量：```A = [1 2 3; 4 5 6; 7 8 9];[eigenvectors, eigenvalues] = eig(A);```在上述代码中，变量eigenvectors将存储A的特征向量，而变量eigenvalues 将存储A的特征值。

通过打印这两个变量的值，我们可以得到矩阵A的特征值和特征向量。

接下来，我们将探讨作业的第二个问题，即求解线性方程组。

假设我们有一个3×3的系数矩阵A和一个3×1的常数向量b，我们需要求解方程组Ax=b。

在Matlab中，可以使用backslash运算符来求解线性方程组。

例如，假设我们有以下方程组：```A = [1 2 3; 4 5 6; 7 8 9];b = [10; 20; 30];x = A \ b;```在上述代码中，变量x将存储方程组的解。

通过打印变量x的值，我们可以得到方程组的解。

此外，作业的第三个问题要求学生使用Matlab绘制函数图像。

在Matlab中，可以使用plot函数来实现这一功能。

例如，假设我们要绘制函数y=sin(x)，其中x的取值范围为0到2π，可以按照以下方式绘制函数图像：```x = 0:0.1:2*pi;y = sin(x);plot(x, y);```在上述代码中，变量x将存储x的取值范围，变量y将存储对应的函数值。

通过调用plot函数，我们可以将函数y=sin(x)的图像绘制出来。

matlab编程与工程应用课后答案-第2章

matlab编程与工程应用课后答案-第2章1%计算空心球体积r0 = 3;ri = 2;volume = 4*pi/3*(r0^3-ri^3) —————————————————————————————————————2%计算过氧化氢分子量H = 15.9994;O = 1.0079;H2O2 = 2*H + 2*O —————————————————————————————————————3%计算字符创长度string = input('enter string: ','s');stringlength = length(string) —————————————————————————————————————%4 输入实数，fprintf输出这个变量，格式为两位小数num = input('enter a num: ');fprintf('The num is %.2f\n',num) —————————————————————————————————————6fprintf('不指定宽度: %f\n',12345.6789)不指定宽度: 12345.678900fprintf('10个字符宽度，4个小数位: %10.4f\n',12345.6789)10个字符宽度，4个小数位: 12345.6789fprintf('10个字符宽度，2个小数位: %10.2f\n',12345.6789)10个字符宽度，2个小数位: 12345.68fprintf('6个字符宽度，4个小数位: %6.4f\n',12345.6789)6个字符宽度，4个小数位: 12345.6789fprintf('2个字符宽度，4个小数位: %2.4f\n',12345.6789)2个字符宽度，4个小数位: 12345.6789 —————————————————————————————————————7fprintf('不指定宽度：%int16\n',12345)不指定宽度：12345nt16fprintf('不指定宽度：%d\n',12345)不指定宽度：12345fprintf('5个字符宽度：%5d\n',12345)5个字符宽度：12345fprintf('8个字符宽度：%8d\n',12345)8个字符宽度：12345fprintf('3个字符宽度：%3d\n',12345)3个字符宽度：12345 —————————————————————————————————————8x = 12.34;y = 4.56;fprintf('x is %.3f\n',x)x is 12.340fprintf('x is %2.0f\n',x)x is 12fprintf('y is %.1f\n',y)y is 4.6fprintf('y is %-6.1f!\n',y)y is 4.6 ! —————————————————————————————————————9%计算矩形面积x = input('enter length:');y = input('enter borad:');area = x*y;fprintf('area is %.2f\n',area) —————————————————————————————————————10string = input('what is your name?','s');fprintf('WoW,your name is %s\n',string) —————————————————————————————————————11string = input('enter your string:','s');fprintf('your string was:''%s''\n',string) %输出单引号是应采用''格式—————————————————————————————————————12v = input('enter the flow in m^3/s:');fprintf('a flow rate of %.3f meters per sec\nis equivalent to %.3f feet per sec\n',v,v/0.028) —————————————————————————————————————13incomeyear = input('enter your income per year: ');fprintf('the range of food expenditure per year %f~%f\nthe range of food expenditure per mounth %f~%f\n',incomeyear*0.08,incomeyear*0.1,incomeyear* 0.08/12,incomeyear*0.1/12)—————————————————————————————————————14wight = input('wight of plant: ');area = input('area of wing: ');fprintf('W/A: %f kg/m^2\n',wight/area)—————————————————————————————————————15x= 10;y =22;plot(x,y,'g+')—————————————————————————————————————16x = -2:0.1:2;plot(x,exp(x))xlabel('x')ylabel('y')title('y=e^x') —————————————————————————————————————17x = 1:5:100;y = sqrt(x);figure(1)plot(x,y)figure(2)bar(x,y)—————————————————————————————————————18略—————————————————————————————————————19x1 = linspace(0,pi,10);figure(1)plot(x1,sin(x1))x2 = linspace(0,pi,100);figure(2)plot(x2,sin(x2)) —————————————————————————————————————20mat = [1000 2000 3000 5000 10000;288 281 269 256 223]';x = mat(:,1);y = mat(:,2);plot(x,y)xlabel('high')ylabel('tempture')title('high-tempture') —————————————————————————————————————21mat1 = randi([50,100],3,6)save randfile.dat mat1 -ascii;mat2 = randi([50,100],2,6)save randfile.dat mat2 -ascii -append;load randfile.dat;;randfile —————————————————————————————————————22mat = rand(2,3)*4-1;save testtan.dat mat -ascii;load testtan.dat;mattan = tan(testtan)—————————————————————————————————————23mat = [89 42 49 55 72 63 68 77 82 76 67;90 45 50 56 59 62 68 75 77 75 66;91 44 43 60 60 60 65 69 74 70 70]save hightemp.dat mat -ascii;load hightemp.dathightemp(:,1) = hightemp(:,1)+1900;hightempsave y2ktemp.dat hightemp -ascii—————————————————————————————————————24%24 Calculates y as a function of xfunction y = fn(x)y = x^3-4*x^2+sin(x); —————————————————————————————————————25%25 Converts from MWh to GJfunction gj = mwh_to_gj(mwh)gj = 3.6*mwh;—————————————————————————————————————26%26 converta from inch/h to meter/sfunction meter_sec = converta(inch_hour)meter_sec = inch_hour*5280*0.3048/3600;—————————————————————————————————————27function Tn = fn27(P,i,n)Tn = P*(1+i)*n;—————————————————————————————————————28略—————————————————————————————————————29function V = fn29(Pt,Ps)V = 1.016*sqrt(Pt-Ps);—————————————————————————————————————30function THR = fn30(A)THR = (220-A)*0.6;—————————————————————————————————————31function outdate = fn31(n)outdate = sqrt(2*pi*n)*(n/exp(1))^n;—————————————————————————————————————32%32 脚本n = input('enter the number of units: ');Cn = costn(n);fprintf('the cost for %d units will be $%.2f\n',n,Cn)%32 mygcost函数function Cn = costn(n)Cn = 5*n^2-44*n+11; —————————————————————————————————————33%33 脚本rain = input('enter the rain: ');snow = fn33(rain);fprintf('the snow is %f\n',snow)%33 调用函数function snow = fn33(rain)snow = rain*6.5; —————————————————————————————————————34%34 脚本s = input('enter long_side: ');V = fn34(s);fprintf('volume is %.2f\n',V)%34 函数function V = fn34(s)V = sqrt(2)*s^3/12; —————————————————————————————————————35%35 pickone(x) returns a random element from vector xfunction outdate = pickone(x)outdate = randi([x(1) x(end)],1,1); —————————————————————————————————————36function outdate = vecout(x)outdate = x:1:x+5;—————————————————————————————————————37%37b = input('enter the first side: ');c = input('enter the second side: ');alpha = input('enter the angle between them: ');a = sqrt(b^2+c^2-2*b*c*cosd(alpha)); %切记这里用sindfprintf('the third side is %.3f\n',a)—————————————————————————————————————38略—————————————————————————————————————39%39mat = [90.5792 27.8498 97.0593;12.6987 54.6882 95.7167;91.3376 95.7507 48.5376;63.2359 96.4889 80.02809.7540 15.7613 14.1886];save floatnums.dat mat -ascii;load floatnums.dat;floatnums = round(floatnums')save intnums.dat floatnums -ascii—————————————————————————————————————40%40costssales = [1100 800;1233 650;1111 1001;1222 1300;999 1221];fprintf('there were %d quarters in the file\n',length(costssales))x = costssales(:,1);y = costssales(:,2);plot(x,'ko')hold onplot(y,'b*')mat = fliplr(costssales)';save newfile.dat mat -ascii。

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

matlab习题2答案MATLAB习题2答案在MATLAB习题2中，我们学习了如何使用MATLAB编程语言解决各种数学问题和数据分析任务。

在本篇文章中，我们将探讨一些常见的习题，并给出相应的答案和解释。

第一道习题是关于矩阵运算的。

我们需要编写一个MATLAB程序，计算两个矩阵的乘积。

在MATLAB中，可以使用矩阵乘法运算符“*”来实现这个功能。

例如，如果我们有两个矩阵A和B，我们可以使用以下代码来计算它们的乘积：```matlabC = A * B;```接下来，我们需要解决一个关于插值的问题。

我们需要编写一个程序，使用插值方法来估算给定点的函数值。

在MATLAB中，可以使用interp1函数来实现这个功能。

例如，如果我们有一个已知的函数y和一组离散的点x，我们可以使用以下代码来进行插值计算：```matlabx_interp = 0:0.1:10;y_interp = interp1(x, y, x_interp, 'spline');```最后，我们需要解决一个关于数据可视化的问题。

我们需要编写一个程序，绘制给定数据的散点图和拟合曲线。

在MATLAB中，可以使用scatter和plot函数来实现这个功能。

例如，如果我们有一组数据点x和y，我们可以使用以下代码来进行数据可视化：```matlabscatter(x, y);hold on;p = polyfit(x, y, 1);y_fit = polyval(p, x);plot(x, y_fit, 'r');```通过以上习题的练习，我们可以更好地掌握MATLAB编程语言的基本语法和常用函数。

希望本篇文章对大家有所帮助，也希望大家能够在日常工作和学习中充分利用MATLAB的强大功能，提高工作效率和解决问题的能力。

（完整版）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语言及应用大作业02答案

MATLAB语言及应用大作业姓名：学号：班级：1.利用plot函数在区间[0,2π]同时绘制x=sin(t)和y=cos(t),z=sin(t)+cos(t)的图形。

要求：1)对图像x轴和y轴分别标注“时间”和“函数值”2)对曲线加图例标注,图例位置自动定位（12分）。

>> t=[0:0.1*pi:2*pi];>> x=sin(t);>> y=cos(t);>> z=sin(t)+cos(t);>> plot(t,x,t,y,t,z)>> xlabel('时间')>> ylabel('函数值')>> legend('x=sin(t)','y=cos(t)','z=sin(t)+cos(t)')2. 求2!+4!+6!+8!+10! （12分）sum=0;for i=2:2:10pdr=1;for k=2:ipdr=pdr*k;endsum=sum+pdr;end在M文件编辑窗口编写上述程序，保存为exam02.m文件；在MATLAB命令窗口执行命令：>>exam02>>sumsum =36698663．编写m 函数，计算函数值（12分）221,10002,,0≥<≤<≤<⎪⎪⎩⎪⎪⎨⎧-=x x x x x x y 编写函数代码function f=exam03(t) if(t<0) f=0;elseif(t>=0)&(t<1) f=t;elseif(t>=1)&(t<=2) f=2-t; else f=0; end在M 文件编辑窗口编写上述程序，保存为exam03.m 文件在MATLAB 命令窗口输入： >> x=1;>> y=exam07(x) y = 14．计算下面矩阵的特征值、特征向量、迹和秩。

matlab习题二答案

matlab习题二答案Matlab习题二答案Matlab是一种强大的数值计算和科学计算软件，广泛应用于工程、科学、金融等领域。

在学习和使用Matlab的过程中，习题是一种很好的练习和巩固知识的方式。

本文将为大家提供一些Matlab习题二的答案，希望能对大家的学习有所帮助。

1. 编写一个函数，输入一个矩阵A，输出矩阵A的转置矩阵。

```matlabfunction B = transposeMatrix(A)B = A';end```2. 编写一个函数，输入一个矩阵A，输出矩阵A的每一行的平均值。

```matlabfunction avg = rowAverage(A)avg = mean(A, 2);end```3. 编写一个函数，输入一个矩阵A和一个标量k，输出矩阵A中大于k的元素个数。

```matlabfunction count = countGreaterThanK(A, k)count = sum(A(:) > k);end```4. 编写一个函数，输入一个矩阵A和一个标量k，输出矩阵A中大于k的元素的索引。

```matlabfunction indices = findGreaterThanK(A, k)indices = find(A > k);end```5. 编写一个函数，输入一个矩阵A，输出矩阵A的每一列的方差。

```matlabfunction variances = columnVariance(A)variances = var(A);end```6. 编写一个函数，输入一个矩阵A和一个标量k，将矩阵A中小于k的元素替换为0。

```matlabfunction B = replaceLessThanK(A, k)B = A;B(A < k) = 0;end```7. 编写一个函数，输入一个矩阵A和一个标量k，将矩阵A中小于k的元素替换为k。

```matlabfunction B = replaceLessThanK(A, k)B = A;B(A < k) = k;end```8. 编写一个函数，输入一个矩阵A和一个标量k，将矩阵A中大于k的元素替换为k。

matlab课后习题解答第二章

matlab课后习题解答第⼆章第2章符号运算习题2及解答1 说出以下四条指令产⽣的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+)〖⽬的〗不能从显⽰形式判断数据类型，⽽必须依靠class指令。

〖解答〗c1=3/7+c2=sym(3/7+c3=sym('3/7+')c4=vpa(sym(3/7+)Cs1=class(c1)Cs2=class(c2)Cs3=class(c3)Cs4=class(c4)c1 =c2 =37/70c3 =c4 =Cs1 =doubleCs2 =symCs3 =symCs4 =sym2 在不加专门指定的情况下，以下符号表达式中的哪⼀个变量被认为是⾃由符号变量.sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')〖⽬的〗理解⾃由符号变量的确认规则。

〖解答〗symvar(sym('sin(w*t)'),1)ans =wsymvar(sym('a*exp(-X)'),1)ans = asymvar(sym('z*exp(j*th)'),1) ans = z5求符号矩阵=333231232221131211a a a a a a a a a A 的⾏列式值和逆，所得结果应采⽤“⼦表达式置换”简洁化。

〖⽬的〗理解subexpr 指令。

〖解答〗A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')DA=det(A) IA=inv(A);[IAs,d]=subexpr(IA,d) A =[ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d* (a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)8（1）通过符号计算求t t y sin )(=的导数dtdy。

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 语言与应用》实验课程任务书一、实验教学目标与基本要求上机实验是本课程重要的实践教学环节；实验的目的不仅仅是验证理论知识，更重要的是通过上机实验，加强学生的实验手段与实践技能，掌握应用MATLAB 语言求解问题的方法，培养学生分析问题、解决问题、应用知识的能力和创新精神，全面提高学生的综合素质。

上机实验共8学时。

主要实验内容是基于理论课所学知识对课后典型习题进行MATLAB 求解，基本掌握常见数学问题的求解方法与命令调用，更深入地认识和了解MATLAB 语言强大的计算功能。

上机实验最终以书面报告的形式提交，并作为期末成绩考核内容的一部分。

二、实验内容（8学时）第一部分MATLAB 语言编程、科学绘图与基本数学问题求解（4学时）主要内容：掌握MATLAB 语言编程基础、科学绘图方法、微积分问题、线性代数问题等基本数学问题的求解与应用。

练习题：1、安装MATLAB 软件，应用demo 命令了解主要功能，熟悉基本功能，会用help 命令。

2、用MATLAB 语句输入矩阵A 和B⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1423143212344321A ， ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++++++++++++++++=4j 11j43j 22j34j 11j 42j 33j 24j 13j 22j 31j 41j 42j 33j 24j 1B 前面给出的是44⨯矩阵，如果给出5)6,5(=A 命令将得出什么结果？代码：A=[1 2 3 4;4 3 2 1;2 3 4 1;3 2 4 1]B=[1+4j 2+3j 3+2j 4+1j;4+1j 3+2j 2+3j 1+4j;2+3j 3+2j 4+1j 1+4j;3+2j 2+3j 4+1j 1+4j] A(6,5)=53、假设已知矩阵A，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给B矩阵，用magic(8)A 命令生成A矩阵，用上述命令检验一下结果是不是正确。

matlab习题二答案

习题二1.选择题(1)下列变量名中A是合法的。

A. char_1,i,jB. x*y,a.1C. x\y,a1234D. end,lbcx说明：end是关键字，变量名中不能有运算符、标点符号，可以有下划线。

内置常量可以作为变量名，但一般不提倡这样做。

(2)下列C是合法常量。

A. 3*e10B. 1e500C. -1.85e-56D. 10-2说明：10-2是表达式，1e500超过了realmax，3*e10也是表达式。

(3)x=uint8(2.3e10)，则x所占的字节是D个。

A. 1B. 2C. 4D. 8(4)已知x=0:10，则x有B个元素。

A. 10B. 11C. 9D. 12(5)产生对角线上为全1其余为0的2行3列矩阵的命令是C。

A. ones(2,3)B. ones(3,2)C. eye(2,3)D. eye(3,2)(6)已知数组123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，则a(:,end)是指C 。

A. 所有元素B. 第一行元素C. 第三列元素D. 第三行元素(7)已知数组123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，则运行a(:,1)=[]命令后C。

A. a变成行向量B. a数组为2行2列C. a数组为3行2列D. a数组中没有元素3(8)已知数组123456789a⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，则运行mean(a)命令是B。

A. 计算a每行的平均值B. 计算a每列的平均值C. a数组增加一行平均值D. a数组增加一列平均值(9)已知x为一个向量，计算ln(x)的MATLAB命令是计算B。

A. ln(x)B. log(x)C. Ln(x)D. lg10(x)(10)当a=2.4，使用取整函数计算得出3，则该取整函数名为C。

A. fixB. roundC. ceilD. floor(11)已知a=0:4，b=1:5，下面的运算表达式出错的为D。

A. a+bB. a./bC. a'*bD. a*b(12)命令day(now)是指C。

MatLab语言及应用--习题及参考答案

MatLab语言及应用--习题1.产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！1.A = [21 2 3; 4 35 6; 7 8 49];B = [3; 2; 1];设求AB、BTA、A-1与A-1B，求A2、A的2次数幂，并求矩阵A的秩、转置和行列式。

2.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5]4. 请修改下面的程序，让他们没有for循环语句！A=[1 2 3; 4 5 6; 7 8 9];[r c]=size(A);for i=1:1:rfor j=1:1:cif (A(i,j)>8 | A(i,j)<2)A(i,j)=0;endendend5.a=[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB中最简单的方法计算b，注意最简单。

6. 求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3x +8z- 5w=9-6x+5y-2z+10w=-8求系数矩阵的秩；求出方程组的解。

7.设 y=cos[0.5+((3sinx)/(1+x^2))]把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线；8. 设f(x)=x^5-4x^4+3x^2+2x+6取x=[-2,8]之间函数的值（取100个点），画出曲线，看它有几个零点。

9. 在[-10，10；-10，10]范围内画出函数的三维图形。

10.编写一个函数，使其能够产生如下的分段函数：⎪⎩⎪⎨⎧≤≤<-≤=x x x x x x f 65.0620.251.525.0)(，，，，并调用此函数，绘制曲线2)()(范围的2]，[0在+•+=x f x f x 。

（10分）11.将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。

MATLAB程序设计与应用第二版课后题答案2(可编辑修改word版)

第二章3.设矩阵A 为A=[24 23 9 21 6;65 74 24 11 21;34 5 98 75 21;8 42 42 53 121;43 21 45 64 21];（1）B=A(2:5,1:2:5)B =65 24 2134 98 218 42 12143 45 21（2）A(7)=[]A =24 65 34 8 43 23 5 42 21 9 24 98 42 45 21 11 75 53 64 6 21 21 121 21（3）A+30（4）size(A);ndims(A)（5）题目有误（6）reshape(x,3,4)（7）abs(x)（8）char(x)4. L1 =0 0 0 0 1 0 0 0 0L2 =1 1 1 1 1 0 0 0 0L3 =0 0 0 1 1 1 0 0 0L4 =4 5 65.(1)B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*CB =23.0000 10.0000 -0.7780 041.0000 -45.0000 65.0000 5.000032.0000 5.0000 0 32.0000C =23.0000 10.000041.0000 -45.000032.0000 5.00006.0000 -9.5400D =65.0000 5.00000 32.000054.0000 3.1400E =1.0e+003 *0.9141 -0.22391.20802.71231.1330 -0.2103(2)E<Dans =0 10 00 1E&Dans =1 10 11 1E|Dans =1 11 11 1~D|~Eans =0 01 00 0find(A>=10&A<25)ans =156.all(A)ans =any(A)ans =1isnan(A)ans =0 1 0 0 0 0 0isinf(A)ans =0 0 1 1 0 0 0isfinite(A)ans =1 0 0 0 1 1 17.A(1).x1=’学号’;A(1).x2=’姓名’;A(1).x3=’专业’;A(1).x4.x41=’成绩1’;………. A(2).x1=’学号’;A(2).x2=’姓名’;A(2).x3=’专业’;A(2).x4.x41=’成绩1’;………. A(3).x1=’学号’;A(3).x2=’姓名’;A(3).x3=’专业’;A(3).x4.x41=’成绩1’;………. A(4).x1=’学号’;A(4).x2=’姓名’;A(4).x3=’专业’;A(4).x4.x41=’成绩1’;………. A(5).x1=’学号’;A(5).x2=’姓名’;A(5).x3=’专业’;A(5).x4.x41=’成绩1’;……….8.(1)size(B)ans =2 2ndims(B)ans =2(2)B(2)ans =[3x3 doubleB(4)ans ={3x3 cell}(3)B(3)=[]B =[1][3x3 double] {3x3 cell}B{3}=[]B =[1] [3x3 double] []第三章1.(1)A=eye(3)(2) A=100+100*rand(5,6)(3)A=1+sqrt(0.2)*randn(10,50)(4)B=ones(size(A))(5)A+30*eye(size(A))(6)B=diag(diag(A))2.B=rot90(A)C=rot90(A,-1)3.B=inv(A) ;A 的逆矩阵C=det(A) ;A 的行列式的值D=A*BE=B*AD=E 因此A 与A-1 是互逆的。

MATLAB及应用参考答案

《MATLAB及应用》上机作业学院名称：机械工程学院专业班级：测控1201学生姓名：学生学号：201 年 4 月《MATLAB及应用》上机作业要求及规范一、作业提交方式：word文档打印后提交。

二、作业要求：1．封面：按要求填写学院、班级、姓名、学号，不要改变封面原有字体及大小。

2．内容：只需解答过程（结果为图形输出的可加上图形输出结果），不需原题目；为便于批阅，题与题之间应空出一行；每题答案只需直接将调试正确后的M文件内容复制到word 中（不要更改字体及大小），如下所示：%作业1_1clcA=[1 2 3 4;2 3 5 7;1 3 5 7;3 2 3 9;1 8 9 4];B=[1+4*i 4 3 6 7 7;2 3 3 5 5 4+2*i;2 6+7*i 5 3 4 2;1 8 9 5 4 3];C=A*BD=C(4:5,4:6)三、大作业评分标准：1．提交的打印文档是否符合要求；2．作业题的解答过程是否完整和正确；3．答辩过程中阐述是否清楚，问题是否回答正确；4．作业应独立完成，严禁直接拷贝别人的电子文档，发现雷同者都以无成绩论处。

作业11、用MATLAB 可以识别的格式输入下面两个矩阵12342357135732391894A ⎛⎫⎪⎪ ⎪= ⎪⎪ ⎪⎝⎭，144367723355422675342189543i i B i +⎛⎫⎪+⎪= ⎪+ ⎪⎪⎝⎭再求出它们的乘积矩阵C ，并将C 矩阵的右下角23⨯子矩阵赋给D 矩阵。

赋值完成后，调用相应的命令查看MATLAB 工作空间的占有情况。

解：A=[1 2 3 4;2 3 5 7;1 3 5 7 ;3 2 3 9 ;1 8 9 4;]B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] C=A*B D=C(4:5,4:6);2、设矩阵16231351110897612414152⎛⎫⎪⎪ ⎪ ⎪⎝⎭，求A ，1A -，3A ，12A A -+，1'3A A --，并求矩阵A 的特征值和特征向量。

MATLAB及其应用课后习题解答

2
5
“MATLAB 及其应用”课程作业
图 4- 2 画出函数图像
第5章
1．已知椭圆的长、短轴 a 4, b 2 ，用“小红点线”画如下图所示的椭圆
x a cos t 。 y b sin t
（提示：参量 t ；点的大小；axis equal）答：运行 :clf; a=4;b=2;t=0:pi/80:2*pi;x=a*cos(t);y=b*sin(t);plot(x,y,’r’,’markersize’,15);axis equal;xlabel(‘x’);ylabel(‘y’);shg;结果如图 5-1 所示：
1
“MATLAB 及其应用”课程作业
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))); da2 = vpa(a1-a2,30) da3 = vpa(a1-a3,30) da4 = vpa(a1-a4,30) da5 = vpa(a1-a5,30) da6 = vpa(a1-a6,30) da7 = vpa(a1-a7,30) da8 = vpa(a1-a8,30) 结果如图 2-2 所示：只有 a2 是精准的，观察可知，这是由于 sym('Num')和 sym(Num)的区别带来图 2- 2 vpa 观察误差的前者以字符串的形式传给符号运算内核，可以保留完整的精度；而后者经过浮点数运算之后再转换为符号类型，存在精度损失。 3．在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是独立自由变量。 sym('sin(w*t)') , sym('a*exp(-X)' ) , sym('z*exp(j*th)') 答：分别为 w,a,z，具体图 2-3 所示。