请教怎么样用数值方法画Duffing方程的幅频特性曲线

如枫

请教怎么样用数值方法画Duffing方程的幅频特性曲线？
新建 Microsoft Word 文档.doc (15.5 KB, 下载次数: 13)

2011-4-14 09:14 上传

点击文件名下载附件

hsfy919

回复 1 # 如枫 的帖子

先用数值方法解出振幅x，程序论坛里有，再绘制出频率和振幅的曲线。

如枫

利用傅里叶变换求振幅这一部分程序不太清楚，
在论坛里面找到一个程序，有几个地方不明白
ff=90; %(Hz) 1/ff 表示周期
N=1024*16;
n=4; % how many period

grid=[0:N-1]';
incre= 1/ff/(N-1);  %将周期分成N-1等分，incre表示步长
t = incre*grid*n;

% plot the original signal
y=1.75*cos(2*pi*ff*t) + 2.3*sin(4*pi*ff*t)+4.3*sin(5*pi*ff*t);
figure(1);
plot(t,y);
title('Signal');
xlabel('time');
ylabel('signal');

% FFT
Y = fft(y)/N*2;
Pyy = abs(Y);

% plot the FFT
f = 1/incre*grid/n/(length(Pyy)-1); %求各阶谐波幅值？
figure(2);
plot(f(1:N/2),Pyy(1:N/2),'r')
title('Frequency content of y');
xlabel('frequency (Hz)');
ylabel('magnitude (Hz)');

1. Y = fft(y)/N*2; 这里为什么要除以N*2
2.   f = 1/incre*grid/n/(length(Pyy)-1);  横坐标频率为什么要做这样的转换？

请教高手！！！

tudexin

如枫 发表于 2011-4-15 21:45

                               
登录/注册后可看大图

利用傅里叶变换求振幅这一部分程序不太清楚，
在论坛里面找到一个程序，有几个地方不明白
ff=90; %(Hz) 1 ...

你这是幅值谱图程序吧?不是幅频特性图程序好像，呵呵 个人观点，不知道对错

octopussheng

FFT后得到的是某一组解的幅值-频率特性，和你题目中说 的不一样吧

如枫

上面的程序是和标题不一样，继续请教
怎么样用数值方法从原方程得到频幅特性曲线？

伤痕累累

先用傅里叶转化成级数叠加（类似y的多个简谐函数叠加），得到的结果将原程序的y替换一下就可以得到了吧，个人意见

伤痕累累

今天发现新的方法，发现我上面说的不妥。
将duffing方程数值解之后于将对应的数值快速值傅里叶转化。我给个例子
w=3000*pi;  
T=2*pi; ts=[0:T/512:128*T];
y=z(24577:end,1);  %取出的数值解
fs=4096; N=length(y);
zz=fft(y); zz=fftshift(zz);
f=(0:N-1)*fs/N;
plot(f,zz)