試了很多次
希望能夠做迭代法
目標是想要讓得到的y值和輸入的x值最為接近
但是卻一直出錯,和error
不知道是哪裡出問題了QQ
希望板上的大大可以幫忙看一下QQQQQ
clear all
close all
for i=1:1:100
x(1)=18;
O=40.6;
TO=tand(O);
SeO=1/sind(O);
b=10.42/200;
h=100/1000;
a=2.81/200;
p=1.64*9.81;
qu=10550;
up=(b*h/(a^2))*cosd(O);
low=(x(i)*(1-sind(O))/2)-sind(O);
aa(i)=acotd(TO+SeO*(1+(up/low(i)))^0.5);
K(i)=tand(2*aa(i)+O);
qt(i)=p/(pi*(b*h*K(i)-a^2));
y(i)=qu/qt(i)
differ(i)=abs(y(i)-x(i))
if differ(i)<0.5
x(i)
else
x(i+1)=x(i)+0.05;
end
end