matlab课后习题答案 附图.docx
- 文档编号:393091
- 上传时间:2022-10-09
- 格式:DOCX
- 页数:25
- 大小:265.79KB
matlab课后习题答案 附图.docx
《matlab课后习题答案 附图.docx》由会员分享,可在线阅读,更多相关《matlab课后习题答案 附图.docx(25页珍藏版)》请在冰豆网上搜索。
matlab课后习题答案附图
习题2.1
画出下列常见曲线的图形
(1)立方抛物线
命令:
symsxy;
ezplot('x.^(1/3)')
(2)高斯曲线y=e^(-X^2);
命令:
clear
symsxy;
ezplot('exp(-x*x)')
(3)笛卡尔曲线
命令:
>>clear
>>symsxy;
>>a=1;
>>ezplot(x^3+y^3-3*a*x*y)
(4)蔓叶线
命令:
>>clear
>>symsxy;
>>a=1
ezplot(y^2-(x^3)/(a-x))
(5)摆线:
命令:
>>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)对数螺线
命令:
clear
theta=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.2
1.求出下列极限值
(1)
命令:
>>symsn
>>limit((n^3+3^n)^(1/n))
ans=
3
(2)
命令:
>>symsn
>>limit((n+2)^(1/2)-2*(n+1)^(1/2)+n^(1/2),n,inf)
ans=
0
(3)
命令:
symsx;
>>limit(x*cot(2*x),x,0)
ans=
1/2
(4)
命令:
symsxm;
limit((cos(m/x))^x,x,inf)
ans=
1
(5)
命令:
symsx
>>limit(1/x-1/(exp(x)-1),x,1)
ans=
(exp
(1)-2)/(exp
(1)-1)
(6)
命令:
symsx
>>limit((x^2+x)^(1/2)-x,x,inf)
ans=
1/2
练习2.4
1.求下列不定积分,并用diff验证:
(1)
>>Clear
>>symsxy
>>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)
clear
symsxy
y=1/(1+exp(x));
f=int(y,x)
f=
-log(1+exp(x))+log(exp(x))
symsxy
y=-log(1+exp(x))+log(exp(x));
yx=diff(y,x);
y1=simple(yx)
y1=
1/(1+exp(x))
(3)
symsxy
y=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^2
clear
symsxyy=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)
symsxy
y=sec(x)^3;
f=int(y,x)
f=
1/2/cos(x)^2*sin(x)+1/2*log(sec(x)+tan(x))
clear
symsxy
y=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)^3
2.求下列积分的数值解
1)
clear
symsx
y=int(x^(-x),x,0,1)
y=
int(x^(-x),x=0..1)
vpa(y,10)
ans=
1.291285997
2)
clear
symsx
y=int(exp(2*x)*cos(x)^3,x,clear
symsx
y=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)
>>clear
>>symsx
>>y=int(1/(2*pi)^(1/2)*exp(-x^2/2),0,1);
>>vpa(y,14)
ans=
.34134474606855
2(4)
>>clear
>>symsx
>>y=int(x*log(x^4)*asin(1/x^2),1,3);
Warning:
Explicitintegralcouldnotbefound.
>Insym.intat58
>>vpa(y,14)
ans=
2.4597721282375
2(5)
>>clear
>>symsx
>>y=int(1/(2*pi)^(1/2)*exp(-x^2/2),-inf,inf);
>>vpa(y,14)
ans=
.99999999999999
练习2.5
1判断下列级数的收敛性,若收敛,求出其收敛值。
1)symsn
s1=symsum(1/n^(2^n),n,1,inf)
s1=
sum(1/(n^(2^n)),n=1..Inf)
vpa(s1,10)
ans=
1.062652416
因此不收敛
2)symsn
s1=symsum(sin(1/n),n,1,inf)
s1=
sum(sin(1/n),n=1..Inf)
vpa(s1,10)
ans=
sum(sin(1/n),n=1..Inf)
不收敛
(3)
>>clear
>>symsn
>>s=symsum(log(n)/n^3,n,1,inf)
s=
-zeta(1,3)
收敛
(4)symsn
s1=symsum(1/(log10(n))^n,n,3,inf)
s1=
sum(1/((log(n)/log(10))^n),n=3..inf)
不收敛
(5)symsn
s1=symsum(1/n*log10(n),n,2,inf)
s1=
sum(1/n*log(n)/log(10),n=2..Inf)
不收敛
(6)
>>clear
>>symsn
>>s=symsum((-1)^n*n/n^2+1,n,1,inf)
s=
sum((-1)^n/n+1,n=1..Inf)
不收敛
习题3.1
1)clear;
[x,y]=meshgrid(-30:
0.3:
30);
z=10*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2);
>>meshc(x,y,z)
clear
>>[x,y]=meshgrid(-30:
0.1:
30);
>>z=10*sin((x^2+y^2)^(1/2))/(1+x^2+y^2)^(1/2)
mesh(x,y,z)
1.
2.取适当的参数绘制下列曲面的图形。
(1)
clear
>>a=-2:
0.1:
2;
>>b=-3:
0.1:
3;
>>[x,y]=meshgrid(a,b);
>>z=(1-(x.^2)/4-(y.^2)/9).^(1/2);
>>mesh(x,y,z)
>>holdon
mesh(x,y,-z)
(2)
clear
>>a=-1:
0.1:
1;
>>b=-2:
0.1:
2;
[x,y]=meshgrid(a,b);
>>z=(4/9)*(x.^2)+(y.^2);
>>mesh(x,y,z)
(4)
clear
>>[x,y]=meshgrid(-1:
0.1:
1);
>>z=(1/3)*(x.^2)-(1/3)*(y.^2);
>>mesh(x,y,z)
习题3.2
P49/例3.2.1
命令:
symsxy
limit(limit((x^2+y^2)/(sin(x)+cos(y)),0),pi),
ans=
-pi^2
limit(limit((1-cos(x^2+y^2))/((x^2+y^2)),0),0),
ans=
0
P49/例3.2.2
命令:
clear;symsxyzdxdydzzxzzyzxxzxy
z=atan(x^2*y)
z=
atan(x^2*y)
zx=diff(z,x),zy=diff(z,y)
zx
2*x*y/(1+x^4*y^2)
zy=
x^2/(1+x^4*y^2)
dz=zx*dx+zy*dy,
dz=
2*x*y/(1+x^4*y^2)*dx+x^2/(1+x^4*y^2)*d
zxx=diff(zx,x),zxy=diff(zx,y)
zxx=
2*y/(1+x^4*y^2)-8*x^4*y^3/(1+x^4*y^2)^2
zxy=
2*x/(1+x^4*y^2)-4*x^5*y^2/(1+x^4*y^2)^2
3.2.1作图表示函数z=x*exp(-x^2-y^2)(-1 clear >>a=-1: 0.1: 1; >>b=0: 0.1: 2; >>[x,y]=meshgrid(a,b); >>z=x.*exp(-x.^2-y.^2); >>[px,py]=gradient(z,0.1,0.1); contour(a,b,z),holdon, >>quiver(a,b,px,py),holdoff 习题3.4 1.解下列微分方程 (1)y=dsolve('Dy=x+y','y(0)=1','x') y= -x-1+2*exp(x) x=[123] x=123 -x-1+2*exp(x) ans= 3.436611.778136.1711 (2)x'=2*x+3*y,y'=2*x+y,x(0)=-2,y(0)=2.8,0 新建M函数 functiondy=weifen1(t,y) dy=zeros(2,1); dy (1)=2*y (1)+
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab课后习题答案 附图 matlab 课后 习题 答案