| 查看: 477 | 回复: 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变化的图像，该怎么编程 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
johnliu1983
至尊木虫 (著名写手)
- 应助: 43 (小学生)
- 金币: 15285.3
- 红花: 8
- 帖子: 2474
- 在线: 150.8小时
- 虫号: 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
