Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

wokman1890 发表于 2006-5-27 16:18

菜鸟请教个菜鸟级的问题

小弟是这方面的菜鸟。我从某处粘贴了一个正确的程序到Matlab里运行，结果运行不出来，它提示说 Error: Missing variable or function. 这是哪里有错误呢？应该怎么改原来的程序才能正确运行呢？ 谢谢！！！

NVHht 发表于 2006-5-27 16:27

能不能把程序贴出来看看。

wokman1890 发表于 2006-5-27 16:40

可以哈来了 FREQZ_M.M function=freqz_m(b,a) %db为相对振幅(dB)；mag为绝对振幅；pha为相位响应；grd为群延时；w为频率样本点向量。 =freqz(b,a,1024, 'whole'); %b和a分别为H(z)分子和分母多项式系数(对FIR而言，b=a，a=1)． H=(H(1:501)),; W=(w(1:501)); mag=abs(H); db=20*logl0((mag+eps)/max(mag)); pha=angle(H); grd=grpdelay(b,a,w); IDEAL_LP.M function hd=ideal_lp(wc,M) %理想低通滤波器；hd为0到M一1之间的理想脉冲响应；wc为截止频率；M为理想滤波器的长度 alpha=(M-1)/2; n=0:M-1; m =n-alpha+eps; hd=sin(we*m)./(pi*m);
ACTUAL_LP.m wp=0.10 *pi; ws=0.20 *pi; tr_ width=ws-wp； M=ceil(6.6 *pi/tr_width)+1; disp(['滤波器的长度为',hum2str(M)]); n=0:M -1; wc=(ws+wp)/2; hd=IDEAL_LP(wc,M); W_han=(harming(M)); h=hd.* w_han; =FREQZ_M(h,); delta_w =2 *pi/lO00; Rp=-(min(db(1:1:wp/delta_w+1))); disp(['实际通带波动为',num2str(Rp)]); As= -round(max(db(ws/delta_w+1:1:501))); disp(['最小阻带衰减为', num2str(As)]); subplot(221) stem(n,hd); tide('理想冲激响应'); axis(); ylabel(hd(n)'); subplot(222) stem(n,w_han); title('汉宁窗'); axis(); ylabel( 'w(n)'), subplot(223) stem(n,h); title('实际冲激响应'); axis(); xlabel('n'); ylabel('h(n)'); subplot(224) plot(w/pi,db); title('幅度响应(dB)'); axis()， grid; xlabel('似pi为单位的频率'); ylabel('分贝数'); 这是个低通滤波器的例子我粘贴进去就运行不起

lxq 发表于 2006-5-27 18:25

是不是你存储的名字没有用英文呢？

wokman1890 发表于 2006-5-27 18:54

是要把存储程序的文件名改成FREQZ_M.M吗？

ericlin 发表于 2006-5-28 12:27

要新建一个M文件

你这个东西是程序，不是直接贴到MATLAB的运行窗口的，要新建一个M文件，以freqz_m程序名保存后再在命令窗口中运行。
[此贴子已经被作者于2006-5-28 12:30:06编辑过]

lixiang02 发表于 2006-5-28 16:17

你的m文件的文件名里面没有汉字吧?或者是不是用数字开头了?

wokman1890 发表于 2006-5-29 13:19

谢谢了哈那个问题解决了不过它又提示 =freqz(b,a,1024, 'whole');有错 怎么改呢？

NASA 发表于 2006-5-29 14:50

回复：（wokman1890）谢谢了哈那个问题解决了不过它...

<DIV class=quote>以下是引用wokman1890在2006-5-29 13:19:56的发言： 
谢谢了哈那个问题解决了不过它又提示 =freqz(b,a,1024, 'whole');有错 怎么改呢？ </DIV>
 freqz的调用格式如下，你的调用格式不对阿 = freqz(Hq,n) h = freqz(Hq,w) = freqz(Hq,n,'whole') = freqz(Hq,...) = freqz(Hq,n,fs) h = freqz(Hq,f,fs) = freqz(Hq,n,'whole',fs) = freqz(Hq,...) = freqz(Hq,...,fs) freqz(Hq,...)

页: [1]

声振论坛's Archiver

菜鸟请教个菜鸟级的问题

要新建一个M文件

回复：（wokman1890）谢谢了哈 那个问题解决了 不过它...

回复：（wokman1890）谢谢了哈那个问题解决了不过它...