轴承静力学分析的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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(4)))/41 - 1) + (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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(4)))/41 - 1))/(cos(x(4) - 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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(3)))/41 + 1)))^2)/42501298345826806923264;

f(5)=100000000*x(1) + 100000000*x(2) - (332087940869622973*pi^2*(2*sin(x(3)) - (3883*x(1)*sin(x(4))*(((4*cos(x(3)))/41 + 1)/((4*cos(x(4)))/41 - 1) - 1))/(3676*x(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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(4)))/41 - 1) + ((4*cos(x(4)))/41 - 1)/cos(x(4) - 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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41))))*(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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(4)))/41 - 1) + (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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(4)))/41 - 1))/(cos(x(4) - 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))/(2*cos(x(3)) + 3883/37679*x(1)/x(2)*(cos(x(4) - x(3)) + (4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1)) - 3883/3676*x(1)/x(2)*sin(x(4))*((4/41*cos(x(3)) + 1)/(4/41*cos(x(4)) - 1) - 1) + 8/41)))*((4*cos(x(3)))/41 + 1))))/(5570730176784211237046059008*((2*sin(x(3)) - (3883*x(1)*sin(x(4))*(((4*cos(x(3)))/41 + 1)/((4*cos(x(4)))/41 - 1) - 1))/(3676*x(2)))^2/(2*cos(x(3)) + (3883*x(1)*(cos(x(4) - x(3)) + ((4*cos(x(3)))/41 + 1)/((4*cos(x(4)))/41 - 1)))/(37679*x(2)) - (3883*x(1)*sin(x(4))*(((4*cos(x(3)))/41 + 1)/((4*cos(x(4)))/41 - 1) - 1))/(3676*x(2)) + 8/41)^2 + 1)^(1/2)*(2*cos(x(3)) + (3883*x(1)*(cos(x(4) - x(3)) + ((4*cos(x(3)))/41 + 1)/((4*cos(x(4)))/41 - 1)))/(37679*x(2)) - (3883*x(1)*sin(x(4))*(((4*cos(x(3)))/41 + 1)/((4*cos(x(4)))/41 - 1) - 1))/(3676*x(2)) + 8/41));

轴承静力学分析的matlab代码,求球轴承的静力学与拟静力学的matlab程序相关推荐

【数字信号处理】卷积编程实现 ( Matlab 卷积和多项式乘法 conv 函数 | 使用 matlab 代码求卷积并绘图 )
文章目录一.Matlab 卷积和多项式乘法 conv 函数二.使用 matlab 代码求卷积并绘图一.Matlab 卷积和多项式乘法 conv 函数 Matlab 文档地址 : https:// ...
matlab圆周率计算,matlab代码求圆周率的简单算法
说起圆周率的算法很多人都会想起一大堆的无穷级数等各种表达式,但是这样的算法需要比较高的数学推理水品,而且对于很多初学者而言很难理解.程序员不能仅仅是机械的写程序,必须要真正的理解程序中每个代码的意义, ...
arima模型matlab代码_PSTR面板平滑转换模型简介（附Matlab代码分享）
写论文的时候用到的~相关的资料太少了,做一些简单内容和资料的分享.(PSTR模型的Matlab代码分享在最后)本文主要为简单理论和粗暴实操~ 有用的话可以点个赞哟(知乎小白卑微求赞) 嘻嘻下面进入正题 ...
a*算法matlab代码_NSGAII多目标优化算法讲解（附MATLAB代码）
小编今天为大家讲解NSGA-II多目标优化算法,提到多目标优化,大家可能第一个就想到NSGA-II算法,今天小编就带领大家解开NSGA-II的神秘面纱. NSGA-II全称是快速非支配排序遗传算法,这 ...
认知无线电matlab代码详解,认知无线电频谱感知之功率检测matlab代码.docx
认知无线电频谱感知之功率检测matlab代码能量检测仿真实验代码:clear all;clc;n = 5;ps = 1;SNR1 = -5;SNR2 = -8;SNR3 = -10;% Sim_Ti ...
matlab中求分配系数的命令,汽车理论课后习Matlab程序.doc
汽车理论课后习Matlab程序 1.3 确定一轻型货车的动力性能(货车可装用4挡或5挡变速器,任选其中的一种进行整车性能计算): 1)绘制汽车驱动力与行驶阻力平衡图. 2)求汽车最高车速,最大爬坡度 ...
matlab怎么求一个三元一次方程组的解,用MATLAB求解一个带参数的三元一次方程组，求大神指点!...
想求一个方程组,改了很多遍都还是出错,请求大神指点 . 代码如下 i=[1 2 3]; %编号为1的机械臂 r=50; %动平台半径 R=210; 想求一个方程组,改了很多遍都还是出错,请求大神指点 ...
领域平均法matlab代码实验,数字图像处理邻域平均法滤波实验报告matlab实现.doc...
数字图像处理邻域平均法滤波实验报告matlab实现数字图像处理实验报告实验三邻域平均法滤波学号姓名实验三邻域平均法滤波一.实验内容选取噪声较明显的图像,分别采用3*3.5*5.7* ...
matlab中求三维中的多个体积,用matlab计算由下面2个几何体围成的体积： x^2+y^2+z^2=36,((x-4)/5)^2+((y-1)/3)^2+((z-2)/5)^2=1...
答:>> triplequad(@(x,y,z)1*(x.^2+y.^2+z.^2 答:首先建立一个m文件我取的名字叫 syfs0000 function y=syfs0000(x) y ...

轴承静力学分析的matlab代码,求球轴承的静力学与拟静力学的matlab程序

轴承静力学分析的matlab代码,求球轴承的静力学与拟静力学的matlab程序相关推荐

最新文章

热门文章