声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1524|回复: 4

[综合讨论] 请帮我看一下这个程序

[复制链接]
发表于 2005-12-21 18:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
clear;wp1=0.2*pi;ws1=0.25*pi;ws2=0.30*pi;wp2=0.35*pi;<BR>rp=2;rs=50;<BR>tr_width=min((ws1-wp1),(wp2-ws2));<BR>M=ceil(6.6*pi/tr_width)+1<BR>n=[0:1:M-1];<BR>wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2;<BR>hd=ideal_lp(wc1,M)+ideal_lp(pi,M)-ideal_lp(wc2,M);<BR>w_bla=(blackman(M));<BR><FONT color=#0909f7>h=hd.*w_bla</FONT>;<BR>[db,mag,pha,grd,w]=freqz_m(h,[1]);<BR>delta_w=2*pi/1000;<BR>rp=-min(db(wp1/delta_w+1:1:wp2/delta_w))<BR>rs=-round(max(db(ws2/delta_w+1:1:ws1/delta_w)))<BR>subplot(2,3,1);stem(n,hd);title('理想脉冲响应')<BR>axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('gd(n)');<BR>subplot(2,3,2);stem(n,w_box);title('矩形窗')<BR>axis([0 M-1 0 1.1]);xlabel('n');ylabel('w(n)');<BR>subplot(2,3,3);stem(n,h);title('实际脉冲响应')<BR>axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('h(n)');<BR>subplot(2,3,4);plot(w/pi,db);axis([0 1 -100 10]);<BR>title('幅度响应(单位:dB)');grid;<BR><BR>有错,但我不知如何改

[ 本帖最后由 eight 于 2008-5-7 15:46 编辑 ]
回复
分享到:

使用道具 举报

发表于 2005-12-21 19:15 | 显示全部楼层
楼上,matlab提示是什么错误?这有助于别人调试,对比,解决你的问题
发表于 2005-12-21 19:42 | 显示全部楼层
<P>你运行是这个错误吗?<BR>“<BR>&gt;&gt; <BR>M =</P>
<P>   133</P>
<P>??? Undefined function or variable 'ideal_lp'.</P>
<P>Error in ==&gt; d:\MATLAB6p5\work\haha.m<BR>On line 7  ==&gt; hd=ideal_lp(wc1,M)+ideal_lp(pi,M)-ideal_lp(wc2,M);</P>
<P>&gt;&gt; <BR>”<BR>如果是的话就是<BR>1。因为你调用了一个未定义的子函数'ideal_lp';<BR>2。因为虽然你定义了,但未放到同一目录下<BR></P>
发表于 2008-5-7 14:13 | 显示全部楼层
那数'ideal_lp'应该怎么定义呢?
发表于 2008-5-7 15:45 | 显示全部楼层
这个应该是你自己定义的函数,或者你从哪里拷来的程序,你再好好找找。

[ 本帖最后由 eight 于 2008-5-7 16:31 编辑 ]
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-18 04:39 , Processed in 0.063357 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表