【信号处理】基于matlab模拟频谱混叠、频谱泄漏和栅栏现象

% 定义梁的几何参数
L = 10; % 梁的长度（单位：m）
b = 0.3; % 梁的宽度（单位：m）
h = 0.5; % 梁的高度（单位：m）
A = b * h; % 梁的截面面积（单位：m^2）
I = (b * h^3) / 12; % 梁的惯性矩（单位：m^4）

% 定义材料参数
E = 30e9; % 弹性模量（单位：Pa）
fpu = 1700e6; % 预应力钢筋抗拉强度（单位：Pa）
fp0 = 1300e6; % 初始预应力（单位：Pa）

% 定义荷载
P = 50000; % 荷载大小（单位：N）

% 计算预应力钢筋的应力
Aps = 0.001; % 预应力钢筋的截面面积（单位：m^2）
fp = fp0 * (Aps / A); % 预应力钢筋的应力（单位：Pa）

% 计算梁的弯矩和曲率
x = linspace(0, L, 100); % 梁的长度离散点
M = zeros(size(x)); % 弯矩
curvature = zeros(size(x)); % 曲率

for i = 1:length(x)
% 计算弯矩
if x(i) <= L/2
M(i) = -P * x(i) + fp * Aps * (L - x(i));
else
M(i) = -P * (L - x(i)) - fp * Aps * x(i);
end

% 计算曲率
curvature(i) = M(i) / (E * I);

end

% 绘制弯矩曲率图
figure;
subplot(2,1,1);
plot(x, M, ‘b-’, ‘LineWidth’, 2);
xlabel(‘位置（m）’);
ylabel(‘弯矩（N·m）’);
title(‘弯矩分布’);
grid on;

subplot(2,1,2);
plot(x, curvature, ‘r-’, ‘LineWidth’, 2);
xlabel(‘位置（m）’);
ylabel(‘曲率（1/m）’);
title(‘曲率分布’);
grid on;