'struct' 类型的操作数不支持运算符 '>'。 出错 fudongfudong (第 45 行) ee = double(sol(sol>0)); % 找到正根还是错的啊
时间: 2023-12-03 16:04:43 浏览: 63
这个错误提示是因为在 `sol` 变量中,可能没有找到正的根,导致 `sol>0` 的结果是一个 `struct` 类型的空数组,不支持使用运算符 ">" 进行比较。可以在使用 `sol` 前,先判断一下 `sol` 是否为空数组,如果是的话,可以给 `ee` 赋一个默认值。
以下是修改后的代码:
```matlab
clc
clear
% 数值法
% 初值
tmin = 0;
tmax = 100;
% 精度
d_doc = 1;
doc = (tmax-tmin)/d_doc;
% 参数直接在后面改
Pf = 10;
m = 700;
ii = 0.03;
% 记得改
i0 = 0.02;
nx = 45;
r = 0.7*0.01;
E = 1;
theta = 0.1;
d = -0.01;
gamma = 1;
kc = 20;
aerfa = 0.7;
lamuda = 0.8;
fai = 10;
beita = 1;
w1 = 2;
w2 = 1;
n = 0.13;
P0 = 25;
huibig = 25;
iworld = 0.025;
miu = 33600;
syms ee dp p
P1 = -m*beita*(i0+d)*huibig*Pf/(((-fai*theta-(w1-w2)*E-log(n)+i0*beita+d*beita)...
*(kc-huibig)*Pf*((-fai*theta-(w1-w2)*E-log(n)+i0*beita)/beita/aerfa)^(aerfa/(aerfa-1)))-beita*m*(i0+d)*E)
eqn = miu*(-fai*theta-(w1-w2)*ee-log(n)-iworld*beita)/beita- m*(i0+d)...
*(E*p-kc*Pf)*beita/p/(-fai*theta-(w1-w2)*E-log(n)+d*beita)/(kc-huibig)/Pf==0;
sol = solve(eqn, E);
if isempty(sol)
ee = 1; % 如果没有找到正的根,则默认为1
disp('未找到正根,已将 ee 赋默认值为 1');
else
ee = double(sol(sol>0)); % 找到正根
end
disp(ee);
T = linspace(tmin,tmax,doc);
dt = T(2)-T(1);
for i = 1:doc
result_p(i) = P0;
p = P0;
eqn = (-fai*theta-(w1-w2)*ee-log(n))/beita+i0-dp/p...
-aerfa*(beita*m*(ee*p-huibig*Pf)*(i0+d)/p/(-fai*theta-(w1-w2)*ee-log(n)+i0*beita+d*beita)...
/(kc-huibig)/Pf)^((aerfa-1)/aerfa)==0;
temp_dp = solve(eqn,dp);
temp_dp = double(min(real(temp_dp)));
dp1(i) = temp_dp;
P0 = P0 + temp_dp*dt;
disp(["计算中...",string(i/doc*100)," %"]);
end
figure
plot(T,result_p)
xlabel("t")
ylabel("p")
figure
plot(T,dp1);
xlabel("t")
ylabel("dp")
dp_p = dp1./result_p;
figure;
plot(T,dp_p)
xlabel("t")
ylabel("dp/p")
```
希望能够帮到你!