在命令窗输入如下内容就得到图形了: [t,x]=ode45( godhua ,[0,10],[0,0,0,1]) subplot(131) plot(t,x(:,4)) xlabel( t )
在命令窗输入如下内容就得到图形了: [t,x]=ode45( godhua ,[0,10],[0,0,0,1]) subplot(131) plot(t,x(:,4)) xlabel( t ) ylabel( y ) subplot(132) plot(t,x(:,1)) xlabel( ...
在命令窗输入如下内容就得到图形了: [t,x]=ode45( godhua ,[0,10],[0,0,0,1]) subplot(131) plot(t,x(:,4)) xlabel( t ) ylabel( y ) subplot(132) plot(t,x(:,1)) xlabel( ...
%radon transform clear all % N=800 n=1:N fs=200 t=n/fs x1=exp(j*2*pi*(5*t+0.5*5*t.^2)) x2=exp(j*2*pi*(5*t+0.5*15*t.^2)) x=x1+x2 ...
y3k=fft(u,(m+n-2)/4) i=1:(m+n-2)/4 subplot(5,2,9) stem(i,u) title( 滤波后上采样 ) k=1:(m+n-2)/4 subplot(5,2,10) stem(k,y3k) title( 上采样频谱 )...
%球体 close all; G=6.67e-11; R=2;%球体半径 p=4.0;%密度 D=10.0;%深度 M=(4/3)*pi*R^3*p;%质量 x=-20:1:20; g=G*M*D./((x.^2+D^2).^(3/2)); Vxz=-3*G*M*D.*x./((x.^2+D^2...
x=[1,2,0,-1,3,2];h=[1,-1,1]; y1=x*h(1); y2=x*h(2); y3=x*h(3); Y1=[0,0,y1]; Y2=[0,y2,0]; Y3=[y3,0,0]; y=Y1+Y2+Y3; L=-2:1:5; figure(1); subplot(211);ste...