试了很多次
希望能够做迭代法
目标是想要让得到的y值和输入的x值最为接近
但是却一直出错,和error
不知道是哪里出问题了ＱＱ
希望板上的大大可以帮忙看一下ＱＱＱＱＱ
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