杜哈梅积分的matlab程序Word文档格式.doc

杜哈梅积分的matlab程序Word文档格式.doc

《杜哈梅积分的matlab程序Word文档格式.doc》由会员分享，可在线阅读，更多相关《杜哈梅积分的matlab程序Word文档格式.doc（2页珍藏版）》请在冰点文库上搜索。

t=bb:

bb:

aa;

t1=t;

%不用改

theta=1;

%输入荷载频率

w=2;

%输入自振频率

m=1;

%输入质量

p0=4;

%输入荷载幅值

%%%%

p0=p0*ones（1,aa/bb）;

p=p0.*sin（theta*t）.*（theta*t<

=pi）+0.*（theta.*t>

pi）;

%荷载函数

%%%%%%%%修改参数完毕，接下来的就不用管了。

%y2=3/16*（1/（1-0.25））*（sin（theta*t1）-0.5*sin（w*t1））;

fori=1:

（aa/bb）

forj=1:

i

canshu1（j）=p（j）/（m*w）*bb*sin（w*（t（i）-t1（j）））;

%杜哈梅积分中的被积函数

%canshu2（j）=p（j）*b*cos（w*t1（j））;

%速度的A

%canshu3（j）=p（j）*b*sin（w*t1（j））;

%速度的B

end

%v（i）=cos（w*t（i））/m*sum（canshu2）+sin（w*t（i））/m*sum（canshu3）;

%%速度值

y（i）=sum（canshu1）;

%%位移值

end

aa/bb-1

v1（i）=（y（i+1）-y（i））/bb;

%计算速度

（aa/bb-2）

a（i）=（v1（i+1）-v1（i））/bb;

%计算加速度

holdon

%%plot（t1,y2）

plot（t,y,'

linewidth'

3）%画位移图

plot（t（1:

aa/bb-1）,v1,'

-r'

'

1.8）%画速度图

%plot（t,v,'

k'

）

aa/bb-2）,a,'

m'

）%画加速度图

holdoff

legend（'

\fontsize{9}\fontname{黑体}位移'

\fontsize{9}\fontname{黑体}速度'

\fontsize{9}\fontname{黑体}加速度'

%i=find（y==（max（y）））;

%disp（'

时间'

%t（i-3:

i+3）'

荷载'

%p（i-3:

位移'

%y（i-3:

速度'

%v1（i-3:

加速度'

%a（i-3: