以下是我从网上下载的一段程序，用MATLAB运行后。得不到别人所说的结果，请高手指点。

发布网友发布时间：2023-01-01 21:49

共3个回答

热心网友时间：2023-10-08 20:02

不知道你所谓别人所说的结果是什么，我运行了结果如下图。

另外如果你把上面的整段全部在一个脚本文件运行会有如下报错：

??? Error: File: F:\Matlab7\work\test57.m Line: 54 Column: 1

Function definitions are not permitted at the prompt or in scripts.

意思是函数不能在命令窗口或者脚本文件中定义。

只要把函数freqz_m定义的部分剪切在另一个M文件中保存就可以了（文件明须为freqz_m.m）

%所使用的M文件函数

function[db,mag,pha,grd,w]=freqz_m(b,a)

%滤波器的幅值响应(相对、绝对)、相位响应及群延迟

%Usage:[db,mag,pha,grd,w]=freqz_m(b,a) %500点对应[0,pi]

%w 采样频率； b系统函数H(z)的分子项(对FIR，b=h)

%a 系统函数H(z)的分母项(对FIR，a=1)

[H,w]=freqz(b,a,500); %500点的复频响应

mag=abs(H); %绝对幅值响应

db=20*log10(mag/max(mag)); %相对幅值响应

pha=angle(H); %相位响应

grd=grpdelay(b,a,w); %群延迟响应

热心网友时间：2023-10-08 20:02

clear Ac=4e-6; f=3.62e4; v=1/60; t=0.00035:0.00000001:0.size(xt,2)有错，xt是一维数据要输出sx和sy，程序中少写 sx sy 结果全

热心网友时间：2023-10-08 20:02

我的问题和你一样，请问你现在解决了吗？