| 查看: 693 | 回复: 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 |
回复此楼» 猜你喜欢
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
假如你的研究生提出不合理要求
已经有6人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
pdl9527
专家顾问 (小有名气)
-
专家经验: +8 - 应助: 100 (初中生)
- 金币: 2107.8
- 红花: 23
- 帖子: 282
- 在线: 111.9小时
- 虫号: 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