Matlab 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Matlab
标  题: [合集]怎么进行加窗分析？
发信站: 哈工大紫丁香 (2004年04月10日08:54:28 星期六), 站内信件


────────────────────────────────────────
 XIMENGNI (丁丁)                      于 Wed Apr  7 13:46:13 2004 说道:


t=0:.1:6*pi;
x1=4*sin(2*pi*202.5*t) ;
x2=4*sin(2*pi*502.5*t)+4*sin(2*pi*514*t);
n=50;window=boxcar(n);[x2,w]=freqz(window,1);subplot(2,1,1);stem(window);subpl
ot(2,1,2);plot(w/pi,20*log(abs(x2)/abs(x2(1))));
我这么错在那里？

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2004年04月07日14:05:16 星期三 说道:

你应该把窗口加到信号上去

────────────────────────────────────────
 XIMENGNI (丁丁)                      于 Wed Apr  7 14:07:15 2004 说道:


怎么加阿

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2004年04月07日14:09:39 星期三 说道:

比如说对x2加窗：
xx=x2(1:50).*window;

────────────────────────────────────────
 Skyblue (原来的我)                   于 2004年04月07日14:12:08 星期三 说道:

                        ~~~~~~~~~~~~~~~这步做什么？
★另外，加矩形窗等于没有加窗！

────────────────────────────────────────
 XIMENGNI (丁丁)                      于 Wed Apr  7 14:18:47 2004 说道:


这样也不行。
 t=0:.1:6*pi;
x1=4*sin(2*pi*202.5*t) ;
x2=4*sin(2*pi*502.5*t)+4*sin(2*pi*514*t);
n=50;window=boxcar(n);xx1=x1(1:50)*window;xx2=x2(1:50)*window;
subplot(2,1,1);plot(t,xx1);subplot(2,1,2);plot(t,xx2);

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2004年04月07日14:19:59 星期三 说道:

注意那个小点点，不写可不行啊

────────────────────────────────────────
 XIMENGNI (丁丁)                      于 Wed Apr  7 14:27:36 2004 说道:


不行，那个点什么意思。不加就是两条线。

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2004年04月07日14:29:10 星期三 说道:

看看例子：
Hd = zeros(11,11); Hd(4:8,4:8) = 1;
[f1,f2] = freqspace(11,'meshgrid');
mesh(f1,f2,Hd), axis([-1 1 -1 1 0 1.2]), colormap(jet(64))
h = fwind1(Hd,hamming(11));
figure, freqz2(h,[32 32]), axis([-1 1 -1 1 0 1.2])

────────────────────────────────────────
 Skyblue (原来的我)                   于 2004年04月07日14:31:28 星期三 说道:

加点表示对应的项相乘，否则为矩阵相乘！
另外，应将window转置，t也应该等长。
最好加别的窗，矩形窗等于sample！
t=0:.1:6*pi;
x1=4*sin(2*pi*202.5*t) ;
x2=4*sin(2*pi*502.5*t)+4*sin(2*pi*514*t);
n=50;window=boxcar(n);xx1=x1(1:50).*window';xx2=x2(1:50).*window';
t=t(1:50);
subplot(2,1,1);plot(t,xx1);subplot(2,1,2);plot(t,xx2); 

────────────────────────────────────────