Numerical solution

Question

0 个投票

Hello

I´m working on a numerical solution for an equation. Therefore I wrote the following function:

function test(dBL,i)

 for n=1:i
    BL=0;

while (1+cos(BL)*cosh(BL))*(1+cos(BL+dBL)*cosh(BL+dBL))>0

   BL=BL+dBL;
fc=1+cos(BL)*cosh(BL);
end

vfc(n,1)=fc;

vBL(n,1)=BL;

end

vfc

vBL

end

My aim is that the matlab function produce as many solutions as n and the important output is vBL, which is a vector with the solutions.

I want the while loop to start with BL=0 for n=1, but for n=2 I want the loop to start with the last BL from the previous while loop and so on. Anybody can help??

Thanks in advance

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

Answer 1

Dr. Seis 2012-3-6

在 MATLAB Online 中打开

0 个投票

Unless I misunderstood something, I think all you need to do is change:

BL = 0;

to

if n == 1
    BL = 0;
end

Then for n > 1, "BL" will be whatever value it was at the end of the previous iteration.

In your case (now that I am looking at it more closely), if "BL" starts off as the same value as the previous iteration, then the "while" loop probably will not execute... since "BL" and "dBL" are not changing, (1+cos(BL)*cosh(BL))*(1+cos(BL+dBL)*cosh(BL+dBL))>0 will also not change. Should that expression be some sort of function of "n" as well?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Ali Ihsan 2012-3-6

Hello Elige

I found out a solution. I did exactly what you mentioned, what are the odds :)

The function is now:

function test(dBL,i)

for n=1:i

if n==1

BL=0;

elseif n>1

BL=result+dBL;

end

while (1+cos(BL)*cosh(BL))*(1+cos(BL+dBL)*cosh(BL+dBL))>0

BL=BL+dBL;

fc=1+cos(BL)*cosh(BL);

end

result=BL;

vBL(n,1)=BL;

end

vBL

end

By adding dBL to the result, the loop starts in the other side of the BL axes and it works now.

Thank you very much.

请先登录，再进行评论。

Numerical solution

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

更多回答（0 个）

类别

标签

Community Treasure Hunt

Numerical solution

0 个评论 显示 -2更早的评论 隐藏 -2更早的评论

采纳的回答

1 个评论 显示 -1更早的评论 隐藏 -1更早的评论

更多回答（0 个）

类别

标签

另请参阅

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论