轴承静力学分析的matlab代码,求球轴承的静力学与拟静力学的matlab程序 - 机械 - 小木虫 - 学术 科研 互动社区...

function solvenonl2_Eqs

clear all;

clc

format short

x0=[0.001 0.001  0.2  0.2];

options = optimset('TolFun',1e-5,'Display','iter');

[x ,favl]=fsolve(@Eqs,x0,options)

function f=Eqs(x)

x=x

%主方程组

f(1)=sin(x(3))*(x(1) + 3/50) + sin(x(4))*(x(2) + 3/50) - 3789537307181861/72057594037927936;

f(2)=cos(x(3))*(x(1) + 3/50) + cos(x(4))*(x(2) + 3/50) - 5566447865454253/36028797018963968 ;

f(3)=100000000*x(2)*cos(x(4)) - 100000000*x(1)*cos(x(3)) - 10000000000*x(1)*sin(x(3)) + 10000000000*x(2)*sin(x(4));

f(4)=10000000000*x(2)*cos(x(4)) - 10000000000*x(1)*cos(x(3)) - 100000000*x(1)*sin(x(3)) - 100000000*x(2)*sin(x(4)) + (6334074799911937*pi^2*(cos(x(3) - atan((2*sin(x(3)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1))/