| 查看: 2371 | 回复: 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为止。 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
【答案】应助回帖
|
数据会影响程序的，估计你要的是这个 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
25