| 查看: 812 | 回复: 1 | |||||
[求助]
将计算值不断输入一个矩阵当中，求最新的输入矩阵的负值应该怎么表示？ 已有1人参与
|
|
求教大家一个问题，我在编一个很简单的小程序，其实就是二分法求解。现在生成一个矩阵p，将计算的p值输到矩阵p中，现在我将最新输入到矩阵p中的正值（或负值）取出，请问该怎么办？谢谢大家。 另外附上我的二分法部分的程序麻烦大家看一下： cn(i)=A-B if (((cn(i)>=0 & cn(i-1)>=0)) p=(p+pmax)/2 elseif (cn(i)<=0 & cn(i-1)<=0) p=(0+p)/2 else p=(pn(i)+pn(i-1))/2 end 其实已经知道这个程序哪里不太对，应该改成： if ((cn(i)>=0 & cn(i-1)>=0)) p=(p+上一个取正值的cn(i)时对于的p)/2 elseif (cn(i)<=0 & cn(i-1)<=0) p=(0+上一个取负值的cn(i)时对于的p)/2 else p=(pn(i)+pn(i-1))/2 end |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有19人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2110.3
- 红花: 23
- 帖子: 282
- 在线: 112小时
- 虫号: 1227333
- 注册: 2011-03-09
- 性别: GG
- 专业: 分离过程
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
哈格不受伤: 金币+20, ★★★★★最佳答案 2016-12-16 16:34:11
感谢参与，应助指数 +1
哈格不受伤: 金币+20, ★★★★★最佳答案 2016-12-16 16:34:11
|
function [k,x,wuca,yx]=erfen(a,b,abtol) a(1)=a; b(1)=b; ya=fun(a(1)); yb=fun(b(1)); %程序中调用的fun.m 为函数 if ya* yb>0, disp('注意：ya*yb>0,请重新调整区间端点a和b.'), return end max1=-1+ceil((log(b-a)- log(abtol))/ log(2)); % ceil是向+ 方向取整 for k=1: max1+1 a;ya=fun(a); b;yb=fun(b); x=(a+b)/2; yx=fun(x); wuca=abs(b-a)/2; k=k-1; [k,a,b,x,wuca,ya,yb,yx] if yx==0 a=x; b=x; elseif yb*yx>0 b=x;yb=yx; else a=x; ya=yx; end if b-a< abtol , return, end end k=max1; x; wuca; yx=fun(x); |
2楼2016-12-02 04:01:03
