利用随机信号进行离散系统分析的问题
连续系统微分方程为：y''(t)+2ζΩ y'(t)+Ω^2*y(t)=x(t)现需要将系统离散化，并加随机噪声输入，用Y(w)/X(w)的方法求频响函数H(w)，不知道程序对不对，这块儿概念比较模糊，请各位前辈指点一下：
f=50; %系统频率
ts=0.01; %采样周期
w=2*3.14159*f; %角频率
j=0.1; %阻尼比
R=256;%输入信号长度
N=512;%FFT变换点数
num=;
den=;
sys=tf(num,den); %连续系统传递函数
=c2dm(num,den,ts,'tustin'); %双线性变换得离散方程
x=rand(1,R)-0.5; %产生（-0.5，0.5）均匀分布的随机噪声
y=filter(numd,dend,x); %随机噪声信号通过系统的响应
fx=fft(x,N);
fx_pos=fx(1:N/2);
fy=fft(y,N);
fy_pos=fy(1:N/2);
Hw1=fy_pos./fx_pos;
figure
Hw1f=abs(Hw1);
f_pos=;
plot(f_pos,Hw1f);
[ 本帖最后由 xyhan 于 2008-1-7 02:15 编辑 ]
页:
[1]