MATLAB习题参考答案
matlab课后习题答案(1-9章)
1 数字1.5e2，1.5e3 中的哪个与1500相同吗？1.5e32 请指出如下5个变量名中，哪些是合法的？abcd-2xyz_33chan a 变量ABCDefgh 2、5是合法的。
3 在MATLAB 环境中，比1大的最小数是多少？ 1+eps4 设 a = -8 , 运行以下三条指令，问运行结果相同吗？为什么？w1=a^(2/3) w2=(a^2)^(1/3) w3=(a^(1/3))^2w1 = -2.0000 + 3.4641i ;w2 = 4.0000 ;w3 =-2.0000 + 3.4641i 5 指令clear, clf, clc 各有什么用处？clear 清除工作空间中所有的变量。
clf 清除当前图形。
clc 清除命令窗口中所有显示。
第二章1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象？3/7+0.1双; sym(3/7+0.1)符; sym('3/7+0.1') 符;; vpa(sym(3/7+0.1)) 符;2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') symvar(sym('sin(w*t)'),1) w a z3 （1）试写出求三阶方程05.443=-x 正实根的程序。
注意：只要正实根，不要出现其他根。
（2）试求二阶方程022=+-a ax x 在0>a 时的根。
（1）reset(symengine)syms x positive solve(x^3-44.5) ans =(2^(2/3)*89^(1/3))/2（2）求五阶方程022=+-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 positivesolve(x^2-a*x+a^2) ans =a/2 + (3^(1/2)*a*i)/2 a/2 - (3^(1/2)*a*i)/24 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。
MATLAB-实用教程-课后习题标准答案
MATLAB-实用教程-课后习题标准答案第二章1.计算复数3+4i与5-6i的乘积。
a=3+4ib=5-6ic=a*b2.构建结构体Students，属性包含Name、age和Email，数据包括｛’Zhang’，18，[‘’，’’]｝、｛’Wang’，21，[]｝和｛’Li’，[],[]｝，构建后读取所有Name属性值，并且修改’Zhang’的Age 属性值为19。
Students(1).Age=18Students(1).Email='',''Students(2).Name='Wang'Students(2).Age=21Students(2).Email=[]Students(3).Name='Li'Students(3).Age=[]Students(3).Email=[]/doc/805835364.html,Student(1).Age(1)=19Student.Age3.用满矩阵和稀疏矩阵存储方式分别构造下属矩阵：A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0]A=[0 1 0 0 0;1 0 0 0 0;0 0 0 0 0;0 0 0 1 0]S=sparse(A)S=sparse([2,1,4],[1,2,4],[1,1,1],4,5)4.采用向量构造符得到向量[1,5,9....,41].A=1:4:415.按水平和竖直方向分别合并下述两个矩阵：A=[1 0 0;1 1 0;0 01],B=[2 3 4;5 6 7;8 9 10]A=[1 0 0;1 1 0;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]6.分别删除第五题两个结果的第2行。
A=[1 0 0;1 1 0;0 0 1]B=[2 3 4;5 6 7;8 9 10]C=[A B]D=[A;B]C(2,:)=[]D(2,:)=[]7.分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。
Matlab习题答案
参考答案： （1） >> (3-5*i)*(4+2*i) ans =
22.0000 -14.0000i
（2） >> sin(2-8*i) ans =
1.3553e+003 +6.2026e+002i
5．判断下面语句的运算结果。 (1) 4 < 20 (2) 4 <= 20 (3) 4 == 20 (4) 4 ~= 20 (5) 'b'<'B' 参考答案： （1） >> 4<20 ans =
y_nearest(i)=interp1(x,y,scalar_x(i),'nearest'); y_linear(i) =interp1(x,y,scalar_x(i),'linear'); y_spline(i) =interp1(x,y,scalar_x(i),'spline'); y_cubic(i) =interp1(x,y,scalar_x(i),'cubic'); end subplot(2,2,1),plot(x,y,'*'),hold on,plot(scalar_x,y_nearest),title('method=nearest'); subplot(2,2,2),plot(x,y,'*'),hold on,plot(scalar_x,y_linear),title('method=linear'); subplot(2,2,3),plot(x,y,'*'),hold on,plot(scalar_x,y_spline),title('method=spline'); subplot(2,2,4),plot(x,y,'*'),hold on,plot(scalar_x,y_cubic),title('method=cubic'); 得到结果为：
MATLAB习题参考答案(胡良剑,孙晓君)
解：建立 M 文件 pxy 如下： xa=-2:0.05:2;ya=xa; nx=length(xa);ny=length(ya); [x,y]=meshgrid(xa,ya); z=zeros(nx,ny); [a1,b1]=find(x+y>1); %第 a1 列 b1 行对应的 x+y>1 (x 对应列；y 对应行) %第 a1 列对应的 x 值是 xa(a1)；第 b1 行对应的 y 值是 ya(b1) z((a1-1)*ny+b1)=0.5457*exp(-0.75*ya(b1).^2-3.75*xa(a1).^2-1.5*xa(a1)); [a2,b2]=find(x+y<=1&x+y>-1); z((a2-1)*ny+b2)=0.7575*exp(-ya(b2).^2-6*xa(a2).^2); [a3,b3]=find(x+y<=-1); z((a3-1)*ny+b3)=0.5457*exp(-0.75*ya(b3).^2-3.75*xa(a3).^2+1.5*xa(a3)); surf(x,y,z); 命令窗口： >> pxy 运行结果如右图： 或者 M 文件如下： clear;close; xa=-2:0.1:2;ya=-2:0.1:2;[x,y]=meshgrid(xa,ya); z=zeros(size(x)); k1=find(x+y>1); z(k1)=0.5457*exp(-0.75*y(k1).^2-3.75*x(k1).^2-1.5*x(k1)); k2=find(x+y<=1&x+y>-1); z(k2)=0.7575*exp(-y(k2).^2-6*x(k2).^2); k3=find(x+y<-1); z(k3)=0.5457*exp(-0.75*y(k3).^2-3.75*x(k3).^2+1.5*x(k3)); mesh(x,y,z); 6、运行 demo 解：>>demo 7、查询 trapz 的功能、用法、目录、程序结构、相同目录下其它文件 解： >> help trapz ――功能用法 >> type trapz――程序结构，源码 >> which trapz――所在目录 >> help C:\MATLAB6p5\toolbox\matlab\datafun――该目录下其它文件
matlab课后习题答案 (附图)
习题2.1画出下列常见曲线的图形y （1）立方抛物线3x 命令：syms x y;ezplot('x.^(1/3)')(2)高斯曲线y=e^(-X^2); 命令：clearsyms x y;ezplot('exp(-x*x)')(3)笛卡尔曲线命令：>> clear>> syms x y;>> a=1; >> ezplot(x^3+y^3-3*a*x*y)(4)蔓叶线命令：>> clear>> syms x y;>> a=1ezplot(y^2-(x^3)/(a-x))(5)摆线:()()t b y t t a x cos 1,sin -=-=命令：>> clear>> t=0:0.1:2*pi;>> x=t-sin(t);>>y=2*(1-cos(t));>> plot(x,y)7螺旋线命令：>> clear >> t=0:0.1:2*pi; >> x=cos(t); >> y=sin(t); >> z=t;>>plot3(x,y,z)(8)阿基米德螺线命令：clear>> theta=0:0.1:2*pi;>> rho1=(theta);>> subplot(1,2,1),polar(theta,rho1)(9) 对数螺线命令：cleartheta=0:0.1:2*pi;rho1=exp(theta);subplot(1,2,1),polar(theta,rho1)(12)心形线命令：>> clear>> theta=0:0.1:2*pi;>> rho1=1+cos(theta);>> subplot(1,2,1),polar(theta,rho1)练习2.21. 求出下列极限值（1）nn n n 3lim 3+∞→命令：>>syms n>>limit((n^3+3^n)^(1/n))ans = 3（2）)121(lim n n n n ++-+∞→命令：>>syms n>>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf)ans =（3）x x x 2cot lim 0→命令：syms x ；>> limit(x*cot(2*x),x,0)ans =1/2（4）)(cos lim c m xx ∞→ 命令：syms x m ；limit((cos(m/x))^x,x,inf)ans =1（5）)111(lim 1--→e x x x 命令：syms x>> limit(1/x-1/(exp(x)-1),x,1)ans =(exp(1)-2)/(exp(1)-1)（6）)(2lim x x x x -+∞→命令：syms x>> limit((x^2+x)^(1/2)-x,x,inf)ans =1/2练习2.41. 求下列不定积分，并用diff 验证：（1）⎰+xdx cos 1 >>Clear>> syms x y>> y=1/(1+cos(x));>> f=int(y,x)f =tan(1/2*x)>> y=tan(1/2*x);>> yx=diff(y,x);>> y1=simple(yx)y1 =1/2+1/2*tan(1/2*x)^2 (2)⎰+e x dx 1clearsyms x yy=1/(1+exp(x));f=int(y,x)f =-log(1+exp(x))+log(exp(x))syms x yy=-log(1+exp(x))+log(exp(x));yx=diff(y,x);y1=simple(yx)y1 =1/(1+exp(x))(3)dx x x ⎰sin 2syms x yy=x*sin(x)^2;>> f=int(y,x)f =x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2clearsyms x y y=x*(-1/2*cos(x)*sin(x)+1/2*x)-1/4*cos(x)^2-1/4*x^2;yx=diff(y,x);>> y1=simple(yx)y1 =x*sin(x)^2(4) xdx ⎰sec 3syms x yy=sec(x)^3;f=int(y,x)f =1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x))clearsyms x yy=1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x));yx=diff(y,x);y1=simple(yx)y1 =1/cos(x)^32. 求下列积分的数值解1）dx x x ⎰-10clearsyms xy=int(x^(-x),x,0,1)y =int(x^(-x),x = 0 .. 1)vpa(y,10)ans =1.2912859972）xdx e x cos 3202⎰πclearsyms xy=int(exp(2*x)*cos(x)^3,x, clearsyms xy=int((1/(2*pi)^(1/2))*exp(-x^2/2),x,0,1)y =7186705221432913/36028797018963968*erf(1/2*2^(1/2))*2^(1/2)*pi^(1/0,2*pi) y =22/65*exp(pi)^4-22/65vpa(ans,10)(3) dx x e 210221-⎰π>> clear>> syms x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),0,1);>> vpa(y,14)ans =.341344746068552（4）>> clear>> syms x>> y=int(x*log(x^4)*asin(1/x^2),1,3);Warning: Explicit integral could not be found.> In sym.int at 58>> vpa(y,14)ans =2.45977212823752（5）>> clear>> syms x>> y=int(1/(2*pi)^(1/2)*exp(-x^2/2),-inf,inf);>> vpa(y,14)ans =.99999999999999练习2.51判断下列级数的收敛性，若收敛，求出其收敛值。
matlab教材习题答案
matlab教材习题答案Matlab是一种广泛应用于科学与工程领域的计算机编程语言和环境。
它具备强大的数值计算和数据可视化功能，被广泛用于数据分析、信号处理、图像处理、机器学习等领域。
对于初学者而言，掌握Matlab的基本语法和常用函数是非常重要的，而教材习题则是帮助学生巩固所学知识的重要资源。
本文将为大家提供一些Matlab教材习题的参考答案，以帮助读者更好地学习和应用Matlab。
1. 基本语法练习题1.1 计算并输出1到10的平方for i = 1:10fprintf('%d的平方是：%d\n', i, i^2);end1.2 计算并输出1到10的阶乘for i = 1:10fact = 1;for j = 1:ifact = fact * j;endfprintf('%d的阶乘是：%d\n', i, fact);end2. 数值计算练习题2.1 求解一元二次方程的根a = 1;b = -3;c = 2;delta = b^2 - 4*a*c;x1 = (-b + sqrt(delta))/(2*a);x2 = (-b - sqrt(delta))/(2*a);fprintf('一元二次方程的根为：%f, %f\n', x1, x2);2.2 求解线性方程组的解A = [1 2; 3 4];B = [5; 6];X = inv(A) * B;fprintf('线性方程组的解为：%f, %f\n', X(1), X(2));3. 数据处理练习题3.1 统计一个数组中的最大值、最小值和平均值data = [1, 2, 3, 4, 5];max_value = max(data);min_value = min(data);average_value = mean(data);fprintf('最大值：%f\n最小值：%f\n平均值：%f\n', max_value, min_value, average_value);3.2 对一个矩阵进行排序matrix = [4 2 3; 1 5 6; 9 8 7];sorted_matrix = sort(matrix);fprintf('排序后的矩阵为：\n');disp(sorted_matrix);4. 图像处理练习题4.1 读取并显示一张图片image = imread('image.jpg');imshow(image);4.2 对一张图片进行灰度化处理gray_image = rgb2gray(image);imshow(gray_image);5. 信号处理练习题5.1 生成并绘制正弦信号t = 0:0.01:2*pi;x = sin(t);plot(t, x);5.2 对一段音频信号进行傅里叶变换[y, fs] = audioread('audio.wav');Y = fft(y);plot(abs(Y));通过以上几个例子，我们可以看到Matlab的强大功能和灵活性。
matlab课后习题答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答：A(7)=[](3)将矩阵A的每个元素值加30。
答：A=A+30;(4)求矩阵A的大小和维数。
答：size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成矩阵。
34答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8)求一个ASCII 码所对应的字符。
答：char(49);4.下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5.已知23100.7780414565532503269.5454 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦完成下列操作：(1)取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角子矩阵构32⨯成矩阵D ，B 与C 的乘积构成矩阵E 。
Matlab参考答案
Matlab参考答案参考答案一、填空题（每空1分，20分）1、MATLAB常用操作界面包括命令窗口、工作空间窗口（浏览器）、命令历史窗口、当前目录窗口、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLABInf或inf表示正无穷大、NaN或nan表示非数、nargout表示函数输出参数。
3、MATLAB中逗号主要用作要显示计算结果的指令与其后指令的分隔；用作输入量与输入量之间的分隔符；用作数组元素分隔符号。
4、工作空间浏览器主要用于内存变量的查阅、保存和编辑。
5、MATLAB实现将全下标转换为单下标的指令为Sub2ind、据单下标换算出全下标的指令为Ind2sub。
6、二维数组的标识有“全下标”标识、“单下标”标识、“逻辑1”标识。
7、在一个元胞数组A中寻访第2行第3列元胞元素用A(2,3)；寻访数组第2行第3列元胞中的内容用A{2,3}。
8、MATLAB中clf用于清除图形窗、clc用于清除指令窗中显示内容、clear用于清除MATLAB工作空间中保存的变量。
二、简答题（每题5分，共20分）1、简述MATLAB历史指令窗的主要作用。
历史指令窗记录着用户在MATLAB指令窗中所输入过的所有指令。
历史记录包括：每次开启MATLAB的时间，每次开启MATLAB后在指令窗中运行过的所有指令。
应用功能有单行或多行指令的复制和运行、生成M文件等。
2、简述空数组的功用。
“空”数组的功用：在没有“空”数组参与的运算时，计算结果中的“空”可以合理地解释“所得结果的含义”；运用“空”数组对其他非空数组赋值，可以改变数组的大小，但不能改变数组的维数。
3、简述MATLAB函数的基本结构。
典型M函数文件的结构：函数申明（定义）行(Function declaration line)、H1行(The first help text line)、在线帮助文本(Help text)区、编写和修改记录、函数体(Function body)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第1章 MATLAB简介1、MA TLAB的主要特点有：①语言简洁，编程效率高。
②人机界面友善，交互性好。
③绘图功能强大，便于数据可视化。
④学科众多、领域广泛的MATLAB工具箱。
⑤源程序的开放性。
MATLAB的典型应用领域有：①自动控制②汽车③电子④仪器仪表⑤生物医学⑥信号处理⑦通信等。
2、填空题⑴命令窗口、命令历史窗口、当前目录窗口⑵查阅、保存、编辑⑶清除图形窗、清除命令窗口中显示内容、清除MATLAB工作空间中保存的变量。
3、如果想查看某一变量具体内容或者对其修改操作，可以在工作空间中双击该变量名称，可以打开数组编辑器，在数组编辑器中可以查看变量的具体内容，也可以对其修改。
如果想删除MATLAB内存中的变量，可以在工作空间中选中该变量，然后利用工作空间窗口的菜单命令或工具条中的快捷图标进行删除。
4、1+2+3+4+5+... （+ 后面可以直接跟...，也可以在+和...中加一个空格。
）6+7+8+91+2+3+4+5 ... （5后面必须跟一个空格，不能直接跟...，否则报错，这在预置一个+6+7+8+9 大数组时很重要。
）第2章矩阵与数值数组1、填空题：⑴非数、无穷大、机器零阈值，浮点数相对精度,eps= 2.2204e-016。
⑵全下标、单下标。
2、阅读程序题：（本题主要考察数组的寻访、赋值和简单运算，提示：带；的语句不显示结果）⑴ans =2 3 7Sa =102030A =1 20 5 30 910 4 6 8 10⑵ A =1 3 5 7 9 11 13 152 4 6 8 10 12 14 16ans =1 5 9 132 6 10 143 7 11 154 8 12 16A =0 0 5 7 0 0 13 152 4 0 0 10 12 0 0⑶ ans =-1 -46 4ans =3 05 -2ans =3 69 12ans =3 69 123、A=magic(4)；L=A<10L =0 1 1 01 0 0 11 1 1 01 0 0 14、省略。
第3章字符串、元胞和构架数组1. 直接创建法：S1=['Where there is life,''there is hope. '] %注意第2行要加入6个空格使其与第一行字符数相等，否则报错。
S2=str2mat('Where there is life,','there is hope.')S3=strvcat('Where there is life,','there is hope.') %后两种方法则不用考虑两行字符的数目2. 填空题：A（2，3）；A{2，3}3. 阅读程序题：⑴subch =ABc12revch =4d321cBA⑵ （提示：本题考察元胞数组两种不同的寻访方法） A{1,1} = this is cell A{2,1} =1.0000 +2.0000i A{1,2}{1} = 1 2 3 4 A{2,2} = 6 7 8 94. % Homework3.4student=cell(2,2); %利用 cell 函数预先分配存储空间 student(1,1)={['20090102001';'20090102002']}; student(2,1)={{'张三';'李小红'}};student(1,2)={{'数学1' '英语1' '代数' '制图';'数学2' '英语2' '电路' '信号'}}; student(2,2)={{[80 90 65 96;92 75 80 78];[78 69 80 74;50 78 68 80]}}; celldisp(student)5. % Homework3_5.mstudent(1)=struct('class',{'电子技术'},'name',{'无敌'},'gender', {'男'},'number',{'200801001'}, ...'course',{{'英语' '数学' '计算机' }}, 'score', {[75 85 67]}) student(2)=struct('class',{'电子技术'},'name',{'王思远'},'gender', {'男'},'number',{'200801002'}, ... 'course',{{'英语' '数学' '计算机' }}, 'score', {[99 78 0]}) student(3)=struct('class',{'通信工程'},'name',{'李蕾'},'gender', {'女'},'number',{'200802001'}, ...'course',{{'英语' '数学' '计算机' }}, 'score', {[80 72 53]}) student(4)=struct('class',{'通信工程'},'name',{'王伟'},'gender', {'男'},'number',{'200802002'}, ... 'course',{{'英语' '数学' '计算机' }}, 'score', {[97 82 62]})第4章 数值计算1．rand(2,3)的含义为 生成2×3随机矩阵 。
2.应用LU 分解和矩阵除法求方程组⎪⎩⎪⎨⎧=+=+-=-10421554832z x z y x y x 的解。
方法一：A=[2,-3,0;4,-5,1;2,0,4]; [L,U]=lu(A);B=[8;15;10]; x=U\(L\B)方法二：A=[2,-3,0;4,-5,1;2,0,4]; B=[8;15;10]; A\B3．求一元非线性方程0sin =-xxx 的根。
程序：%计算函数数值 x=[-10:0.1:10]; y=x-sin(x)./x; %绘制函数图形 plot(x,y);%锁定当前图形并添加网格线 %hold on; grid on% 绘制水平线line([-10 10],[0 0]); %添加坐标名称xlabel('x');ylabel('y(x)');[xx,yy]=ginput(1) %在MATLAB 指令窗中运行，用鼠标获2个零点猜测值。
%第三步，计算xx(1) 和xx(2)的精确零点y=inline('x-sin(x)./x','x'); %构造内联函数求解 [x,y]=fzero(y,xx)运行结果：-10-8-6-4-20246810-10-551015xy (x )xx =0.8525 yy =-0.0219 x =0.8767 y = 04．求f(x)=3x 2-5x ＋1极小值。
%绘制函数曲线 x=[-5:0.1:5]; f=3*x.^2-5.*x+1; plot(x,f);%估计极小值点的坐标。
由图可以看出，在[0 2]间有极小值 [xx,ff]=ginput(1)-5-4-3-2-1012345-20020406080100120%编写函数文件，保存文件名为myfun4_4.mfunction f=myfun4_4(x) f=3*x.^2-5.*x+1;%求极小值>> [x,fval] = fminbnd (@myfun4_4,0,2) x =0.8333 fval =-1.08335．求积分⎰⎰+=πππ22cos(dxdy y x x QQ = dblquad(@(x,y) (x.*cos(x+y^2)), pi, 2*pi, 0, pi)7．求常微分方程的数值解。
⎪⎪⎩⎪⎪⎨⎧-='==++-1)0(1)0(02)1(22y y y dx dyx dx y d dsolve('D2y*(x-1)+2*Dy+y=0','y(0)=1,Dy(0)=1','x')第5章符号运算5.1 C5.2 B5.3 C5.4syms x;f=x^4-5*x^3+5*x^2-6;factor(f)5.5 syms x;f=cos(x)+sqrt(-sin(x)^2);f1=simple(f)f2=simple(f1)5.6 syms a t x;f=[a*x^2,2*t^4;exp(t)*sin(x),t*log(x)];dfdt2=diff(f,t,2)dfdxdt=diff(diff(f,x),t)5.7 syms x y z;f=int(int(int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2),x,1,2) vpa(f)5.8 A=sym([1 1/2 1/2 -1;1 1 -1 1;1 -1/4 -1 1;-8 -1 1 1]);B=sym([0;10;0;1]);X=A\B5.9 S=dsolve('D2x+2*Dy-x=0,Dx+y=0','x(0)=1,y(0)=0');S.xS.y5.10 syms t s;syms a b positive; % a,b为正Dt=sym('Dirac(t)');Ut=sym('Heaviside(t-a)');T=[Dt,Ut;t^2*sin(2*t),exp(-a*t)*cos(b*t)];S=laplace(T,t,s)第6章 MATLAB图形绘制6.1 C6.2 A6.3 D6.4 x=0:0.01:6;y1=sin(2*x);y2=x.*cos(x);y3=exp(2*cos(x));plot(x,y1,'r:',x,y2,'b-.',x,y3,'g--')6.5 x=-10:0.15:10;y=sin(x)./x;plot(x,y)xlabel('x');ylabel('y=sinx/x');title('门函数的频谱')6.6 x=[30.5,26.5,21,10,12];subplot(1,2,1)pie(x,[1 1 0 0 0])subplot(1,2,2)pie3(x,[1 0 1 0 0])legend({'1','2','3','4','5'})6.7 income=[1.02 1.35 1.66 2.01];x=2005:2008;expenses=[0.85 1.06 1.21 1.48];area(x,income,'facecolor',[0.5 0.9 0.6],'edgecolor',...'b','linewidth',2) % 设置填充色、边界色和边界宽度hold onarea(x,expenses,'facecolor',[0.9 0.8 0.7],'edgecolor','r','linewidth',2)hold offset(gca,'xtick',[2005:2008]) % x轴为2005到2008xlabel('年份','fontsize',10) % x轴标注字号为10号ylabel('万元','fontsize',10)gtext('支出') % 通过光标在图形上添加注释gtext('收入')gtext('\rightarrow利润')6.8 x=[0.5 1.5 3];y=[2 5 6 4 13 4 5 2 21 3 82 1];bar3(x,y,'detached')xlabel('x轴')ylabel('y轴')zlabel('z轴')6.9略.6.10略.第7章 M文件和MATLAB编程1.脚本文件的扩展名为m ，函数文件的扩展名为m 。