Matlab7.0 运行问题

来源:百度知道 编辑:UC知道 时间:2024/07/01 04:09:04
刚学matlab

%该程序用于计算SDOF体系的加权假设位移迭代分析
%使用格式为displacement(M,Cz,K,P,u(0),v(0))

syms M Cz K P t i n

i=0
T=10,t=0.02,n=T/t
a=cell(n,1)
b=cell(n,1)
c=cell(n,1)
d=cell(n,1)
e=cell(n,1)
f=cell(n,1)
u=cell(n,1)
v=cell(n,1)
ac=cell(n,1)

u{0,1}=0
v{0,1}=5

M=40,K=4000,C=0,P=0

while i<=n-1
solve('f{i,1}-u{i,1}=0'...
'e{i,1}-v{i,1}=0'...
'd{i,1}-(P-Cz*v{i,1}-K*u{i,1})/(2*M)'...
'-int(Cz*v{i,1},t,0,T)-int(K*u{i,1},t,0,T)+int(P,t,0,T)=M*(v{i+1,1}-v{i,1})'...
'int(Cz*v{i,1}*v{i,1},t,0,T)+int(k*u{i,1}*v{i,1},t,o,T)-int(P*v{i,1},t,0,T)=1/2*M*(v{i+1,1}^2-v{i,1}^2)'...
'a{i,1},b{i,1},c{i,1},d{i,1},e{i,1},f{i,1}')

u{i+1,1}=a{i,1}*t^5+b{i,1}*t^4+c{i,1}*t^3+d{i,1}*t^2+e{i,1}*t+f{i,1}
v{i+1,1}=diff(u{i,1},t)
ac{i+1,1}=diff(u{i,

displacement这个M文件没有
这么简单也不知道?
还写这么好的程序代码,不知道是不是抄的?

不太清楚具体函数的意思,不过你把solve函数后的“... ”全部删除试试。