完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
扫一扫,分享给好友
本帖最后由 zwj108 于 2012-11-27 18:22 编辑
n=0:3; xb=impesq(0,0,3); %用于产生xb(n)=δ(n)信号 hb=impesq(0,0,3)+2.5*impesq(1,0,3)+2.5*impesq(2,0,3)+impesq(3,0,3); %用于产生hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) k=-200:200; w=(pi/100)*k; Xb=xb*(exp(-1i*pi/100)).^(n'*k); %信号xb的傅氏变换 Hb=hb*(exp(-1i*pi/100)).^(n'*k); %信号hb的傅氏变换 n=0:3; subplot(3,2,1); %做序列xb的时域和频域的图形 stem(n,xb); axis([-2 2 0 2]); xlabel('n'); ylabel('xb(n)'); title('xb(n)'); subplot(3,2,2); plot(w/pi,abs(Xb)); axis([-2 2 0 2]); xlabel('w/pi'); ylabel('xb(|(jw)|'); title('[xb(ejw)]'); subplot(3,2,3); %做序列hb的时域和频域的图形 stem(n,hb); axis([0 4 0 3]); xlabel('n'); ylabel('Hb'); title('Hb(n)'); subplot(3,2,4); plot(w/pi,abs(Hb)); axis([-2 2 0 8]); xlabel('w/pi'); ylabel('Hb(|(jw)|'); title('[Hb(ejw)]'); n=0:6; y=conv(xb,hb); %序列xb,hb的卷积 Y=y*(exp(-1i*pi/100)).^(n'*k); %卷积的傅氏变换 subplot(3,2,5); %响应y的时域和频域的图形 stem(n,y); %卷积作图 axis([0 7 0 3]); xlabel('n'); ylabel('y(n)'); title('xb*hb'); subplot(3,2,6); plot(w/pi,abs(Y)); %傅氏变换作图 axis([-2 2 0 8]); xlabel('w/pi'); ylabel('|Y(jw)|'); title('[Y(ejw)]'); ??? Error using ==> conv2 First and second arguments must be single or double. Error in ==> conv at 39 c = conv2(a(:),b(:),shape); |
|
相关推荐
3个回答
|
|
搞不懂这种高级东西,帮顶
|
|
|
|
|
|
|
|
|
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-8-23 23:31 , Processed in 0.412296 second(s), Total 48, Slave 41 queries .
Powered by 电子发烧友网
© 2015 bbs.elecfans.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191