用matlab画超高斯激励下的二自由度车辆非线性系统的幅频曲线

由于没有给出具体的车辆非线性系统模型，因此以下仅提供绘制超高斯激励下二自由度车辆系统的幅频曲线的一般步骤：

1. 构造车辆非线性系统模型，包括质量、阻尼、刚度等参数。

2. 确定超高斯激励的参数，包括时间常数、幅值等。

3. 根据车辆系统模型和超高斯激励参数，利用matlab中的ode45函数求解系统的时域响应。

4. 对时域响应进行傅里叶变换，得到系统的频域响应。

5. 绘制系统的幅频曲线。

下面是一个简单的示例程序：

% 车辆系统参数
m1 = 1000;  % 质量1
m2 = 1500;  % 质量2
k1 = 50000; % 刚度1
k2 = 80000; % 刚度2
c1 = 1000;  % 阻尼1
c2 = 2000;  % 阻尼2

% 超高斯激励参数
t0 = 0.1;   % 时间常数
A = 1;      % 幅值

% 定义系统模型
f = @(t,y) [
    y(2);
    (-k1*y(1)-c1*y(2)+k2*(y(3)-y(1))-c2*(y(4)-y(2)))/m1;
    y(4);
    (k2*(y(1)-y(3))+c2*(y(2)-y(4)))/m2+A*exp(-t^2/t0^2)
];

% 定义初始状态和时间范围
y0 = [0;0;0;0];
tspan = [0,10];

% 求解系统的时域响应
[t,y] = ode45(f,tspan,y0);

% 对时域响应进行傅里叶变换，得到系统的频域响应
Fs = 1000;       % 采样频率
N = length(y);   % 采样点数
f = Fs*(0:(N/2))/N;
Y = fft(y(:,1));
P2 = abs(Y/N);
P1 = P2(1:N/2+1);
P1(2:end-1) = 2*P1(2:end-1);

% 绘制系统的幅频曲线
plot(f,P1);
xlabel('Frequency (Hz)');
ylabel('Amplitude');

需要注意的是，由于车辆系统是非线性的，因此在实际绘制中可能需要进行一些额外的处理，例如使用线性化方法得到近似模型，或者使用更复杂的数值方法求解系统的响应。此外，超高斯激励也是一种比较特殊的激励形式，需要根据具体问题进行选择。