Matlab 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Matlab
标  题: [合集]困扰几天的问题
发信站: 哈工大紫丁香 (2003年05月10日15:15:26 星期六), 站内信件


────────────────────────────────────────
 few (红叶)                           于 2003年05月06日21:51:54 星期二 说道:

一维信号x做如下变换：
x2=cceps(x);
x3=icceps(x);
为什么x和x3相差甚远

────────────────────────────────────────
 few (红叶)                           于 2003年05月06日21:58:02 星期二 说道:

好像是进行了相位调整，但不知道是怎么调整的，可能实现分段了，但又不知分为多少段

────────────────────────────────────────
 Systems (Control Nerd)               于 2003年05月06日22:00:08 星期二 说道:

许多人都遇到过类似的问题.俺猜想是由于误差造成的.
俺明天做个实验分析一下再说

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2003年05月06日22:01:57 星期二 说道:

这样做呢？
[x2,nd]=cceps(x);
x3=icceps(x2,nd);
我试了一下，
x和x3不是差很远呀

────────────────────────────────────────
 Systems (Control Nerd)               于 2003年05月06日22:05:59 星期二 说道:

nod,这可能也是一个原因,就是反变换取的点的稠密程度不一样.
那你再试试手动设置点数如何?
你看,最近很多人对类似的问题有困惑,应该做个实验看看

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2003年05月06日22:11:00 星期二 说道:

可能是这样，做如下变动就不行了
[x2,nd]=cceps(x);
x3=icceps(x2,nd+1);
x和x3就有差别了

────────────────────────────────────────
 Systems (Control Nerd)               于 2003年05月06日22:19:00 星期二 说道:

你误会俺的意思了,俺是说你最好找找根本原因.
俺手头没有MATLAB,只能凭经验猜.你可以做做实验:)

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2003年05月06日22:23:30 星期二 说道:

呵呵，我再看看，
只是我对这个不熟

────────────────────────────────────────
 few (红叶)                           于 2003年05月06日22:25:42 星期二 说道:

cceps是求一位信号的复倒谱，即先做复立叶变换然后取对数最后再做复里叶逆变换，

────────────────────────────────────────
 Systems (Control Nerd)               于 2003年05月06日22:26:41 星期二 说道:

所以俺猜是两个原因
1正向和反向变换的点的稠密度不一样,正如你所说
2.MATLAB由于截断误差,造成的变换的误差积累,如果
  遇到除法什么的,可能产生奇怪的值.
你不如找个简单的序列试试,比如做做8个点的cos的FFT
正反变换,这样容易看出毛病

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2003年05月07日09:01:13 星期三 说道:

我觉得应是第一个原因，至于截断误差应该很小，
对于输入信号应该可以忽略截断误差

────────────────────────────────────────
 zjliu (秋天的萝卜)                   于 2003年05月07日09:15:12 星期三 说道:

我又得到这样一个结果
比如对于这样的两个信号
x=cos(0:.1:7);
x=cos(0:.1:7)+tan(0:.1:7);
对它们分别进行
如下变换和逆变换：
[x2,nd]=cceps(x);
x3=icceps(x2,nd);
结果前者的结果是x和x3相等，后者存在误差

────────────────────────────────────────
 SELL (赛尔)                          于 2003年05月07日15:48:54 星期三 说道:

复倒谱变换中,存在两个重要的问题,1,线性相位项问题,2,相位解卷问题.
我不知道你原来的信号是不是存在噪声等问题 .
复倒谱对噪声敏感

────────────────────────────────────────
 SELL (赛尔)                          于 2003年05月07日16:04:34 星期三 说道:

问一个问题,你用复倒谱做什么,我也在用

────────────────────────────────────────
 few (红叶)                           于 2003年05月07日16:15:59 星期三 说道:

你是无系的吗

────────────────────────────────────────