实变函数 MATLAB 程序

%实变函数中的阶段函数及其原函数的图像， %以y1=n*sin(x)为例，画出两个函数在0-10*pi之间的函数图像 n1=5;x00=5*pi; x0=0:0.01:x00; y0=n1*sin(x0); plot(x0,y0,'r.','markersize',4);hold on syms n;syms k; [x1,y1]=solve('y=n*sin(x)','y=k','x','y'); x1=subs(x1,n,n1);x1=subs(x1,k,3/4*n1);y1=subs(y1,k,3/4*n1); temp1=x1;temp2=y1; i=1;x(i,1)=x1;x(i,2)=y1; while x(i,1)<x00 x(i+1,1)=x(i,1)+2*pi; x(i+1,2)=y1; i=i+1; end [x1,y1]=solve('y=n*sin(x)','y-k=0','x','y'); x1=subs(x1,n,n1);x1=subs(x1,k,-3/4*n1);y1=subs(y1,k,-3/4*n1); temp3=x1;temp4=y1; x(i,1)=x1+2*pi;x(i,2)=y1; while x(i,1)<x00 x(i+1,1)=x(i,1)+2*pi; x(i+1,2)=y1; i=i+1; end x1=pi-temp1; x(i,1)=x1;x(i,2)=temp2; while x(i,1)<x00 x(i+1,1)=x(i,1)+2*pi; x(i+1,2)=temp2; i=i+1; end x1=pi-temp3; x(i,1)=x1;x(i,2)=temp4; while x(i,1)<x00 x(i+1,1)=x(i,1)+2*pi; x(i+1,2)=temp4; i=i+1; end x=x(1:(size(x,1)-1),:);