| 查看: 2833 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
askformore木虫 (小有名气)
|
[求助]
matlab画图求助
|
||
|
需要用matlab画出下面这个图形【即第一个图形】： 表达式都是很清楚的，只是需要用matlab画出来，我写的程序是下面这样，很简单的： for y=18:0.001:21; k=[2:0.01:25]; a1=pi*(k-2).^4/12./(1+0.8048.*(k-2)+0.1459.*(k-2).^2+1.137*10^(-3)*(k-2).^3-3.879*10^(-6)*(k-2).^4)./(exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))-1); b1=trapz(k,a1); da1=pi*(k-2).^4/12./(1+0.8048.*(k-2)+0.1459.*(k-2).^2+1.137*10^(-3)*(k-2).^3-3.879*10^(-6)*(k-2).^4).*(-3./(exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))-1)+0.511*938.272.*k.* exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))./(20*10^(y-18)*2.7*8.6170)./(exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))-1).^2); db1=trapz(k,da1); c=[25:0.1:10000]'; a2=c.*(-86.07+50.96.*log(c)-14.45.*log(c).^2+8/3.*log(c).^3)./(1-(2.91./c+78.35./c.^2+1837./c.^3))./(exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))-1); b2=trapz(c,a2); da2=c.*(-86.07+50.96.*log(c)-14.45.*log(c).^2+8/3.*log(c).^3)./(1-(2.91./c+78.35./c.^2+1837./c.^3)).*(-3./(exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))-1)+0.511*938.272.*c.*exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))./(20*10^(y-18)*2.7*8.6170)./ (exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))-1).^2); db2=trapz(c,da2); dp=(db1+db2)./(b1+b2); plot(y,dp,'r'); axis([18 21 -1 3]); hold on; end 得到的图形为【第二个图形，即红色曲线的图形】： 很明显，和目标图形相比，matlab只画出了整个图形横坐标在大约y>19.2之后的图形，而且，这一段图形和目标图形是一致的【或者说一样的】。但是，为什么y<19.2的部分画不出来呢？ 请各位同学帮助哈，如果能帮助修改，用matlab获得目标图形，可以再追加金币，多谢！  目标图形  为什么只显示后半段呢？ |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有18人回复
-
基金系统什么内容也没有
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有12人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
面上合作单位盖章
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 继续求助Matlab绘图！
已经有8人回复
- 求助如何用matlab取出直线上的点
已经有10人回复
- 求助matlab三维画图中：Inner matrix dimensions must agree.
已经有4人回复
- 求助：谁会用matlab画图？
已经有13人回复
- matlab 画图求助
已经有5人回复
- 【求助】matlab话纵向图
已经有5人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】初学matlab
已经有11人回复
- 【求助】柱形统计图matlab怎么画？
已经有4人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
askformore: 金币+50, ★★★★★最佳答案, 非常感谢您耐心的修改，运行您修改的程序，已得到目标图形。辛苦辛苦！ 2012-05-28 18:45:38
askformore: 金币+50, ★★★★★最佳答案, 非常感谢您耐心的修改，运行您修改的程序，已得到目标图形。辛苦辛苦！ 2012-05-28 18:45:38
|
楼主的问题出现在inf/inf=nan，如楼上所说，exp(700)左右就达到了inf,解决的方法是将两个inf进行约分，概括起来即对exp(x)/(exp(x)-1)^2进行约分。以下是我修改后的程序： clc clear for y=18:0.001:21; k=[2:0.01:25]; a1=pi*(k-2).^4/12./(1+0.8048.*(k-2)+0.1459.*(k-2).^2+1.137*10^(-3)*(k-2).^3-3.879*10^(-6)*(k-2).^4)./(exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))-1); b1=trapz(k,a1); da1=pi*(k-2).^4/12./(1+0.8048.*(k-2)+0.1459.*(k-2).^2+1.137*10^(-3)*(k-2).^3-3.879*10^(-6)*(k-2).^4).*(-3./(exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))-1)+0.511*938.272.*k.* exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))./(20*10^(y-18)*2.7*8.6170)./(exp(0.511*938.272.*k./(20*10^(y-18)*2.7*8.6170))-1).^2); db1=trapz(k,da1); c=[25:0.1:10000]'; a2=c.*(-86.07+50.96.*log(c)-14.45.*log(c).^2+8/3.*log(c).^3)./(1-(2.91./c+78.35./c.^2+1837./c.^3))./(exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))-1); b2=trapz(c,a2); da2=c.*(-86.07+50.96.*log(c)-14.45.*log(c).^2+8/3.*log(c).^3)./(1-(2.91./c+78.35./c.^2+1837./c.^3)).*(-3./(exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))-1)+0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170)./(exp(0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))+exp(-0.511*938.272.*c./(20*10^(y-18)*2.7*8.6170))-2)); db2=trapz(c,da2); dp=(db1+db2)./(b1+b2); plot(y,dp,'r'); axis([18 21 -1 3]); hold on; end 直接运行就可以，图我已经得到，不好上传，需要的话再联系  |
10楼2012-05-26 13:40:35
3楼2012-05-23 09:47:21
4楼2012-05-23 09:48:48
5楼2012-05-23 12:00:14
投票:
