matlab教程2011a(张志涌)-课后答案

ezplot(subs(S(1),'C3',1),[-2,2 -2,2],1), hold on
ezplot(subs(S(2),'C3',1),[-2,2 -2,2],1) %用此两条指令绘圆，在y=0处有间隙
ezplot(subs(y^2-(S(1))^2, 'C3', 1),[-2,2 -2,2],2) %用椭圆方程绘图不产生间隙
L=length(t)
for k=1:L
yy(k)=1-exp(-0.5*t(k))*cos(2*t(k))
end
figure(1),plot(t,yy)
%(2)
y=1-exp(-0.5.*t).*cos(2.*t),
figure(2),plot(t,y)
figure(3),subplot(1,2,1),plot(t,y)
ss=d*(cumtrapz(y)) %计算梯形法累计积分并绘积分曲线
plot(t,y,t,ss,'r'),hold on
%用find指令计算y(4.5),并绘出该点
y4_5=ss(find(t==4.5))
%插值法计算y(4.5),并绘出该点
yi=interp1(t,ss,4.5), plot(4.5,yi,'r+')
subplot(1,2,2),plot(t,yy)
%习题3-6
clear,format long,
rand('twister',1),A=rand(3,3),
B=diag(diag(A)),
C=A.*(~B)
%或C=A-B%或C=triu(A,1)+tril(A,-1)
%习题3-7
clear,
x=-3*pi:pi/15:3*pi; y=x;
N=20
diff_y1=(diff(y(1:N:end)))./diff(t(1:N:end));
gradient_y1=(gradient(y(1:N:end)))./gradient(t(1:N:end));
% plot(t(1:end-1),diff_y,t,gradient_y)
t1=t(1:N:end);
%习题2-9
clear all, syms x;
f=exp(-abs(x)).*abs(sin(x));
fint=int(f,x,-5*pi,1.7*pi), digits(64), vpa(fint)
class(fint)
%习题2-10
clear all,syms x y,f=x.^2+y.^2,
fint=(int(int(f,y,1,x.^2),x,1,2)), double(fint)
%习题2-11
clear all, syms t x; f=sin(t)/t, yx=int(f,t,0,x), ezplot(yx,[0 2*pi])
fint=subs(yx,x,4.5),%或yxd=int(f,t,0,4.5),fint=double(yxd)
[I1,J1]=find(A>0.5)
subindex_A=sub2ind(size(A),I1,J1)
subindex_A=find(A>0.5)
[I,J]=ind2sub(size(A),subindex_A)
index_A=[I J]
%122_2
rand('twister',0),A=rand(3,5),B=A>0.5,C=find(B)
% (2)绘图
surf(X,Y,Z); shading interp
% (3)“无裂缝”图形的全部指令：
x=-3*pi:pi/15:3*pi;y=x';
[X,Y]=meshgrid(x,y); % X=ones(size(y))*x;Y=y*ones(size(x));
X=X+(X==0)*eps;
G =
piecewise([1/2 < y, y^2/4], [y <= 1/2, y])
%习题2-1
class(vpa(sym(3/7+0.1),4))
%习题2-2
findsym(sym('sin(w*t)'),1)
findsym(sym('a*exp(-X)'),1)
findsym(sym('z*exp(j*th)'),1)
hold on, plot(4.5,fint,'*r')
%习题2-12
% clear all, syms x n; f=(sin(x))^n; yn=int(f,'x',0,pi/2), class(yn)
clear all, syms x n; syms n positive ; f=(sin(x))^n; yn=int(f,'x',0,pi/2), class(yn)
[ci,cj]=ind2sub(size(A),C) %此法太繁
[Ci,Cj]=find(B)
%122_5
t=linspace(1,10,100)
(1) y=1-exp(-0.5.*t).*cos(2.*t),plot(t,y)
(2) L=length(t)
for k=1:L, yy(k)=1-exp(-0.5*t(k))*cos(2*t(k)), end, plot(t,yy)
%122_6
clear,format long, rand('twister',1),A=rand(3,3),
B=diag(diag(A)),C=A.*(~B)%或C=A-B%或C=triu(A,1)+tril(A,-1)
%习题3-3
% s=sign(randint(1,1000,[],123)-.5);
[X,Y]=meshgrid(x,y); % X=ones(size(y))*x;Y=y*ones(size(x));
warning off;
Z=sin(X).*sin(Y)./X./Y;
% (1)“非数”数目
m=sum(sum(isnan(Z))); %k=Z(isnan(Z));m=length(k)
colormap([0 0 1]) %用ezplot(fun)绘图时，如果fun中只有一个参数，绘图的颜色是蓝色；如果fun中有两个参数，绘图的颜色是绿色，此指令设置图形颜色为蓝。
grid on
S1=subs(S(1),'C3',1)
subs(S1,{x},{1.6^(1/2)})
y=double(solve(S1))
A=reshape(1:10*k,k,10);
if k==1
Sa(k,:)=A;
else
Sa(k,:)=sum(A);
end
end
Sa
第四章
%习题4-1
clear all
load prob_401;
diff_y=diff(y)./diff(t);
gradient_y=gradient(y)./gradient(t);
ezplot(subs(S(1),'C3',1),[-2,2 -2,2],1), hold on
ezplot(subs(S(2),'C3',1),[-2,2 -2,2],1)
%解为S =
% 2^(1/2)*(C3 - (5*x^2)/8)^(1/2)
% -2^(1/2)*(C3 - (5*x^2)/8)^(1/2)
length(t1)
figure(2)
subplot(1,2,1),plot(t,y,t1(1:end-1),diff_y1)
subplot(1,2,2),plot(t,y,t1,gradient_y1)
%习题4-2
d=0.5; tt=0:d:10; t=tt+(tt==0)*eps; y=sin(t)./t; s=d*trapz(y) %计算出积分值
evalin(symengine,'anames(Properties)')
evalin(symengine,'getprop(x)')
%习题2-4
7/3, pi/3 , pi*3^(1/3)
a = pi*3^(1/3) ,
b = sym( pi*3^(1/3) ),
c = sym( pi*3^(1/3),'d'),
t=linspace(y(2),y(1),100)
S2=subs(S(2),'C3',1)
figure
plot(t,subs(S1,x,t)), hold on
plot(t,subs(S2,x,t))
axis([-2,2 -2,2])
%习题2-24
x=dsolve('Dx=a*t^2+b*t','x(0)=2','t')
fk =2/(2*k+1)*(((x-1)/(x+1))^(2*k+1))
s=symsum(fk,k,0,inf)
s1=simple(s)
%习题2-8
clear all, syms t
y=abs(sin(t))
df=diff(y),class(df)
df1=limit(df,t,0,'left')
df2=subs(df,'t',sym(pi/2))
%习题2-25
[f,g]=dsolve('Df=3*f+4*g','Dg=-4*f+3*g','f(0)=0','g(0)=1','x')