邢哈哈哈 发表于 2018-11-15 21:43

关于数值解画频响曲线的matlab程序，希望和大家一起讨论下

function dx=funpin(t,x,w)
P=100;L=0.1;H=3e-3;miu1=18e3;miu2=42e3;
J1=110; J2=55;KEXI=3.9825e-11;ROU=960;
dx=zeros(3,1);
dx(1)=x(2);
dx(2)=((3*P)/(4*(L^3)*H*ROU))+((3*KEXI*((10000+10000*cos(2*pi*w*t))^2)*x(1)^3)/(4*H^2*L^2*ROU))-...
    (1/(L^2*ROU))*((3*miu1*(x(1)-x(1)^-5))/(1-((2*x(1)^2+(x(1)^-4)-3)/J1))+...
    (3*miu2*(x(1)*(x(3)^-2)-(x(1)^-5)*x(3)^4))/(1-((2*x(1)^2*x(3)^-2+(x(1)^-4)*x(3)^4-3)/J2)));
dx(3)=((1/6000)*x(3))*(miu2*((x(1)^2)*(x(3)^-2)-(x(1)^-4)*(x(3)^4)))/(1-((2*x(1)^2*(x(3)^-2)+(x(1)^-4)*(x(3)^4)-3)/J2));
      


邢哈哈哈 发表于 2018-11-15 21:44

clc;
clear all;
close all;
tspan=(0:0.0001:10);

x0=;

%options = odeset('RelTol',1e-9,'AbsTol',1e-11);

B=[];
A=[];
for w=0:0.1:50   
   options = odeset('RelTol',1e-9,'AbsTol',1e-11);
B=;   
=ode45(@(t,x)funpin(t,x,w),tspan,x0,options);
FU=max(x(end-50000:end,1))-min(x(end-50000:end,1));
A=;
end
%B
%A

plot(B,A)

邢哈哈哈 发表于 2018-11-15 21:47

我就想知道，为啥得到的幅频曲线的幅值在初始阶段这么大？与文献对不上。然而在t >10以后就基本可以对上。困扰我好久了。。。。。

邢哈哈哈 发表于 2018-11-15 21:51

E:\

邢哈哈哈 发表于 2018-11-15 21:52

E:\图

邢哈哈哈 发表于 2018-11-15 21:54

这是文献的图，最上方的那条线就是我要得到的线

邢哈哈哈 发表于 2018-11-16 11:05

大神在哪里啊？求助啊

邢哈哈哈 发表于 2018-11-27 09:32

大神们来啊

xcy3508 发表于 2018-11-27 19:20

楼主，你问题解决了吗？我最近也遇到类似的问题，可否一块交流一下？

邢哈哈哈 发表于 2018-12-3 16:19

xcy3508 发表于 2018-11-27 19:20
楼主，你问题解决了吗？我最近也遇到类似的问题，可否一块交流一下？

没有啊，你遇到的问题是什么？

xcy3508 发表于 2018-12-3 16:24

邢哈哈哈 发表于 2018-12-3 16:19
没有啊，你遇到的问题是什么？

我做齿轮动力学的扫频分析，也是用ode45算频响曲线，也是对不上文献上的数据

邢哈哈哈 发表于 2018-12-3 16:44

xcy3508 发表于 2018-12-3 16:24
我做齿轮动力学的扫频分析，也是用ode45算频响曲线，也是对不上文献上的数据

你的幅频曲线是都对不上还是一部分对不上

邢哈哈哈 发表于 2018-12-3 16:44

xcy3508 发表于 2018-12-3 16:24
我做齿轮动力学的扫频分析，也是用ode45算频响曲线，也是对不上文献上的数据

你的幅频曲线是都对不上还是一部分对不上

xcy3508 发表于 2018-12-3 22:15

邢哈哈哈 发表于 2018-12-3 16:44
你的幅频曲线是都对不上还是一部分对不上

我是都对不上，很奇怪，我是非线性系统，有多值和跳跃现象

邢哈哈哈 发表于 2018-12-5 09:29

xcy3508 发表于 2018-12-3 22:15
我是都对不上，很奇怪，我是非线性系统，有多值和跳跃现象

哦，你是用什么方法画的幅频？也是数值解法？
页: [1] 2
查看完整版本: 关于数值解画频响曲线的matlab程序，希望和大家一起讨论下