矩阵实验室实验

1、第三章3.5（1） n=100;i=1:1:n;y=sum(1./(i.2);Y第二种方法y=0;n=100;for i=1:1:n y=y+1/(i2);endy3.5（2）y=1;n=100;fori=1:100 y=y*(2*i)*(2*i)/(2*i-1)*(2*i+1);endy第四章4.1y=x-x.3/6;subplot(1,2,1)plot(x,y)subplot(1,2,2)ezplot('x2+2*y2-64')4.2t=-pi:0.1*pi:pi;y=1./(1+exp(-t);subplot(2,2,1);bar(t,y,'b');tit

2、le('条形图');subplot(2,2,2);stairs(t,y,'r');title('阶梯图');subplot(2,2,3);stem(t,y,'y');title('杆图');subplot(2,2,4);semilogy(t,y,'g');title('对数图');t=0:pi/10:2*pi;r=5*cos(t)+4;4.3subplot(1,2,1)polar(t,r,'.')subplot(1,2,2)t1=-pi/3:pi/250:pi/3;r

3、1=5*(sin(t1).2)./cos(t1);polar(t1,r1,'r')4.4（1）t=0:0.1*pi:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t);z=t;plot3(x,y,z)4.4（2）z=5*ones(size(x);x=-5:5;y=-5:5;x,y=meshgrid(x,y);mesh(x,y,z)shading interp4.5maxx=200;m=moviein(maxx);x,y,z=sphere(40);for r=1:maxx X=x*r; Y=y*r; Z=z*r; mesh(X,Y,Z);

4、axis(-maxx,maxx,-maxx,maxx,-maxx,maxx) axis square; m(:,r)=getframe;endmovie(m,2)5.2(2)x=1:1:100;y=sqrt(x);interp1(x,y,x,'cubic')第五章5.2(1)x1=0:pi/12:pi/2;y1=0,0.2588,0.5000,0.7071,0.8660,0.9659,1.000;p1=polyfit(x1,y1,5)t1i=0:pi/12:pi/2;y1i=polyval(p1,t1i)x=0:pi/12:pi/2;y=sin(x);interp1(x,y,x

5、,'linear')5.3N=64;T=5;t=linspace(0,T,N);h=exp(-t);dt=t(2)-t(1);f=1/dt;H=fft(h);F=H(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')xlabel('Frequency');ylabel('|F(k)|')5.4P=2,-3,0,5,13;Q=1,5,8;p1=polyder(P)p2=polyder(P,Q)p,q=polyder(P,Q)5.5P1=1,2,4,0,5;P2=0,0,0,1,2;P3=0,0,1,

6、2,3;P=P1+P2+P3x=roots(P)x=-1,1.2,-1.4;0.75,2,3.5;0,5,2.5;y1=polyval(P,x)y2=polyvalm(P,x) 第八章8.2（2）sym x;limit(sqrt(pi)-sqrt(acos(x)/sqrt(x+1),x,-1,'right')8.5syms x n;s=symsum(x(2*n-1)/(2*n-1),'x',inf)8.8syms x y n;x y=dsolve('x*D2y+(1-n)*Dy+y=0','y=(0)=0','Dy(0)=

7、0','x')第十章计算器（编程）hf=figure('Color',1,0,1,'Name','计算器','Menubar','none');uicontrol(hf,'Style','Text','Units','normalized',. 'Position',0.05,0.8,0.45,0.1,'String','输入框',. 'Horizontal',&#

8、39;center','Back',1,0.1,0.6);uicontrol(hf,'Style','Text','Units','normalized',. 'Position',0.5,0.8,0.45,0.1,'String','输出框',. 'Horizontal','center');h1=uicontrol(hf,'Style','Edit','Units','

9、normalized',. 'Position',0.05,0.6,0.25,0.1,. 'Horizontal','center','Back',1,0.2,0.1);h2=uicontrol(hf,'Style','Edit','Units','normalized',. 'Position',0.05,0.4,0.25,0.1,. 'Horizontal','center','Back',1

10、,0.2,0.1);h3=uicontrol(hf,'Style','Text','Units','normalized',. 'Position',0.5,0.4,0.25,0.25,. 'Horizontal','center','Back',1,0.2,0.1);COMM1='a=str2num(get(h1,''String'');',. 'b=str2num(get(h2,''Strin

11、g'');',. 's=a+b;','set(h3,''String'',s);'COMM2='a=str2num(get(h1,''String'');',. 'b=str2num(get(h2,''String'');',. 's=a-b;','set(h3,''String'',s);'COMM3='a=str2num(get(h1,&

12、#39;'String'');',. 'b=str2num(get(h2,''String'');',. 's=a*b;','set(h3,''String'',s);'COMM4='a=str2num(get(h1,''String'');',. 'b=str2num(get(h2,''String'');',. 's=a/b;','

13、;set(h3,''String'',s);'uicontrol(hf,'Style','Push','Units','normalized',. 'Position',0.10,0.1,0.1,0.12,. 'String','+','Back',1,0.2,0.1,. 'Call',COMM1);uicontrol(hf,'Style','Push','Units

14、9;,'normalized',. 'Position',0.20,0.1,0.1,0.12,. 'String','-','Back',1,0.2,0.1,. 'Call',COMM2);uicontrol(hf,'Style','Push','Units','normalized',. 'Position',0.30,0.1,0.1,0.12,. 'String','*','B

15、ack',1,0.2,0.1,. 'Call',COMM3);uicontrol(hf,'Style','Push','Units','normalized',. 'Position',0.40,0.1,0.1,0.12,. 'String','/','Back',1,0.2,0.1,. 'Call',COMM4);计算器（图形界面）function pushbutton1_Callback(hObject, eventdata

16、, handles)a=str2num(get(handles.edit1,'string');b=str2num(get(handles.edit4,'string');s=a+b;set(handles.edit2,'string',s);10.1screen=get(0,'ScreenSize');figure('color',1,1,1,'MenuBar','none')hm=uimenu(gcf,'Label','File');hm1=uim

17、enu(hm,'Label','New');uimenu(hm1,'Label','New Item');hm2=uimenu(hm,'Label','Open');uimenu(hm2,'Label','Open Item');hm3=uimenu(hm,'Label','Save','Enable','off');uimenu(hm3,'Label','Save Item

18、9;);hm4=uimenu(hm,'Label','Exit','Call','close(gcf)');hn=uimenu(gcf,'Label','help');hn1=uimenu(hn,'Label','About.');uimenu(hn1,'Label','Help Item','Call',. 'set(hm3,''Enable'',''on'

19、;')');10.2x=0:1:10;y=x.*x;h1=plot(x,y);hc=uicontextmenu;h1c=uimenu(hc,'Label','颜色');uimenu(h1c,'Label','Red','Call',. 'set(h1,''Color'',''r'')');uimenu(h1c,'Label','Blue','Call',. 'set

20、(h1,''color'',''b'')');uimenu(h1c,'Label','Yellow','Call',. 'set(h1,''color'',''y'')');set(h1,'UIContextMenu',hc);10.3hpopup=uicontrol(gcf,'Style','popup','String',. &#

21、39;cos|sum|fix|exp|roots','CallBack',. 'cbcol=''C'',''S'',''F'',''E'',''R''',. 'set(gcf,''Function'',cbcol(get(hpopup,''Value'')');hedit=uicontrol(gcf,'Styl

22、e','edit',. 'Unit','normalized','String','输入多项式系数',. 'Position' ,0.6,0.2,0.2,0.6,'Max',10);hpush1=uicontrol(gcf,'Style','push','Unit','normalized',. 'Position',0.8,0.3,0.15,0.15,'String','求根',. 'Call','comm(hedit)');hpush2=uicontrol(gcf,'Style','push','Unit','normalized',. 'Position',0.45,0.45,0.1,0.1,'CallBack',. 'set(hpush2,''Position'',rand,rand,0.1,0.1)');10.4设置字体