我想用vpasolve解决一些方程式。
ct(1,j)==r1+r5+r6+r8,
ct(2,j)==r2+r6+r7,
ct(3,j)==r3+r8+r9+r10+r11+r12,
r5==c(4,j)*r1*10^1.98,
r6==c(4,j)*r1*r2*10^1.08,
cct(4)/v(1,j)-c(4,j)-r5-r6-r8-r9-2*r10==c(4,j)*r2*10^2.25,
r8==c(4,j)*r1*r3*10^2.48,
r9==c(4,j)*r3*10^4.04,
r10==c(4,j)^2*r3*10^5.38,
r11==1.9*10^(-3)*r3/r1,
r12==2.5*10^(-3)*r3/r1^2
到目前为止，我曾经使用index=find(r1>=0)...来消除负面影响（这是不相关的）但现在有些r没有任何正面影响。
我想这样做：
if index=find(r1>=0) is empty
r1=0
...直到我找到所有否定并消除它们。
事情是我的方程式会自行生成（就像将来一样，但现在我不知道该怎么做），我需要一种方法来自动归零r并消除它们方程。
之前我曾问过，但我不确定这是同样的问题...... matlab making a custom vpasolve
谢谢你， 丹。