如何根据已给出的表达式画图（附图）
求助如何根据表达式画幅频曲线我这里有个表达式，就是根据这个表达式画出图形。而根据相关资料，我找到了参考图。
但是我编程后仿真图与参考图相差很大。
所以想请哪位高手帮忙指点下，哪里出了问题。
我是用.m文件编程的。
表达式如下：
而参考图形如下：
对表达式进行说明：
该模型描述了典型的电力线信道的频率响应，概括了500KHz-20MHz的电力线的信道特征。
我的程序如下：echo on
t0=2.*10.^7;
ts=5000;
f=;
m=0.64.*exp(-1.56.*(10.^(-7)).*f).*exp(-2.66.*j.*pi.*f)+0.38.*exp(-1.73.*(10.^(-7)).*f).*exp(-2.97.*j.*pi.*f)-0.15.*exp(-1.9.*(10.^(-7)).*f).*exp(-3.26.*j.*pi.*f)+0.05.*exp(-2.09.*(10.^(-7)).*f).*exp(-3.57.*j.*pi.*f);
absm=abs(m);
plot(f,absm)
xlabel('Frequency')
title('The Discrete MagnitudeSpectrum')
但是这个程序运行后的仿真图与参考图大不一样。。。
所以想请哪位高手帮忙指点下哪里出了问题？
[ 本帖最后由 eight 于 2008-4-29 09:53 编辑 ]
回复 楼主 的帖子
你的坐标和图中的坐标不一样吧？你的要比他大 楼主的matlab程序我看了一下，没有问题，虽然用了很多不必要的'.'，但是结果应该是没问题的。
我觉得可能还是楼主的表达式和图形对应关系有问题——或者说那个图并不是由你给出的那个函数直接得出来的，或许经过了什么变换？ 采样点间隔太大
回复 4楼
我用ts=1000绘图和ts=5000绘图差异不大，而才用ts=500时出现了以下的图形,高频率的振荡，不知道为什么。[ 本帖最后由 sogooda 于 2008-4-29 10:43 编辑 ] 原帖由 sogooda 于 2008-4-29 09:23 发表 http://www.chinavib.com/forum/images/common/back.gif
我用ts=1000绘图和ts=5000绘图差异不大，而才用ts=500时出现了一下的图形,高频率的振荡，不知道为什么。
我画了一下，跟你画的图差不多。看了一下m的值都非常接近，所以我觉得参考图的纵坐标是经过处理的，比如取对数或用分贝表示的，希望楼主再仔细想想
[ 本帖最后由 ch_j1985 于 2008-4-29 09:50 编辑 ] 原帖由 sigma665 于 2008-4-29 08:43 发表 http://www.chinavib.com/forum/images/common/back.gif
你的坐标和图中的坐标不一样吧？
你的要比他大
谢谢你的答复，只是我看到我这个信道的频率范围是500KHz-20MHz，所以我才取到2x10^7 原帖由 sogooda 于 2008-4-29 08:59 发表 http://www.chinavib.com/forum/images/common/back.gif
楼主的matlab程序我看了一下，没有问题，虽然用了很多不必要的'.'，但是结果应该是没问题的。
我觉得可能还是楼主的表达式和图形对应关系有问题——或者说那个图并不是由你给出的那个函数直接得出来的，或许经过了什 ...
我也觉得应该是经过了一定的变换
它这个参考图的仿真工具是COSSAP
而我的是MATLAB
我怀疑它这个是有一个输入信号，因为这个幅度这么大
而这个表达式给出的权重Gi数值只有0-1，根本没那么大的 原帖由 huright 于 2008-4-29 09:09 发表 http://www.chinavib.com/forum/images/common/back.gif
采样点间隔太大
采样点，我是根据它那个参考图来考虑的
只是，原本我取的采样点也是蛮小的
但是很平滑。。。
这个表达式是一个代表“衰减，时延”等影响的一个信道模型 原帖由 sogooda 于 2008-4-29 09:23 发表 http://www.chinavib.com/forum/images/common/back.gif
我用ts=1000绘图和ts=5000绘图差异不大，而才用ts=500时出现了以下的图形,高频率的振荡，不知道为什么。
这个情况我也出现过
可能这个模型只只合于低频吧 原帖由 ch_j1985 于 2008-4-29 09:49 发表 http://www.chinavib.com/forum/images/common/back.gif
我画了一下，跟你画的图差不多。看了一下m的值都非常接近，所以我觉得参考图的纵坐标是经过处理的，比如取对数或用分贝表示的，希望楼主再仔细想想
谢谢你的提议~
其实我也想过取分贝
我怀疑它这个参考图肯定取过分贝
因为我们无论怎么取值，画出的幅値永远都是正值
但是参考图里却是负值...
我想请教下，怎么取分贝？？用什么函数？
这个我们当时没学过。。。 原帖由 lise2001 于 2008-4-30 00:30 发表 http://www.chinavib.com/forum/images/common/back.gif
我也觉得应该是经过了一定的变换
它这个参考图的仿真工具是COSSAP
而我的是MATLAB
我怀疑它这个是有一个输入信号，因为这个幅度这么大
而这个表达式给出的权重Gi数值只有0-1，根本没那么大的
如果是有输入信号的话，就可能了。
同时在不知道输入的前提下，也就没法画出和示范图一样的图了。
[ 本帖最后由 sogooda 于 2008-4-30 08:42 编辑 ]
回复 12楼 的帖子
怪不得啊~~谢谢你的答复~~我等下加个输入信号进去试试看~~ 原帖由 lise2001 于 2008-4-30 00:40 发表 http://www.chinavib.com/forum/images/common/back.gif
谢谢你的提议~
其实我也想过取分贝
我怀疑它这个参考图肯定取过分贝
因为我们无论怎么取值，画出的幅値永远都是正值
但是参考图里却是负值...
我想请教下，怎么取分贝？？用什么函数？
这个我们 ...
《自动控制原理》中分贝是这样求的：
L(w)=20lg|G(jw)|=20lgA(w) 单位(dB),其中G(jw)为频率特性,A(w) 为幅频特性 谢谢~~现在出现了另一个问题。。。
那就是我的输入信号是一个sincx函数
它的周期很小，我希望这个输入信号经过这个信道，但是信道带宽很大
所以想产生一个连续的sincx信号，不知道怎么产生
页:
[1]
2