24小时热门版块排行榜    

CyRhmU.jpeg
查看: 468  |  回复: 1

haitao_lee

金虫 (著名写手)

[求助] 线性系统做一个扫频模拟,不知道问题出在哪里?

我想试着对一个线性系统做一个扫频模拟, 下面代码中他的固有频率显然是1,但是我模拟出的峰值确实在0.5处。 如果按照稳态来算的频率在1处的振幅确实比0.5处要大。怎么来模拟瞬态扫频结果,还是一头雾水。。
clc
clear all
global w0 alpha
w0=0;
alpha=0.001
x0=[0;0];
TSPAN=0:0.01:1000;
[T,Y] = ode45('gongzhen',TSPAN,x0)
ww=w0+alpha.*TSPAN;
figure
plot(ww,Y(:,1))


function dx=gongzhen(t,x)
global w0  alpha
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=-0.01*x(2)-1*x(1)-1*sin((w0+alpha*t)*t);

图太小了传不上去,拜托了,各位老师
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haitao_lee

金虫 (著名写手)

我知道是哪里出错了。。。
2楼2015-09-16 08:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 haitao_lee 的主题更新
信息提示
请填处理意见