结果表明，搜索每个负载值都要花费20秒以上的时间，即使如此，搜索似乎也会失败。
请注意，我以0.002结束，这是我在不超时的情况下对这个在线版本所能做到的最大值。
抽搐
信谊ω
一个= 3.666;
k = 0.467;
s = 100;
r = 1.4;
ef = 0.08483;
N = (* k) / (1 + (a * k));
Eloadvals = 0:0.001:0.002;
为Eloadidx = 1:长度(eloadvals)
Eload = eloadvals(eloadidx);
c1 =(其它)* (((r *) / 2) + (r / 2) 1)——(r1);
c2 = (N * (r / (2 *)) + (r / 2) 1);
Ac2 = ((r *) *(其它))/ 2;
Ah1 = (r * N) / (2 *);
= 8 *(其它)* c2 * sin (pi-atan (s /((其它)*ω)))+ 8 * N * c1 *罪(2 * (pi-atan (1 / (N *ω))))+ (32 * N * ef + 32 * * eload);
B = 16 *(其它)*(1 *(((其它)*ω)/ (2 * (s)) * (((Ah1 * c1 * sin (pi-atan (1 / (N *ω)))* cos (pi-atan (1 / (N *ω))))+ ((((N * s) /((其它)*ω))+ (N * (r / 2 - 1))) * c2 * cos((每股((((其它)* ((r / 2) 1)) / (s)) *ω))+ (pi-atan (s /((其它)*ω))))* cos (pi-atan (s /((其它)*ω)))))))* c2 *罪(2 * (pi-atan (s /((其它)*ω))))+ (16 * (((N *ω)/ 2)*((((((其它)/ (N *ω))+((其它)* (((r / 2) 1) - (r1)))) * c1 * cos((每股((((其它)* ((r / 2)(1)——(r1)) /(其它))* (N *ω)))+ (pi-atan (1 / (N *ω))))* cos(每股((((其它)* ((r / 2)(1)——(r1)) /(其它))* (N *ω))))+ (Ac2 * c2 * sin(π谭(s /((其它)*ω)))* cos (pi-atan (s /((其它)*ω)))))))+ 2 * c1) * N * c1 *罪(2 * (pi-atan (1 / (N *ω))))+ N * c1 * c1 *罪(4 * (pi-atan (1 / (N *ω))));
thissol = vpasolve(A==B);
如果长度(thissol) >= 1
Sol (eloadidx) = thissol(1);
其他的
Sol (eloadidx) = nan;
结束
结束
toc
情节(eloadvals,索尔)
由于vpasolve()失败，在情节中看不到任何东西。
