Chào mọi người ạ.
Mình đang làm một bài tập là đi tìm chuyển vị u và v của hệ thanh, mình làm bằng M-filen nhưng mà nó chỉ chạy được mấy dòng lệnh đầu thôi chứ phần giải phương trình sau của mình không được mọi người xem giúp mình với ạ
code:
function [u,v]=chuyenvi(anpha)
n= input('so thanh:' )
A= input('dien tich thiet dien cac thanh:' )
E= input(' do cung cua cac thanh:' )
L= input(' chieu dai cac thanh:')
for i=1:n
anpha(i)=input( 'goc hop boi thanh va phuong ngang la: ')
end
y=0
z=0
for i=1:n
y=y+(sqrt(u.^2+v.^2)*cosd(180-anpha(i)-atand(v./u))*cosd(anpha(i)))
z=z+(sqrt(u.^2+v.^2)*cosd(180-anpha(i)-atand(v./u))*sind(anpha(i)))
end
[u,v]=solve('y=0','z=F/(A*E)');
disp('chuyen vi theo phuong ngang la:')
disp(u)
disp('chuyen vi theo phuong doc la:')
disp(v)
end
Cảm ơn mọi người!
Mình đang làm một bài tập là đi tìm chuyển vị u và v của hệ thanh, mình làm bằng M-filen nhưng mà nó chỉ chạy được mấy dòng lệnh đầu thôi chứ phần giải phương trình sau của mình không được mọi người xem giúp mình với ạ
code:
function [u,v]=chuyenvi(anpha)
n= input('so thanh:' )
A= input('dien tich thiet dien cac thanh:' )
E= input(' do cung cua cac thanh:' )
L= input(' chieu dai cac thanh:')
for i=1:n
anpha(i)=input( 'goc hop boi thanh va phuong ngang la: ')
end
y=0
z=0
for i=1:n
y=y+(sqrt(u.^2+v.^2)*cosd(180-anpha(i)-atand(v./u))*cosd(anpha(i)))
z=z+(sqrt(u.^2+v.^2)*cosd(180-anpha(i)-atand(v./u))*sind(anpha(i)))
end
[u,v]=solve('y=0','z=F/(A*E)');
disp('chuyen vi theo phuong ngang la:')
disp(u)
disp('chuyen vi theo phuong doc la:')
disp(v)
end
Cảm ơn mọi người!
Comment