我自己写了个工具包，但是有问题，帮忙看一下

石头王石头

  function [X]=dct(xn)
N=length(xn);
sum=0;
for n=1:N
    sum=sum+xn(n);
end
X(0)=sum2*sqrt(N);
for k=1:N-1
    sum2=0;
    for n=1:N
        sum2=sum2+xn(n)*cos((2*n+1)*k*pi/(2*N));
    end
    X(k)=sum2*sqrt(2/N);
end
不胜感激

ChaChing

1.与工具箱有关係?
2.若有报错的,应给齐讯息; 若没报错的就不是编程问题吧

石头王石头

我用上面的程序做工具包，写的下面程序
clc;
clear;
f=50;
fs=1000;
N=1000;
n=0:N-1;
xn=cos(2*pi*f*n/fs);
y=dct(xn);
num=find(abs(y)<5);
y(num)=zeros(size(num));
zn=idct(y);
figure(2);
subplot(211),plot(n,xn);
subplot(212),plot(n,zn);
rp=100-norm(xn-zn)/norm(xn)*100；

结果出现错误

??? Attempt to execute SCRIPT dct as a function:
D:\My Documents\MATLAB\dct.m

Error in ==> dct at 8
y=dct(xn);
是不是我开始写的工具包错了

ChaChing

石头王石头 发表于 2013-7-3 09:52

                               
登录/注册后可看大图

我用上面的程序做工具包，写的下面程序
clc;
clear;

看到报错讯息, 可以确定LZ方向错误!

工具包是什麼!? 不懂
但1F的m文件, 会报错

??? Undefined function or variable "sum2".

Error in ==> dct at 7
X(0)=sum2*sqrt(N);

不过LZ的错误可能还出现在更前端!?
看看下帖6F
常见的程序出错问题整理 http://www.chinavib.com/thread-46001-1-1.html
(from http://home.chinavib.com/blog-63979-18250.html)

石头王石头

ChaChing 发表于 2013-7-3 20:24

                               
登录/注册后可看大图

看到报错讯息, 可以确定LZ方向错误!

工具包是什麼!? 不懂

谢谢了，终于知道哪错了
仔细查找了下原因是因为文件名取错了改为dct1或者其它都可以了，不能用matlab内部自带的函数名来命名。

ChaChing

石头王石头 发表于 2013-7-5 20:24

                               
登录/注册后可看大图

谢谢了，终于知道哪错了
仔细查找了下原因是因为文件名取错了改为dct1或者其它都可以了，不能用matlab内 ...

...不能用matlab内部自带的函数名来命名...

严谨而言, 是能以内部自带的函数名来命名的, 不过若不清楚如何控制的话, 使用不同名是个安全的做法

zswseu

石头王石头 发表于 2013-7-5 20:24

                               
登录/注册后可看大图

谢谢了，终于知道哪错了
仔细查找了下原因是因为文件名取错了改为dct1或者其它都可以了，不能用matlab内 ...

楼上说的很有道理。
我记得如果把函数定义为子函数。函数名和matlab内部函数重名就不要紧了
因为子函数优先级高。你可以试一下。