24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

askformore

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1641.5
散金: 55
红花: 3
帖子: 134
在线: 129.2小时
虫号: 859087
注册: 2009-09-28
专业: 天体中基本物理过程的理论

[求助] 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获得目标图形，可以再追加金币，多谢！

目标图形

为什么只显示后半段呢？

回复此楼

» 猜你喜欢

上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

继续求助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人回复

1楼 2012-05-22 21:03:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hhucy

实习版主 (著名写手)

DOE锅炉工

应助: 53 (初中生)
贵宾: 0.264
金币: 7775.6
散金: 1861
红花: 13
沙发: 7
帖子: 2588
在线: 1275.5小时
虫号: 1079650
注册: 2010-08-22
性别: GG
管辖: 能源

引用回帖:

8楼: Originally posted by 疾凤劲草 at 2012-05-25 16:59:26
楼主啊，你的程序我前几天用我的matlab r2010a算了一下很久都没有反映，那天我也在编程，运行了一个很多循环的matlab程序，以为matlab给弄出问题了。今天我又运行了一下，还是很久都没有反映，但是关闭matlab主窗口 ...

是的，新版本matlab会这样的