| 查看: 420 | 回复: 2 | |||
[求助]
matlab编程
|
|
T=1200:10:1400; a=exp(-(2.98.*log(T)+4.*T/100000+12500./T.^2+47787./T-52.01)./1.987); b=exp(-(5.79.*log(T)+1.75.*T/10000+48150./T.^2+92542./T-116.93)./1.987); c=exp(-(-0.42.*log(T)+1.1.*T/10000-89000./T.^2+114925./T-76.16)./1.987); p=[1 0 -b -2.*(a.*b-c.^2)]; A=roots(p) 求出根后，我想画出其中任意一个根随T变化的图像，该怎么编程 |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
心脉受损
已经有6人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
johnliu1983
至尊木虫 (著名写手)
- 应助: 43 (小学生)
- 金币: 15097.3
- 红花: 8
- 帖子: 2434
- 在线: 149.9小时
- 虫号: 370817
- 注册: 2007-05-14
- 性别: GG
- 专业: 物理学II
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2013-03-02 14:12:10
感谢参与，应助指数 +1
jjdg: 金币+1, 辛苦了 2013-03-02 14:12:10
|
Tvalue=1200:10:1400 for k=1:length(Tvalue) T=Tvalue(k); a=exp(-(2.98.*log(T)+4.*T/100000+12500./T.^2+47787./T-52.01)./1.987); b=exp(-(5.79.*log(T)+1.75.*T/10000+48150./T.^2+92542./T-116.93)./1.987); c=exp(-(-0.42.*log(T)+1.1.*T/10000-89000./T.^2+114925./T-76.16)./1.987); p=[1 0 -b -2.*(a.*b-c.^2)]; A=roots(p); target(k,1)=A(1); target(k,2)=A(2); target(k,3)=A(3); end figure(1) plot(Tvalue',target(:,1),'r') figure(2) plot(Tvalue',target(:,2),'g') figure(3) plot(Tvalue',target(:,3),'k') 图一到三是三个解随T的变化图像 |
2楼2013-03-02 01:20:21
3楼2013-03-02 13:21:21