建设博客网站步骤,做跨境电商需要哪些条件,湘潭seo,专业建站团队附件给出的是齿轮动力学振动模型和微分方程其中x为要求的位移响应#xff0c;Zeta为齿轮副相对阻尼系数#xff0c;Zeta 0.07#xff0c;m为等效质量m3#xff1b;Fn为法向载荷Fn100N。k(t)为齿轮副的时变啮合刚度对k(t)我做了关于时间t的傅里叶变换#xff1a;k750527161…附件给出的是齿轮动力学振动模型和微分方程其中x为要求的位移响应Zeta为齿轮副相对阻尼系数Zeta 0.07m为等效质量m3Fn为法向载荷Fn100N。k(t)为齿轮副的时变啮合刚度对k(t)我做了关于时间t的傅里叶变换k750527161.0161.2466e008*sin(2*pi*f0*t-0.86045)8.3314e007*sin(4*pi*f0*t-0.2292)4.9789e007*sin(6*pi*f0*t0.41062)1.7784e007*sin(8*pi*f0*t1.0182)7.633e006*sin(10*pi*f0*t-1.2305);这个方程是变系数微分方程其中就k(t)一个时变系数k(t)该如何表示是否可把k(t)作为t的函数直接代入方程下面是我编的M文件和求解程序各位大神帮我看看处理的对不对谢谢M文件function Yd fun7( t,x,k,zeta,m,Fn)% FUN7.M 时变微分方程例子Fn100;zeta0.07;m3;n13000;z145;Tz0.00044444;f01/Tz;k750527161.0161.2466e008*sin(2*pi*f0*t-0.86045)8.3314e007*sin(4*pi*f0*t-0.2292)4.9789e007*sin(6*pi*f0*t0.41062)1.7784e007*sin(8*pi*f0*t1.0182)7.633e006*sin(10*pi*f0*t-1.2305);dxzeros(2,1);dx(1)x(2);Yd[x(2);Fn/m-2*zeta*sqrt(k/m)*x(2)-k/m*x(1)];求解程序clearm3;Tz0.00044444;f01/Tz;tspanlinspace(0,0.00044444,120);Y0[1.5040e-007;0];[t,XX]ode45(fun7,tspan,Y0);XXX(:,1);YXX(:,2);plot(t,X),xlabel(经历时间/(s)),ylabel(位移x),figure;plot(t,Y),xlabel(经历时间/(s)),ylabel(速度v),figure;plot(X,Y)
