| 查看: 2575 | 回复: 20 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
跪求matlab的for循环语句的编写
|
|||
|
跪求各位大侠，指导一下下面这个程序的编写。 %需要输入的数组为x数组和T数组 x=[0.07 0.13 0.28 0.44] T=[283 296 312 320] %计算用到的函数如下，这是我随便写的函数，主要写明每个函数里面的未知输入量是哪个 m=x^2+2x n=2x^+3T o=4x^2+5T^2 p=m+n*o q=3/(T^2) %我要得到的结果如下 y(i)=exp(q-p) %最终目的是将计算得到的y(i)与代入的初值x(i)比较 k=y(i)-x(i) %如果k的绝对值abs(k)<=0.00001，则最终的输出结果为xcal=y(i)，否则将得出的y(i)赋值给最初代入计算的x，接着进行计算，直到abs(k)<=0.00001为止。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有13人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有14人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
【答案】应助回帖
|
数据会影响程序的，估计你要的是这个 clc clear x0=[0.07 0.13 0.28 0.44]; T=[283 296 312 320]; for i=1:length(x0) k(i)=1; y(i)=x0(i); while k(i)>0.00001 x(i)=y(i); m=x(i)^2+2*x(i); n=2*x(i)^2+3*T(i); o=4*x(i)^2+5*T(i)^2; p=m+n.*o; q=3./(T(i)^2); y(i)=exp(q-p); k(i)=abs(y(i)-x(i)); end end x k y 但是用这个也行 %需要输入的数组为x数组和T数组 function main clc clear x0=[0.07 0.13 0.28 0.44]; options = optimset('Display','iter','TolFun',1e-5); x=fsolve(@calk,x0,options) end %计算用到的函数如下，这是我随便写的函数，主要写明每个函数里面的未知输入量是哪个 function k=calk(x) T=[283 296 312 320]; m=x.^2+2*x; n=2*x.^2+3*T; o=4*x.^2+5*T.^2; p=m+n.*o; q=3./(T.^2); %我要得到的结果如下 y=exp(q-p); %最终目的是将计算得到的y(i)与代入的初值x(i)比较 k=abs(y-x)-0.00001; end %如果k的绝对值abs(k)<=0.00001，则最终的输出结果为xcal=y(i)，否则将得出的y(i)赋值给最初代入计算的x，接着进行计算，直到abs(k)<=0.00001为止。 |
17楼2012-06-04 10:27:39
2楼2012-06-03 17:10:55
yujia5038
铁杆木虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 8025.6
- 帖子: 86
- 在线: 178.2小时
- 虫号: 511726
- 注册: 2008-02-25
- 专业: 化学计量学与化学信息学
3楼2012-06-03 17:52:39
4楼2012-06-03 17:57:21
