数值计算杜哈梅积分

ion_ura

于 2023-05-11 00:00:00 发布

阅读量1.3k

点赞数 2

文章标签：科技

本文链接：https://blog.csdn.net/ion_ura/article/details/130586767

版权

结构动力学的大作业。记录下。

给出如下系统：

其中,F(t)为一个复杂激励。该系统的运动方程为：

选取如下F（t）随时间变化的规律。

杜哈梅积分是h（t）和F(t)的卷积。h(t)的表达形式为：

为便于计算，可以将杜哈梅积分写成如下形式：

如果用矩阵来表达的话，可以写作：

其中，H矩阵由h矩阵倒序得到。

因此可以写出matlab代码：

t=linspace(0,100,1001);
x=zeros(1,1001);
F_t=6*sin(0.5*t)+5.5*sin(3.6*t/2)+7.6*sin(-2.3*t/2)+4.4*sin(-1.1*t);
h=exp(-0.5*10*t).*sin(8.66*t)/8.66;
for i =1:1001
    x(1,i)=sum(flip(h(1:i)).*F_t(1:i));
end
plot(t,x(1,:),t,F_t(1,:)/5,'linewidth',1.2);
legend('X=X(t)单位：m','F=F(t),单位：5N');

最终结果：