杜哈梅积分的matlab程序(共2页).doc
《杜哈梅积分的matlab程序(共2页).doc》由会员分享,可在线阅读,更多相关《杜哈梅积分的matlab程序(共2页).doc(2页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
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*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
2、)=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 for i=1:aa/bb-1 v1(i)=(y(i+1)-y(i)/bb;%计算速度endfor i=1:(aa/bb-2) a(i)=(v1(i+1)-v1(i)/bb
3、;%计算加速度endhold on%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)%画加速度图hold off legend(fontsize9fontname黑体 位移,fontsize9fontname黑体 速度,fontsize9fontname黑体 加速度) %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) 专心-专注-专业
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 杜哈梅 积分 matlab 程序
限制150内