Bài tập matlab

[darkbaron92]

Chào mọi người ạ.
Mình có làm một bài matlab mà không biết nó sai ở chỗ nào nữa ạ mong mọi người chỉ giúp mình với để mình hoàn thiện bài.
Bài của mình đây ạ:
function chuyenvi(u,v)
syms u v
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 voi phuong ngang cua thanh thu i')
end
syms y=0
syms z=0
for i=1:n
y=y+sqrt(u^2+v^2)*cos(180-anpha(i)-atan(v/u))*cos(anpha(i))
z=z+sqrt(u^2+v^2)*cos(180-anpha(i)-atan(v/u))*sin(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


Mình xin cảm ơn!

 

[mduc]

Ðề: Bài tập matlab

Hàm cos(góc tính theo radian). Ở trên bạn dùng cos(180-anpha(i)-atan(v/u)) là sai. Muốn tính lượng giác góc đơn vị độ dùng hàm cosd

 

[darkbaron92]

Ðề: Bài tập matlab

oh thanks bạn, nhưng mà nó vẫn sai bạn ạ, nó ko chạy ra kết quả đó mình nghĩ nó lỗi code nếu mà lỗi lượng giác thì nó chỉ ra sai thôi còn vẫn ra :1:

 

[mori]

Ðề: Bài tập matlab

function chuyenvi(u,v)
syms u v
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 voi phuong ngang cua thanh thu i')
end
syms y=0
syms z=0
for i=1:n
y=y+sqrt(u^2+v^2)*cos(180-anpha(i)-atan(v/u))*cos(anpha(i))
z=z+sqrt(u^2+v^2)*cos(180-anpha(i)-atan(v/u))*sin(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

Mình chưa hiểu sao chỗ bạn định nghĩa syms sao bạn lại còn ghi thêm "=0" với mục đích gì vậy???
và bạn muốn hiện công thức cho u với v thì bạn không nên cho vào trong nháy đơn nữa ' '

 

[darkbaron92]

Ðề: Bài tập matlab

Bạn ơi cho mình hỏi mình đã làm được cho nó chạy một số cái nhưng một số trường hợp nó lại ko chạy mà nó ra empty ví dụ như là:>> [u,v]=chuyenvi(3,[pi/3 pi/2 135*pi/180])Warning: Explicit solution could not be found.
> In solve at 81
In chuyenvi at 15

u =

[ empty sym ]



v =


[]
ko biết nó là lỗi do cái gì ạ? đây là code của mình đã làm ngắn đi
function [u,v]=chuyenvi(n,anpha)
L=1;%m
E=200*10^9;%N/m2
A=10^(-4);%m2
F=10000;%N
z=0;
y=0;
syms u v
for i=1:n
y=y+(cos(pi-anpha(i)-atan(v/u))*cos(anpha(i)));
z=z+((sqrt(u^2+v^2))*cos(pi-anpha(i)-atan(v/u))*sin(anpha(i)));
end
y1=y;
z1=z-(F*L)/(A*E);
[u,v]=solve(y1,z1);
end