matlab课后习题解答第二章
MATLAB程序设计与应用第二版课后题答案2(最新整理)
-0.6863 0.5621 0.4615
-0.0937 -0.6976 0.7103
D=
-0.0166
0
0
0 1.4801
0
0
0 2.5365
第四章
1.a=input('请输入一个 4 位数:');
while (a<1000|a>9999)
a=input('输入错误，请重新输入一个4位数:');
end
9
function f=factor(n,m); y=0; for k=1:n
y=y+k^m; end
10.(1)S=108 (2)x=4 12 20
y=2 4 6
1. (1) x=-10:0.1:10; y=100./(1+x.^2); plot(x,y)
第五章
(2) x=-10:0.1:10; y=1/(2*pi)*exp(-x.^2/2); plot(x,y)
第二章 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 21 34 98 21 8 42 121 43 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 = 000010000 L2 =
精讲多练MATLAB(第二版)课后习题答案
第1章练习题答案1-2.A=1.2 ;B=-4.6 ;C=8.0 ;D=3.5 ;E=-4.0;T=atan((2*pi*A+E/(2*pi*B*C))/D)1-3clearx=45;(sin(x*pi/180)+sqrt(35))/(72^(1/5))1-4cleara=5.67;b=7.811;exp(a+b)/log10(a+b)1-5clearx=3;y=sqrt(x)-6*(x+1/x)+(x-3.2)^2/(x+7.7)^3 1-6clearr=15;d=2*rc=2*pi*rs=pi*r*r1-7cleara=8.5;b=14.6;c=18.4;s=(a+b+c)/2;area=sqrt(s*(s-a)*(s-b)*(s-c))第2章练习题答案2-1 c learA=[ 3 1 1;2 1 2;1 2 3];B=[1 1 -1;2 -1 0;1 -1 1];2*A+B,4*(A^2)-3*(B^2),A*B,B*A,A*B-B*A2-2 clearA=[1/3 0 0;0 1/4 0;0 0 1/7];B=inv(inv(A)-eye(3))*62-3 clearB=[1 2 -3 -2;0 1 2 -3;0 0 1 2;0 0 0 1]; C=[1 2 0 1;0 1 2 0;0 0 1 2;0 0 0 1];A=(inv(2*eye(4)-inv(C)*B)*inv(C))'2-4 clearA=[2 -1;-1 2];B=[0 -2;-2 0];X=(1/2)*(B+2*A)2-5 clearA=[2 -3 0 2;1 5 2 1;3 -1 1 -1;4 1 2 2];B=[8;2;7;12];X=A\B2-6 A=[3 12 4 7 0 8 1];p=poly2str(A,'x')x=roots(A)2-7 A=[3 12 4 7 0 8 1];B=[1 -3 5 -15 0];[div,rest]=deconv(A,B)第三章3-1.x=sym('x');limit((cos(sqrt(x)))^(pi/x),x,0,'right')ans =exp(-1/2*pi)3-2.x=sym('x');f=sym('(3*sin(x)+x^2*cos(1/x))/((1+cos(x))*log(1+x))');limit(f,x,0)ans =3/23-3.x=sym('x');f=sym('(sqrt(4*x^2+x-1)+x+1)/sqrt(x^2+sin(x))');limit(f,x,-inf)ans =13-4.s yms x y;f=sym('(x^2+y^2)^(x^2*y^2)');limit(limit(f,x,0),y,0)ans =13-5.x=sym('x');y=sym('(tan(sqrt(x+sqrt(x+sqrt(2*x)))))^2');diff(y,x)ans=tan((x+(x+2^(1/2)*x^(1/2))^(1/2))^(1/2))*(1+tan((x+(x+2^(1/2)*x^(1/2 ))^(1/2))^(1/2))^2)/(x+(x+2^(1/2)*x^(1/2))^(1/2))^(1/2)*(1+1/2/(x+2^( 1/2)*x^(1/2))^(1/2)*(1+1/2*2^(1/2)/x^(1/2)))3-6.x=sym('x');y=sym('cos(x^2)*sin(1/x)^2');diff(y,x)ans =-2*sin(x^2)*x*sin(1/x)^2-2*cos(x^2)*sin(1/x)*cos(1/x)/x^23-7.x=sym('x');y=sym('sqrt(sin(x)-sin(x)^3)');int(y,x,0,pi)ans =4/33-8.x=sym('x');y=sym('sqrt((x+1)/(x-1))/x');int(y,x)ans=-((1+x)/(x-1))^(1/2)*(x-1)*(atan(1/(-1+x^2)^(1/2))-log(x+(-1+x^2)^(1/2)))/((1+x)*(x-1))^(1/2)3-9.dsolve('D2y+4*Dy+4*y=exp(-2*x)','x')ans =1/2*exp(-2*x)*(2*C2+2*x*C1+x^2)3-10.dsolve('x^2*Dy+x*y=y^2','y(1)=1')ans =x*exp(1/x)/(exp(1/x)+exp(1/x*t)*x-exp(1/x*t))第四章4-1.plot([1 2 4 0 5 10 11 21 3 1])4-2.t=0:pi/100:2*pi;x=cos(t);y=sin(t);plot(x,y);axis('square');4-3.R=[1 5 10 20];I=0:0.1:10;U=I'*R;plot(I,U);4-4.mon=1:12;temp=[0.2 2.3 8.7 18.5 24.6 32.1 36.8 37.1 28.3 17.8 6.4 -3.2];rainf=[4.6 3.6 2.1 2.9 3.0 2.7 2.2 2.5 4.3 3.4 2.1 3.7];plot(mon,temp,'r-p',mon,rainf,'b:s');xlabel('month');ylabel('temperature and rainfall');%title('Temperature and Rainfall Curve');%text(3,1,'rainfall');%text(9,28.3,'temperature');flag=12;while(flag)text(mon(flag),temp(flag),[num2str(mon(flag)),num2str(temp(flag))]); text(mon(flag),rainf(flag),[num2str(mon(flag)),num2str(rainf(flag))]);flag=flag-1;end4-5.A=[1 1 1 1 1 1 1;1 2 2 2 2 2 1;1 2 2 3 2 2 1;1 2 2 2 2 2 1;1 1 1 1 1 1 1];plot(A);4-6.x=-10:0.5:10;y=x;[X,Y]=meshgrid(x,y);Z=X.^2+Y.^2;surf(X,Y,Z);4-7.clccleartheta=(0:0.1:2)*pi;phi=(0:0.1:2)*pi;x=zeros(length(theta),length(phi));y=x;z=x;for i=1:length(theta)for j=1:length(phi)x(i,j)=2*sin(theta(i))*cos(phi(j));y(i,j)=3*sin(theta(i))*sin(phi(j));endz(i,:)=4*cos(theta(i));endaxis squaresubplot(2,2,1);mesh(x,y,z)title('三维')subplot(2,2,2);mesh(x,y,z)view(90,0)title('X轴视图')subplot(2,2,3);mesh(x,y,z)view(0,0)title('Y轴视图')subplot(2,2,4);mesh(x,y,z)view(0,90)title('Z轴视图')第五章5-1：clear;x=-2:0.01:2;for ii=1:401if(x(ii)<0)y(ii)=x(ii)+1;elseif (x(ii)<1)&(x(ii)>=0)y(ii)=1;elseif(x(ii)>=1)y(ii)=x(ii)^3;endendplot(x,y);5-2：function y=fenduan(x)if(x<0)y=x+1;elseif (x<1)&(x>=0)y=1;elseif(x>=1)y=x^3;end5-3：建立文件arcsin.m如下：function y=arcsin(x)if abs(x)>1disp('输入参数的绝对值必须小于1');return;endn=0;y=0;u=x;while u>epsu=(jiecheng(2*n)*x^(2*n+1))/(2^(2*n)*(jiecheng(n)^2*(2*n+1))); y=y+u;n=n+1;end建立文件jiecheng.m如下：function y=jiecheng(n)if(n<0)disp('输入参数必须为大于或等于0的整数');return;endif n==0 | n==1y=1;elsey=n*jiecheng(n-1);end5-4：建立shuliang.m文件如下：function [tu,ji]=shuliang(tou,jiao)tu=jiao/2-tou;ji=2*tou-jiao/2;在命令窗口中执行：>>[a,b]=shuliang(36,100)5-5：for i=2:999ge=rem(i,10);shi=rem(fix(i/10),10);bai=fix(i/100);if rem((ge+shi+bai),2) & isprime(i)disp(i);endend5-6：function y=issushu(n)y=1;for ii=2:fix(n/2)if(rem(n,ii)==0)y=0;endend习题1（程序段后面的是其运行的结果）例题1.2function lt12 a=1+2*i;b=3-4*i;disp(['a+b=',num2str(a+b)]) disp(['a-b=',num2str(a-b)]) disp(['a*b=',num2str(a*b)]) disp(['a/b=',num2str(a/b)]) >> lt12 a+b=4-2i a-b=-2+6i a*b=11+2ia/b=-0.2+0.4i >>例题1.3function z=lt13x=(-3.5)*pi/180;y=6.7*pi/180;z=sin(abs(x)+abs(y))/sqrt(cos(abs(x+y))); >> lt13ans =0.1772 >>例题1.4function p1=lt14r=0.0107;n=2015-2000;p0=12.9533E8; p1=p0*(1.0+r)^n; >> lt14ans =1.5195e+009 >>例题1.5function lt15 a=1;b=2;c=3;d=sqrt(b^2-4*a*c);x1=(-b+d)/2*a;x2=(-b-d)/2*a; disp('原方程的根是：')disp(['x1=',num2str(x1)]) disp(['x2=',num2str(x2)])>> lt15原方程的根是：x1=-1+1.4142i x2=-1-1.4142i >>。
matlab 部分答案
4. 用三点公式计算 y=f(x)在 x=1.0,1.2 处的导数值， f(x)值由下表给出； x f(x)
5
1.0 0.25
4
1.1
ww w
1.2
2
.k
1.3 0.1890 1.4 0.1736
后
7. 求解下列线性常微分方程的解析解。
答
6．设方程的根为 x = [−3, −5, −8, −9] ，求它们对应的 x 多项式的系数。
y=x(:,1);
20． 求解线性常微分方程 3 y + 4 y + 5 y + 6 y = 3u + 0.5u + 4u ， 在输入 u(t)为单位脉冲
''' '' ' '' '
并单位阶跃信号时的解析解。
第四章 习题
1．若 x(n) = cos(
课
2．求有限长序列 x( n) = 5(0.6) ( 0 ≤ n < 20 )的圆周移位 f (n) = x((n − 10)) 20 R20 ( n) 。
5． 创建 3×4 矩阵魔方阵和相应的随机矩阵，将两个矩阵并接起来，然后提取任意两个列向
6． 创建一个 4×4 单位阵，提取对角线以上部分。
8． 创建一个 5×5 随机阵并求其逆。 9． 利用上题的矩阵，计算矩阵的 5 次方。
求 C = A * B ， D = A.* B 。
B=[5,4,3,-2; 6,-2,3,-8; -1,3,-9,7]; C=A*B
部分习题与解答
第二章 习题及部分解答
1 1． 计算 y = x3 + ( x − 0.98) 2 /( x + 1.35)3 − 5( x + ) ，当 x = 2 和 x = 4 时的值。 x
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）第2章习题答案
7.平均值说明f (x ，y )的平均值等于其傅里叶变换F (u ，v )在频率原点的值F (0,0)。
2-3证明离散傅里叶变换的频率位移和空间位移性质。
证明：)(2101),(1),(NvyM ux j M x N y e y x f MN v u F +--=-=∑∑=π),(),(1),(100)(21010)(2)(21010000v v u u F dxdy ey x f MNe ey x f MN y Nv v x M u u j M x N y N yv M x u j Nvy M ux j M x N y --==-+---=-=++--=-=∑∑∑∑πππ因为()()v u F y x f ,,⇔ 所以 ),(),(00)(200v v u u F e y x f N y v M x u j --⇔+π2-4小波变换是如何定义的？小波分析的主要优点是什么？小波之所以小，是因为它有衰减性，即是局部非零的；而称为波，则是因为它有波动性，即其取值呈正负相间的振荡形式，将)(2R L 空间的任意函数f (t )在小波基下展开，称其为函数f (t )的连续小波变换。
小波变换是时间（空间）频率的局部化分析，它通过伸缩平移运算对信号（函数）逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频信号的要求从而可聚焦到信号的任意细节，解决了Fourier 变换的困难问题。
2-5 在图像缩放中，采用最近邻域法进行放大时，如果放大倍数太大，可能会出现马赛克效应，这个问题有没有办法解决，或者有所改善。
可以利用线性插值法，当求出的分数地址与像素点不一致时，求出周围四个像素点的距离比，根据该比率， 由四个邻域的像素灰度值进行线性插值。
2-6 复合变换的矩阵等于基本变换的矩阵按顺序依次相乘得到的组合矩阵。
即，T=T N T N-1…T 1。
问矩阵顺序的改变能否影响变换的结果。
矩阵顺序的改变不会影响变换的结果。
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教程第二章课堂练习及答案
Matlab教程第二章符号计算课堂练习1 创建符号变量有几种方法？MATLAB提供了两种创建符号变量和表达式的函数：sym和syms。
sym用于创建一个符号变量或表达式，用法如x=sym(‘x’) 及f=sym(‘x+y+z’)，syms用于创建多个符号变量，用法如syms x y z。
f=sym(‘x+y+z’)相当于syms x y zf= x+y+z2 下面三种表示方法有什么不同的含义？（1）f=3*x^2+5*x+2（2）f='3*x^2+5*x+2'（3）x=sym('x')f=3*x^2+5*x+2（1）f=3*x^2+5*x+2表示在给定x时，将3*x^2+5*x+2的数值运算结果赋值给变量f，如果没有给定x则指示错误信息。
（2）f='3*x^2+5*x+2'表示将字符串'3*x^2+5*x+2'赋值给字符变量f，没有任何计算含义，因此也不对字符串中的内容做任何分析。
（3）x=sym('x')f=3*x^2+5*x+2表示x是一个符号变量，因此算式f=3*x^2+5*x+2就具有了符号函数的意义，f也自然成为符号变量了。
3 用符号函数法求解方程a t2+b*t+c=0。
>> r=solve('a*t^2+b*t+c=0','t')[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]4 用符号计算验证三角等式：sin(ϕ1)cos(ϕ2)-cos(ϕ1)sin(ϕ2) =sin(ϕ1-ϕ2) >> syms phi1 phi2;>> y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2)) y =sin(phi1-phi2)5 求矩阵⎥⎦⎤⎢⎣⎡=22211211a a a a A 的行列式值、逆和特征根。
MATLAB运算基础(第2章)答案
实验01讲评、参考答案讲评未交实验报告的同学名单数学：6人（11、12级）信科：12-04, 12-22, 13-47批改情况：问题1：不仔细，式子中出错。
问题2：提交的过程不完整。
问题3：使用语句尾分号(;)不当，提交的过程中不该显示的结果显示。
问题4：截屏窗口没有调整大小。
附参考答案：《MATLAB软件》课内实验王平实验01 MATLAB运算基础（第2章MATLAB数据及其运算）一、实验目的1. 熟悉启动和退出MATLAB 的方法。
2. 熟悉MATLAB 命令窗口的组成。
3. 掌握建立矩阵的方法。
4. 掌握MATLAB 各种表达式的书写规则以及常用函数的使用。
二、实验内容1. 数学表达式计算先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
1.1 计算三角函数122sin 851z e=+（注意：度要转换成弧度，e 2如何给出） 示例：点击Command Window 窗口右上角的，将命令窗口提出来成悬浮窗口，适当调整窗口大小。
命令窗口中的执行过程：1.2 计算自然对数221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦（提示：clc 命令擦除命令窗口，clear 则清除工作空间中的所有变量，使用时注意区别，慎用clear 命令。
应用点乘方） 命令窗口中的执行过程：1.3 求数学表达式的一组值0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。
命令窗口中的执行过程：1.4 求分段函数的一组值2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5提示：用逻辑表达式求分段函数值。
命令窗口中的执行过程：1.5 对工作空间的操作接着显示MATLAB当前工作空间的使用情况并保存全部变量提示：用到命令who, whos, save, clear, load，请参考教材相关内容。
matlab课后习题答案
第2章 MATLAB 矩阵运算基础2.1 在MATLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
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。
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文件中。
matlab-第二章
1，在一个MA TLAB命令中，6+7i和6+7*i有何区别？i和I有何区别？第一个i代表一个变量,6+7i是一个函数,第二个代表虚数,6是实部,7代表虚步.I只代表变量.2，.设A和B是两个同大小的矩阵，试分析A*B和A.*B、A./B和B.\A、A/B和B/A的区别？如果A和B是两个标量数据，结论又如何？A *B是矩阵相乘，A.*B表示A和B单个元素之间对应相乘。
A./B和B.\A值相等，是一样的。
A/B等效于B的逆右乘A的矩阵，B\A等效于B矩阵逆左乘A矩阵。
如果A和B是两个标量数据，结论是都没有区别。
3，写出完成下列操作的命令。
（1）删除矩阵A的第7号元素A(7)=[]（2）将向量t的0元素用机器零来代替t(find(t==0))=eps（3）将含有12个元素的向量x转换成3*4矩阵newA=reshape(A,3,4)(4) 求一个字符串的ASCIIch=['a b c';'1 2 3'];abs ch(5) 产生和A同样大小的幺矩阵ones(size(A))(6) 从矩阵A提取主对角线元素，并以这些元素构成对角阵B。
A=[1 2 3;4 5 6;7 8 9];D=diag(A)B=diag([1,5,9])4, .要生产均值为3，方差为1的500个正态分布的随机序列，写出相应的表达式y=3+sqrt(1)*randn(500)5,(1)主对角元素1 1 5 9上三角矩阵1 -1 2 30 1 -4 20 0 5 20 0 0 9下三角矩阵1 0 0 05 1 0 03 0 5 011 15 0 9逆矩阵-0.1758 0.1641 0.2016 -0.0227-0.1055 -0.1016 -0.0391 0.0664-0.0508 -0.0859 0.1516 0.00230.3906 -0.0313 -0.1813 0.0281行列式的值1280秩4范数A1=norm(A,1)=20 A2=norm(A)= 21.3005 A3=norm(A,inf)=35条件数A1=cond(A,1)= 14.4531 A2=cond(A,2)= 11.1739 A3=cond(A,inf)= 22.0938迹16（2）主对角元素0.43 4上三角矩阵0.4300 43.0000 2.00000 4.0000 21.0000下三角矩阵0.4300 0 0-8.9000 4.0000 0逆矩阵0.0022 -0.01750.0234 -0.0017-0.0035 0.0405行列式的值秩2范数B1=norm(B,1)= 47 B2=norm(B)= 43.4271 B3=norm(B,inf)= 45.4300 条件数B2=cond(B,2)= 1.9354迹4.43006，all(A)=0any(A)=1isnan(A)= 0 1 0 0 0 0 0isinf(A)= 0 0 1 1 0 0 0isfinite(A)= 1 0 0 0 1 1 1。
第2章 MATLAB数据及其运算_习题答案
……………………………………………………………最新资料推荐…………………………………………………第2章 MATLAB数据及其运算习题2一、选择题1．下列可作为MA TLAB合法变量名的是（）。
DA．合计B．123 C．@h D．xyz_2a 2．下列数值数据表示中错误的是（）。
CA．+10 B．1.2e-5 C．2e D．2i3．使用语句t=0:7生成的是（）个元素的向量。
AA．8 B．7 C．6 D．54．执行语句A=[1,2,3;4,5,6]后，A(3)的值是（）。
BA．1 B．2 C．3 D．45．已知a为3×3矩阵，则a(:,end)是指（）。
DA．所有元素B．第一行元素C．第三行元素D．第三列元素6．已知a为3×3矩阵，则运行a (1)=[]后（）。
AA．a变成行向量B．a变为2行2列C．a变为3行2列D．a变为2行3列7．在命令行窗口输入下列命令后，x的值是（）。
B>> clear>> x=i*jA．不确定B．-1 C．1D．i*j8．fix(354/100)+mod(354,10)*10的值是（）。
DA．34 B．354 C．453D．439．下列语句中错误的是（）。
BA．x==y==3 B．x=y=3C．x=y==3 D．y=3,x=y10．find(1:2:20>15)的结果是（）。
CA．19 20 B．17 19C．9 10 D．8 911．输入字符串时，要用（）将字符括起来。
CA．[ ] B．{ } C．' ' D．" "12．已知s='显示"hello"'，则s的元素个数是（）。
AA．9 B．11 C．7 D．181 / 31 / 31 / 31……………………………………………………………最新资料推荐…………………………………………………13．eval('sqrt(4)+2')的值是（）。
MATLAB习题参考答案（胡良剑
第一章 MATLAB 入门4、求近似解解：>> x=-2:0.05:2;y=x.^4-2.^x两个近似解：y1=f(-0.85)= -0.0328; y2=f(1.250)= 0.0630第二章 MATLAB 编程与作图1、 设x 是数组，求均值和方差解：函数文件如下：function [xx,s]=func1(x)n=length(x);xx=sum(x)/n;s=sqrt((sum(x.^2)-n*xx^2)/(n-1));命令窗口：>> x=[1 2 3 4 5];[xx,s]=func1(x)2、求满足的最小m 值 100)1ln(0>+∑=m n n s=0;n=0;while(s<=100)s=s+log(1+n);n=n+1;endn,s3、用循环语句形成Fibonacci 数列,....4,3,,12121=+===−−k F F F F F k k k 。
并验证极限2511+→−k k F F （提示：计算至两边误差小于精度1e-8为止） 解： 求Fibonacci 数列的函数文件：function f=fun(n)if n<=2f=1;elsef=fun(n-1)+fun(n-2);end验证极限的函数文件：function [k,a]=funTest(e)a=abs(1-(1+sqrt(5))/2);k=2;while(a>e)k=k+1;a=abs(fun(k)/fun(k-1)-(1+sqrt(5))/2);end命令行：>> [k,a]=funTest(10^-8)k =21a =9.7719e-009或者M 文件如下：clear; F(1)=1;F(2)=1;k=2;x=0;e=1e-8; a=(1+sqrt(5))/2;while abs(x-a)>ek=k+1; F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1);enda,x,k4、分别用for 和while 循环结构编写程序，求出∑==610123i i K ，并考虑一种避免循环语句的程序设计，比较各种算法的运行时间。
matlab答案第二章
第二章1．设A和B是两个同维同大小的矩阵，问：（1）A*B和 A.*B的值是否相等？（2）A./B和B.\A的值是否相等？（3）A/B和B\A的值是否相等？（4）A/B和B\A所代表的数学含义是什么？答：（1）不等（2）相等（3）不等（4）A/B含义： A乘B的逆，即A* inv(B)B\A含义： B的逆乘A，即inv(B)*A2.写出完成下列操作的命令（1）将矩阵A第2~5行中第1，3，5列元素赋给矩阵B。
（2）删除矩阵A的第7号元素。
（3）将矩阵A的每个元素值加30。
（4）求矩阵A的大小和维数。
（5）将向量t的0元素用机器0来代替。
（6）将含有12个元素的向量x转换成3x4矩阵。
（7）求一个字符串的ASCII码。
（8）求一个ASCII码所对应的字符。
答：设矩阵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=A(2:5,[1 3 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）设向量t=[1,2,0,5,6,0,9];p=find(t==0);t(p)=eps;或者：t(find(t==0))=eps;（6）reshape(x,3,4)（7）abs(x) 或 double(x)（8）char(x)3.下列命令执行后，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 4.已知：⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=14.35454.9632053256545410778.01023A 完成下列操作：（1） 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角3x2子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
matlab 第二章部分答案
1、答：6+7i中的i为虚数单位，是一个复数常量；而6+7*i中的i为一个表达式，是一个运算量。
i是虚数单位，等于sqrt(-1)；I在MATLAB里没有定义，可根据自己的需要定义为各种变量。
2、答：A*B是矩阵和矩阵相乘，A.*B是指它们的对应元素进行乘法运算，要求两矩阵维参数相同。
A./B表示A矩阵除以B矩阵的对应元素，如果A、B两矩阵具有相同的维数，B.\A 等价于A./B。
如果B矩阵是非奇异方阵，A/B等效于B矩阵的逆右乘A矩阵，也就是A*inv(A)；B\A等效于B的逆左乘A矩阵，也就是inv(B)*A。
如果A和B是两个标量数据，则是它们的对应元素进行相关运算。
3、答：（1）A(7)=[ ];（2）find(t==0)t(ans)=eps（3）a=[0:11]b=reshape(a,3,4)（4）abs('a')（5）A=[0:11]ones[A]（6）A=[1,2,3;4,5,6];diag(A)B=diag([1,5])4、答：a=3+sqrt(1)*randn(50,10)5、答：（1）>> A=[1,-1,2,3;5,1,-4,2;3,0,5,2;11,15,0,9];>> diag(A)ans =1159>> triu(A)ans =1 -12 30 1 -4 20 0 5 20 0 0 9>> tril(A)ans =1 0 0 05 1 0 03 0 5 011 15 0 9>> inv(A)ans =-0.1758 0.1641 0.2016 -0.0227 -0.1055 -0.1016 -0.0391 0.0664 -0.0508 -0.0859 0.1516 0.00230.3906 -0.0313 -0.1813 0.0281 >> det(A)ans =1280>> rank(A)ans =4>> norm(A,1)ans =20>> norm(A)ans =21.3005>> norm(A,inf)ans =35>> cond(A,1)ans =14.4531>> cond(A)ans =11.1739>> cond(A,inf)ans =22.0938>> trace(A)ans =16行列式的值：1280秩：4范数：norm(a,1)=20, norm(a)=21.3005, norm(a,inf)=35条件数：cond(a,1)=14.4531,cond(a)=11.1739,cond(a,inf)=22.0938 迹：16（2）>> B=[0.43,43,2;-8.9,4,21];>> diag(B)ans =0.43004.0000>> triu(B)ans =0.4300 43.0000 2.00000 4.0000 21.0000>> tril(B)ans =0.4300 0 0-8.9000 4.0000 0B有逆矩阵的充要条件是|B|≠0，因为B的行列式等于0，所以B没有逆矩阵。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第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。
（2）然后根据此结果，求-=0t dtdy 和2π=t dtdy 。
〖目的〗diff, limit 指令的应用。
如何理解运行结果。
〖解答〗syms td=diff(y) %求dy/dtd0_=limit(d,t,0,'left') %求dy/dt|t=0- dpi_2=limit(d,t,pi/2) %求dy/dt|t=pi/2 y =abs(sin(t)) d =sign(sin(t))*cos(t) d0_ = -1dpi_2 = 09求出dx x e x sin 7.110⎰--ππ的具有64位有效数字的积分值。
〖目的〗符号积分的解析解和符号数值解。
符号计算和数值计算的相互校验。
〖解答〗（1）符号积分syms x clear syms xy=exp(-abs(x))*abs(sin(x)) si=vpa(int(y,-10*pi,*pi),64) y =abs(sin(x))/exp(abs(x)) si =（2）数值计算复验xx=-10*pi:pi/100:*pi;sn=trapz(exp(-abs(xx)).*abs(sin(xx)))*pi/100 sn =10计算二重积分⎰⎰+211222)(x dydx y x 。
〖目的〗变上限二重积分的符号计算法。
〖解答〗syms x y f=x^2+y^2;r=int(int(f,y,1,x^2),x,1,2) r =1006/10511在]2,0[π区间，画出dt ttx y x⎰=sin )(曲线，并计算)5.4(y 。
在符号计算中，经常遇到计算结果是特殊经典函数的情况。
如何应用subs 获得超过16位有效数字的符号数值结果。
初步尝试ezplot 指令的简便。
〖解答〗（1）符号计算syms t x; f=sin(t)/t; y=int(f,t,0,x) % 将得到一个特殊经典函数y5=subs(y,x,sym('')) ezplot(y,[0,2*pi]) y =sinint(x) y5 =12345600.20.40.60.811.21.41.61.8xsinint(x)（2）数值计算复验tt=0::; tt(1)=eps;yn=trapz(sin(tt)./tt)* yn =12在0>n 的限制下，求xdx n y n ⎰=20sin )(π的一般积分表达式，并计算)31(y 的32位有效数字表达。
〖目的〗一般符号解与高精度符号数值解。
syms xsyms n positive f=sin(x)^n;yn=int(f,x,0,pi/2)y3s=vpa(subs(yn,n,sym('1/3'))) y3d=vpa(subs(yn,n,1/3)) yn =beta(1/2, n/2 + 1/2)/2 y3s =y3d =13有序列k a k x =)(，k b k h =)(，（在此0≥k ，b a ≠），求这两个序列的卷积∑=-=kn n k x n h k y 0)()()(。
〖目的〗符号离散卷积直接法和变换法。
〖解答〗 （1）直接法syms a b k n x=a^k; h=b^k;w=symsum(subs(h,k,n)*subs(x,k,k-n),n,0,k) %据定义 y1=simple(w) w =piecewise([a = b, b^k + b^k*k], [a <> b, (a*a^k - b*b^k)/(a - b)]) y1 =piecewise([a = b, b^k + b^k*k], [a <> b, (a*a^k - b*b^k)/(a - b)])（2）变换法（复验）syms zX=ztrans(a^k,k,z); H=ztrans(b^k,k,z); y2=iztrans(H*X,z,k) %通过Z 变换及反变换y2 =piecewise([b <> 0, (a*a^k)/(a - b) - (b*b^k)/(a - b)])〖说明〗符号计算不同途径产生的结果在形式上有可能不同，而且往往无法依靠符号计算本身的指令是它们一致。
此时，必须通过手工解决。
15求0,)(>=-ααt Ae t f 的Fourier 变换。
〖目的〗符号变量限定性定义的作用。
fourier 指令的应用。
〖解答〗syms A t wa=sym('a','positive');f=A*exp(-a*abs(t)); y=fourier(f,t,w) F=simple(y) y =(2*A*a)/(a^2 + w^2) F =(2*A*a)/(a^2 + w^2)17求4633)(23++++=s s s s s F 的Laplace 反变换。
〖解答〗syms s tF=(s+3)/(s^3+3*s^2+6*s+4); f=simple(ilaplace(F,s,t)) f =(3^(1/2)*sin(3^(1/2)*t) - 2*cos(3^(1/2)*t) + 2)/(3*exp(t))19求Tk ke k f )(λ-=的Z 变换表达式。
〖目的〗注意：变换中，被变换变量的约定。
〖解答〗syms lambda k T z;f_k=k*exp(-lambda*k*T); F_z=simple(ztrans(f_k,k,z)) F_z =(z*exp(T*lambda))/(z*exp(T*lambda) - 1)^220求方程2,122==+xy y x 的解。
〖目的〗solve 指令中，被解方程的正确书写，输出量的正确次序。
〖解答〗eq1='x^2+y^2=1'; eq2='x*y=2';[x,y]=solve(eq1,eq2,'x','y') x =(1/2 + (15^(1/2)*i)/2)^(1/2)/2 - (1/2 + (15^(1/2)*i)/2)^(3/2)/2 - (1/2 + (15^(1/2)*i)/2)^(1/2)/2 + (1/2 + (15^(1/2)*i)/2)^(3/2)/2 (1/2 - (15^(1/2)*i)/2)^(1/2)/2 - (1/2 - (15^(1/2)*i)/2)^(3/2)/2 - (1/2 - (15^(1/2)*i)/2)^(1/2)/2 + (1/2 - (15^(1/2)*i)/2)^(3/2)/2 y =(1/2 + (15^(1/2)*i)/2)^(1/2) -(1/2 + (15^(1/2)*i)/2)^(1/2) (1/2 - (15^(1/2)*i)/2)^(1/2) -(1/2 - (15^(1/2)*i)/2)^(1/2)23求微分方程045=+'x y y 的通解，并绘制任意常数为1时解的图形。
〖目的〗理解指令dsolve 的正确使用。
对dsolve 输出结果的正确理解。
ezplot 指令绘图时，如何进行线色控制。
如何覆盖那些不能反映图形窗内容的图名。
〖解答〗 （1）求通解reset(symengine) clear syms y xy=dsolve('*y*Dy+*x=0','x') y =2^(1/2)*(C3 - (5*x^2)/8)^(1/2) -2^(1/2)*(C3 - (5*x^2)/8)^(1/2)（2）根据所得通解中不定常数的符号写出“对其进行数值替代的指令”yy=subs(y,'C3',1) %将通解中的C3用1代替 yy =2^(1/2)*(1 - (5*x^2)/8)^(1/2) -2^(1/2)*(1 - (5*x^2)/8)^(1/2)（3）观察通解中两个分解的平方是否相同 yy(1)^2==yy(2)^2ans = 1（4）于是可考虑函数的平方关系syms Yfxy=Y^2-yy(1)^2 fxy =Y^2 + (5*x^2)/4 - 2（5）根据平方关系式画完整曲线clfezplot(fxy,[-2,2,-2,2]) axis square grid on（6）假如直接用“分解”画曲线，那么将是不完整的ezplot(yy(1)),hold oncc=get(gca,'Children');set(cc,'Color','r')ezplot(yy(2)),axis([-2 2 -2 2])legend('y(1)','y(2)'),hold off;title(' ') %覆盖不完全的图名gridaxis square-2-1.5-1-0.500.51 1.52-2-1.5-1-0.500.511.52xy(1)y(2)24求一阶微分方程2)0(,2=+=x bt at x的解。