求助如何根据表达式画幅频曲线

lise2001

我这里有个表达式，就是根据这个表达式画出图形。

而根据相关资料，我找到了参考图。

但是我编程后仿真图与参考图相差很大。

所以想请哪位高手帮忙指点下，哪里出了问题。

我是用.m文件编程的。
公式如下：

                               
登录/注册后可看大图

而参考图形如下：

                               
登录/注册后可看大图

对表达式进行说明：
该模型描述了典型的电力线信道的频率响应，概括了500KHz-20MHz的电力线的信道特征。


我的程序如下：echo on
t0=2.*10.^7;
fs=0.2;
ts=1/fs;
df=0.3;
f=[0:ts:t0];
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);
plot(f,m(1:length(f)))
xlabel('Frequency')
title('The Discrete Magnitude  Spectrum')


但是这个程序运行后的仿真图与参考图大不一样。。。


所以想请哪位高手帮忙指点下哪里出了问题？

[ 本帖最后由 eight 于 2008-4-28 18:20 编辑 ]

ch_j1985

初到请看置顶帖：
http://forum.vibunion.com/forum/thread-48179-1-1.html
另外，把图片放在附件里面就可以显示啦

lise2001

好的，谢谢楼上的提醒~~~~我会好好去弄下的哦~~

ch_j1985

原帖由 lise2001 于 2008-4-28 18:01 发表

                               
登录/注册后可看大图

好的，谢谢楼上的提醒~~~~我会好好去弄下的哦~~

我看了一下你的程序，画幅频特性时应该不带j的，你最好把图片和公式帖出来，不然代码不好修改:@)

lise2001

你好，很感谢你肯帮我
只是它说发帖时间与现在已经超过30分钟了
所以不能重新编辑
我可能要重新发下帖子。。。。

[ 本帖最后由 eight 于 2008-4-29 09:52 编辑 ]