在命令窗输入如下内容就得到图形了: [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 ) ylabel( v ) subplot(133) plot(t,x(:,2)) xlabel( t ) ylabel( \\theta )

分析双极性不归零码的功率谱特性请输入随机序列的长度N= -ubplot(312) plot(f,lpf) title( 低通滤波器 ) axis([-fs/2 fs/2-0.2 1.2]) subplot(313) plot(f,abs(yr_lpf)/N) title( 经低通滤波后的频谱 ) xlabel( f /Hz )