|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
信号时域表达式为:y=A.*exp(-((2.^1.5).*(t-t0)/t00).^2).*cos(2*pi*fd*(t-t0)),频域表达式为:Y=(A*t0/4).*(pi/4).^0.5.*(exp(-(pi*(f-fd)*t0)/(2.^1.5)).^2) +exp(-(pi*(f+fd)*t0)/(2.^1.5)).^2) ) ,频域包含两个变量,一个是(f-fd),另一个是(f+fd),我用窗函数想把高频部分滤掉不行,不知道怎么才能把高频部分滤掉哪?请高人指教,谢谢!!
程序如下:
clear
sf=2500000; %采样频率
deltat=1/sf; N=16384;
t=0:deltat:(N*deltat-deltat);
deltaf=sf/N;
fd=600000; Nf=20;% 多普勒频率 测量区条纹数
A=1;%B=0.000001; %信号幅值
t0=0.0015; t00=Nf/fd; %t0信号到达时间
fp1=800000; fp2=900000; %高通滤波器
wp=2*fp1/sf; %阻带截止频率
ws=2*fp2/sf; %通带起始频率
deltaw=ws-wp; %确定过度带宽
N0=ceil(6.2*pi/deltaw); %按汉宁窗计算滤波器长度N0,ceil是向上取整函数
N1=N0+mod(N0+1,2); %为实现第一类对偶滤波器,要确保其长度N1为奇数
n=N1-1; %滤波器的阶数
wn=(ws+wp)/2; %理想高通的截止频率
b=fir1(n,wn,'high',hanning(N1));
y=A.*exp(-((2.^1.5).*(t-t0)/t00).^2).*cos(2*pi*fd*(t-t0));
%单个粒子产生的光电流
y1 = filter(b,1,y);
y2=abs(fft(y));
y3=abs(fft(y1));
subplot(211)
plot(0:(sf/N):(sf-sf/N),y2);
subplot(212)
plot(0:(sf/N):(sf-sf/N),y3);
[ 本帖最后由 eight 于 2008-4-1 10:06 编辑 ] |
|