| 查看: 2383 | 回复: 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为止。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有85人回复
-
感觉是下周放榜了
已经有9人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
|
x=[0.07 0.13 0.28 0.44]; T=[283 296 312 320]; xcal=ones(1,4); t=find(abs(xcal)==1); >> l=length(t); >> while (l>0) for i=1:4 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=abs(y(i)-x(i)); if k<=0.01 xcal(i)=y(i); else x(i)=y(i); end end t=find(abs(xcal)==1); l=length(t); end x y xcal |
21楼2012-06-05 08:45:18
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
