matlab循环选择语句怎么写,请问这个循环语句怎么写?

Q =0.0025375;%('输入体积流率(m^3/s)');

L =609;%('水平井筒长度(m)');

qw=Q/L;%延水平井筒方向单位长度上的注入流量(m^3/s*m);

Cw = 4.2*10^6;%('注入流体定压比热容 (J/m^3*K)');

Cr =2.35*10^6;%('储集层岩石定压比热容 (J/m^3*K)');

pw =1.08;%('注入流体密度(g/m^3)');

pr =2.15;%('储集层岩石密度(g/m^3)');

O =0.2;%('孔隙度(小数)');

KT =2.7778;%('储集层平均导热系数(J/m*K*s)');

T1=100;

Te=100;

re =3;%('油藏径向半径(m)');

rw =0.06858;%('水平井井筒半径(m)');

Aw=pi*rw^2;%水平井筒横截面积(m^2)

x =10;%('水平井筒水平网格数');

Xi=L/x;%(水平井筒水平步长(m);

z=2;%(水平井筒加相邻地层在j方向的网格数)

Zi=1;

ti =5;%('时间增量min');

NR=10;%('径向网格数');

a=(re/rw)^(1/(NR-1));%('径向网格转化为线性网格系数')

A0=(Cw*Aw*pw*Xi*Zi+Cw*pw*qw*Xi*Zi*ti-Cw*pw*qw*Zi*ti+(2*pi*KT*a*Xi*ti)/log(a))/(Cw*Aw*pw*Xi*Zi);%差分方程T(i,j)系数

B0=(Cw*pw*qw*Zi*ti)/(Cw*Aw*pw*Xi*Zi);        %差分方程T(i+1,j)系数

C0=((2*pi*KT*a*Xi*ti)/log(a))/(Cw*Aw*pw*Xi*Zi);%差分方程T(i,j+1)系数

%稀疏矩阵的创建

V0=ones([(x+1)*2,1])*A0;%稀疏矩阵主对角线上的系数,+1是在趾部地层多加一个网格;

V1=ones([(x+1)*2-1,1])*B0;%稀疏矩阵泛1对角线系数,+1是在 趾部地层多加一个网格;

V2=ones([(x+1)*2-10,1])*C0;%稀疏矩阵泛10对角线系数,+1是在趾部地层多加一个网格;

V=diag(V0)+diag(V1,1)+diag(V2,10);%线性方程组稀疏矩阵

jj=1;

T=rand((x+1)*2,1);

D=rand((x+1)*2,1);%D:线性方程组的右端项

D=D((x+1)*2,1);

for ii=1:10;

D(ii,jj)=T1;

end

for ii=11;

D(ii,jj)=Te;

end

for ii=11:21;

D(ii,jj)=(T1+Te)/2;

end

for ii=22;

D(ii,jj)=Te;

end

for i=1

T(i)=V\D;

for i=2:42

T(i)=V\T(i-1);

end

end

运行程序说是

In an assignment  A(I) = B, the number of elements in B and I must be the same.

Error in haole (line 48)

T(i)=V\D;

我这里的V,和D都是列向量,维数也是相同的,不明白为什么会出错?

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值