杜哈梅积分的matlab程序.doc

《杜哈梅积分的matlab程序.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、clcclear%输入数据，开始可以改参数了。aa=10;%输入时间长度bb=0.01;%输入精度%%%%%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*

2、sin(w*t1));fori=1:(aa/bb)forj=1:icanshu1(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));%速度的Bend%v(i)=cos(w*t(i))/m*sum(canshu2)+sin(w*t(i))/m*sum(canshu3);%%速度值y(i)=sum(canshu1);%%位移值endfori=1:aa/bb-1v1(i)=(y(

3、i+1)-y(i))/bb;%计算速度endfori=1:(aa/bb-2)a(i)=(v1(i+1)-v1(i))/bb;%计算加速度endholdon%%plot(t1,y2)plot(t,y,'linewidth',3)%画位移图plot(t(1:aa/bb-1),v1,'-r','linewidth',1.8)%画速度图%plot(t,v,'k')plot(t(1:aa/bb-2),a,'m')%画加速度图holdofflegend('fontsize{9}fontname{黑体}位移','fontsize{9}fontname{黑体}速度

4、','fontsize{9}fontname{黑体}加速度')%i=find(y==(max(y)));%disp('时间')%t(i-3:i+3)'%disp('荷载')%p(i-3:i+3)'%disp('位移')%y(i-3:i+3)'%disp('速度')%v1(i-3:i+3)'%disp('加速度')%a(i-3:i+3)'