用于对矩阵求幂的维度不正确.请检查并确保矩阵为方阵并且幂为标量.要执行按元素矩.pdfVIP

用于对矩阵求幂的维度不正确.请检查并确保矩阵为方阵并且幂为标量.要执行按

元素矩

matlab报错：错误使用^用于对矩阵求幂的维度不正确。请检查并确保矩阵

为方阵并且幂为标量。要执行按元素矩阵求幂，请使用.^

源代码是：

ift1==2%x轴为Db

fori=1:100

Db=0.01+0.0002*i;

symsyx%y=ubvx=cdv

ifrhobrhol

%默认第一个条件

[x,

y]=solve(x==24/129*abs(y-ulv)*196.848*39.3701*Db*rhol/1000/mu,y==

sqrt(2*m*g*(1-rhol/rhob)/x/A/rhol)+ulv);

x=vpa(x(1),5);

y=vpa(y(1),5);

Re=129*abs(y-ulv)*196.848*39.3701*Db*rhol/1000/mu;

ifRe0.1Re1000%第二个

symsyx

[x,

y]=solve(x==24/Re*(1+0.14*Re^0.7),y==sqrt(2*m*g*(1-rhol/rhob)/x/A/r

hol)+ulv);

x=vpa(x(1),5);

y=vpa(y(1),5);

Re=129*abs(y-ulv)*196.848*39.3701*Db*rhol/1000/mu;

elseifRe=1000Re=35000%第三个

symsyx

[x,

y]=solve(x==0.445,y==sqrt(2*m*g*(1-rhol/rhob)/x/A/rhol)+ulv);

x=vpa(x(1),5);

y=vpa(y(1),5);

Re=129*abs(y-ulv)*196.848*39.3701*Db*rhol/1000/mu;

elseifRe=35000Re2.9*10^5%第四个

symsyx

[x,

y]=solve(x==24/Re*(1+0.1806*Re^0.6459)+0.4251/(1+6880.95*Re^(-1))

,y==sqrt(2*m*g*(1-rhol/rhob)/x/A/rhol)+ulv);

x=vpa(x(1),5);

y=vpa(y(1),5);

Re=129*abs(y-ulv)*196.848*39.3701*Db*rhol/1000/mu;

elseifRe=2.9*10^5Re=3.85*10^5%第五个

symsyx

[x,

y]=solve(x==-4.16*10^(-6)*Re+1.67,y==sqrt(2*m*g*(1-rhol/rhob)/x/A/r

hol)+ulv);

x=vpa(x(1),5);

y=vpa(y(1),5);

Re=129*abs(y-ulv)*196.848*39.3701*Db*rhol/1000/mu;