level 1
12ss40
楼主
解非线性方程组问题,代码如下:
Parameter c1,c2,c3,c4,t1,t2;
Constant k1=2.76e3,F0=25,V=1e-6,k2=((pi*206e9*5e-3)/(4*(1-0.3^2))),w=2*pi*34000,a=sqrt((k1+k2)/(4.8/9.8)),b=(k2*V)/(k1+k2-(4.8/9.8)*w^2),c=(F0+4.8)/(k1+k2),d=sqrt(k1/(4.8/9.8)),M=(F0+4.8)/k1,n=(k2*V*w)/(k1+k2-(4.8/9.8)*w^2);
Function c1*cos(a*t1)+c2*sin(a*t1)+b*sin(w*t1)-c-V*sin(w*t1)=0;
c3*cos(d*t1)+c4*sin(d*t1)-M-V*sin(w*t1)=0;
-c1*a*sin(a*t1)+c2*a*cos(a*t1)+n*cos(w*t1)+c3*d*sin(d*t1)-c4*d*cos(d*t1)=0;
c1*cos(a*t2)+c2*sin(a*t2)+b*sin(w*t2)-c-V*sin(w*t2)=0;
c3*cos(d*t2)+c4*sin(d*t2)-M-V*sin(w*t2)=0;
-c1*a*sin(a*t2)+c2*a*cos(a*t2)+n*cos(w*t2)+c3*d*sin(d*t2)-c4*d*cos(d*t2)=0;
牛顿法求解后t1与t2的值一样,但应该是不一样的,限定t2的值范围后多次求解的值不一样,求助如何能求出最优解,初接触1stopt,很多地方不懂,代码也是照搬的,
2020年03月16日 07点03分
1
Parameter c1,c2,c3,c4,t1,t2;
Constant k1=2.76e3,F0=25,V=1e-6,k2=((pi*206e9*5e-3)/(4*(1-0.3^2))),w=2*pi*34000,a=sqrt((k1+k2)/(4.8/9.8)),b=(k2*V)/(k1+k2-(4.8/9.8)*w^2),c=(F0+4.8)/(k1+k2),d=sqrt(k1/(4.8/9.8)),M=(F0+4.8)/k1,n=(k2*V*w)/(k1+k2-(4.8/9.8)*w^2);
Function c1*cos(a*t1)+c2*sin(a*t1)+b*sin(w*t1)-c-V*sin(w*t1)=0;
c3*cos(d*t1)+c4*sin(d*t1)-M-V*sin(w*t1)=0;
-c1*a*sin(a*t1)+c2*a*cos(a*t1)+n*cos(w*t1)+c3*d*sin(d*t1)-c4*d*cos(d*t1)=0;
c1*cos(a*t2)+c2*sin(a*t2)+b*sin(w*t2)-c-V*sin(w*t2)=0;
c3*cos(d*t2)+c4*sin(d*t2)-M-V*sin(w*t2)=0;
-c1*a*sin(a*t2)+c2*a*cos(a*t2)+n*cos(w*t2)+c3*d*sin(d*t2)-c4*d*cos(d*t2)=0;
牛顿法求解后t1与t2的值一样,但应该是不一样的,限定t2的值范围后多次求解的值不一样,求助如何能求出最优解,初接触1stopt,很多地方不懂,代码也是照搬的,