杜哈梅积分的matlab程序(共2页).doc

精选优质文档-倾情为你奉上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*tpi);%荷载函数 %修改参数完毕，接下来的就不用管了。 %y2=3/16*(1/(1-0.25)*(sin(theta*t1)-0.5*sin(w*t1);for i=1:(aa/bb) for j=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