资源描述:
《答案(可通过左边的bookmark迅速查找答案).pdf》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、答案(可通过左边的bookmark迅速查找答案)第一次练习1.在主窗口输入:a=exp(1)A=piπeexp(pi)-pi^exp(1)%比较e和π的大小2.A=rand(3,3)%产生随机数矩阵B=inv(A)%求A的逆矩阵C=B*AD=B.*AE=sqrt(C)formatshort%设置数据的形式E第二次练习1.%%质点运动图形clear;t=0:0.05:1000;w1=2;w2=[1,2,4,3,5,7,sqrt(3),sqrt(5),exp(1)];%w2可以选不同的值x=sin(w1*t);forn=1:9;y=sin(w2(n)*t);s
2、ubplot(3,3,n);plot(x,y);end2.%%本征函数静态%画出四种本征函数系sin(nπx/l),cos(nπx/l),sin((n+1/2)πx/l),cos((n+1/2)πx/l)的图形,其中n取1,2,3,4。clear;clf;n=1:3;x=0:0.01:1;[XN]=meshgrid(x,n)Y1=sin(pi*N.*X);Y2=cos(pi*N.*X);Y3=sin(pi*(N+1/2).*X);Y4=cos(pi*(N+1/2).*X);Y=cat(3,Y1,Y2,Y3,Y4);fork=1:4subplot(2,2,k
3、)plot(x,Y(:,:,k),'b')axis([0,1,-1,1])end第三次练习1.%%本征函数动态clear;clf;clc;n=1:3;%x=0:0.04:1;t=0:0.02*pi:0.8*pi;[XNT]=meshgrid(x,n,t);fork=1:50subplot(2,1,1)Y1=sin(pi*N(:,:,1).*X(:,:,1)).*cos(pi*N(:,:,1).*T(:,:,k));plot(x,Y1);axis([0,1,-1,1])m1(:,k)=getframe;subplot(2,1,2)Y2=-cos(pi*N(:
4、,:,1).*X(:,:,1)).*sin(pi*N(:,:,1).*T(:,:,k));plot(x,Y2);axis([0,1,-1,1])m2(:,k)=getframe;endmovie(m1,1);movie(m2,1);2.%画复变函数图形1/2%(z-0.5)z=cplxgrid(50);cplxmap(z,sqrt(z-0.5))z%画e的图形z=cplxgrid(50);cplxmap(z,exp(z))%画lnz的图像z=cplxgrid(50);cplxmap(z,log(z))3.%两端固定弦的级数解%初速度为零,初位移不为零cle
5、ar;clc;functionjxjN=50t=0:0.005:2.0;x=0:0.001:1;ww=wfun(N,0);ymax=max(abs(ww));h=plot(x,ww,'linewidth',3);axis([0,1,-ymax,ymax])sy=[];forn=2:length(t)ww=wfun(N,t(n));set(h,'ydata',ww);drawnow;sy=[sy,sum(ww)];endfunctionwtx=wfun(N,t)x=0:0.001:1;a=1;wtx=0;forI=1:NifI~=7wtx=wtx+0.05*
6、((sin(pi*(7-I)*4/7)-sin(pi*(7-I)*3/7)).../(7-I)/pi-(sin(pi*(7+I)*4/7)-sin(pi*...(7+I)*3/7))/(7+I)/pi)*cos(I*pi*a*t).*sin(I*pi*x);elsewtx=wtx+0.05/7*cos(I*pi*a*t).*sin(I*pi*x);endend%初速度不为零,初位移为零functionpsiN=50;t=0:0.005:2.0;x=0:0.001:1;ww=psi1fun1(N,0);h=plot(x,ww,'linewidth',3);a
7、xis([0,1,-0.08,0.08])sy=[];forn=2:length(t)ww=psi1fun1(N,t(n));set(h,'ydata',ww);drawnow;sy=[sy,sum(ww)];endfunctionwtx=psi1fun1(N,t)x=0:0.001:1;a=1;wtx=0;fork=1:NBk=2/(k*k*pi*pi)*(cos(3*k*pi/7)-cos(4*k*pi/7));wtx=wtx+Bk*sin(k*pi*t)*sin(k*pi*x);end第四次练习1.%mandelbrot集的绘制dxy=0.001;N
8、1=2600;N2=2700;N=70;RGB=zeros(N1,