专门查大学的网站,wordpress4.9.8升级包,推广自己的店铺推广语,网站建设服务器是什么意思附件给出的是齿轮动力学振动模型和微分方程其中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)
