杜哈梅积分的matlab程序.doc
- 文档编号:1653174
- 上传时间:2022-10-23
- 格式:DOC
- 页数:2
- 大小:17KB
杜哈梅积分的matlab程序.doc
《杜哈梅积分的matlab程序.doc》由会员分享,可在线阅读,更多相关《杜哈梅积分的matlab程序.doc(2页珍藏版)》请在冰豆网上搜索。
clc
clear
%输入数据,开始可以改参数了。
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*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
fori=1:
aa/bb-1
v1(i)=(y(i+1)-y(i))/bb;%计算速度
end
fori=1:
(aa/bb-2)
a(i)=(v1(i+1)-v1(i))/bb;%计算加速度
end
holdon
%%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')%画加速度图
holdoff
legend('\fontsize{9}\fontname{黑体}位移','\fontsize{9}\fontname{黑体}速度','\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)'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 杜哈梅 积分 matlab 程序