for loopの計算中にタイムステップの変更

Takashi Miyatsuji

2018 9 6

0 个回答

3 次查看（30 天）

0 个投票

計算途中にタイムステップが流動的に変更されるコードを作成したいのですがどうすればいいでしょうか。具体的には、時間の変数を指定し、１ステップを50秒としてfor loopを回し、ある変数（C）が0未満になってしまったとき、そのステップで変数Cが０になる時間変数で計算をやり直す　というプログラムを作成したいと考えております。どうかお力添えをお願いいたします。

4 个评论
显示 2更早的评论隐藏 2更早的评论

Takashi Miyatsuji 2018-9-7

在 MATLAB Online 中打开

一部ですが、記載します。
dt =  50;
    for i = 1:18
      eval(['C',num2str(i),' = zeros(M+1,N+1);'])
      end
for j  = 2 : M
     for i = 2 : N  
 C1(j,1) = 5.55     ; 
 C2(j,1) = 168.75*sin(j*2*-pi/(2*M))+143 ; 
 if C2(j,1)<C2min
     C2(j,1)=C2min;
 end
 C3(j,1) = 0.155*sin(j*2*pi/(2*M))+0.1655; 
 C4(j,1) = 0;
 C5(j,1) = 0.0000000035*dt;
 C6(j,1) = 0;
 C7(j,1) = 0;
 C8(j,1) = 0.000000205*dt;
 C9(j,1) = 0.00000015*dt;  
 C10(j,1) = 0;
 C11(j,1) = 0;
 C12(j,1) = 0.8;  %mikawa
 C13(j,1) = 25510;
 C14(j,1) = (0.00001592*sin(j*2*pi/(2*M))+0.00002267)*dt/10; 
 C15(j,1) = (0.00001592*sin(j*2*pi/(2*M))+0.00002267)*dt*5/10; 
 C16(j,1) = (0.00001592*sin(j*2*pi/(2*M))+0.00002267)*dt*4/10; 
 C17(j,1) = 0.000000070*dt; 
 C18(j,1) = 0.000000648*dt;
C2(j,i) = (as2(j-1,i+1)*C2(j-1,i+1)+(1-2*as2(j-1,i))*C2(j-1,i)+as2(j-1,i-1)*C2(j-1,i-1)-bs2(1,i+1)*C2(j-1,i+1)+bs2(1,i-1)*C2(j-1,i-1)+phi*(-2*k1*C1(j-1,i)*C2(j-1,i)-k4*C2(j-1,i)*C6(j-1,i)-k6*C2(j-1,i)*C4(j-1,i)-2*k11*C2(j-1,i)*C7(j-1,i)-2*k12*C2(j-1,i)*C9(j-1,i)-7*k13*C2(j-1,i)*C11(j-1,i)-((kOMf*C14(j-1,i)+kOMs*C15(j-1,i)+kOMr*C16(j-1,i))*(C2(j-1,i)/(C2(j-1,i)+Ks1))))*dt)  ;
      C2(j,N+1) = C2(j,N);