[matlab一些代码]好像是关于拟合曲线的＜忘记了＞

clear all
i1 = 1.841;
i2 = 5.331;
i3 = 8.536;
i4 = 11.706;
i5 = 14.864;
i6 = 18.016;
i7 = 21.164;
i8 = 24.311;
i9 = 27.457;
i10 = 30.602;
r = 0.14;
p2 = 1000;
p1 = 862;% 油的密度
detal1 = 0.0344;  %线性阻尼率
%list_w = [0.5672,0.9922,1.2556,1.4704,1.6569,1.8241,1.9771,2.119,2.2519,2.3774]; % 2号油水界面

%list_w = [0.5395,0.9435,1.194,1.3982,1.5756,1.7346,1.88,2.015,2.1414,2.2607]; % 15号油i油水界面

list_w = [0.479,0.8374,1.0597,1.2409,1.3983,1.5395,1.6685,1.7883,1.9005,2.0064]; % 68号油油水界面

%list_w = [1.8063,3.0761,3.8924,4.5582,5.1364,5.6548,6.129,6.5689,6.981,7.37];%自由液面

w1 = 2*pi*list_w(1);   %1-10阶固有频

w2 = 2*pi*list_w(2);
w3 = 2*pi*list_w(3);
w4 = 2*pi*list_w(4);
w5 = 2*pi*list_w(5);
w6 = 2*pi*list_w(6);
w7 = 2*pi*list_w(7);
w8 = 2*pi*list_w(8);
w9 = 2*pi*list_w(9);
w10 = 2*pi*list_w(10);
a = 0.056;  %激励幅值
j = 1;
for i = 0.9:0.01:1.1
t = 0:0.01:30;   
%fanwei = 0.5:0.01:3.0
w = 2*pi*list_w(1)*i;  %激励频率
% d1 = ((p2 - p1)*a*r*w*w*w1^2*(w1^2 - w^2))/((p2 - p1)*9.81*(i1^2 - 1)*((w1^2-w^2)^2+4*w^2*detal1^2));
% d2 = ((p2 - p1)*a*r*w*w*w2^2*(w2^2 - w^2))/((p2 - p1)*9.81*(i2^2 - 1)*((w2^2-w^2)^2+4*w^2*detal1^2)); 
% d3 = ((p2 - p1)*a*r*w*w*w3^2*(w3^2 - w^2))/((p2 - p1)*9.81*(i3^2 - 1)*((w3^2-w^2)^2+4*w^2*detal1^2));
% d4 = ((p2 - p1)*a*r*w*w*w4^2*(w4^2 - w^2))/((p2 - p1)*9.81*(i4^2 - 1)*((w4^2-w^2)^2+4*w^2*detal1^2));
% d5 = ((p2 - p1)*a*r*w*w*w5^2*(w5^2 - w^2))/((p2 - p1)*9.81*(i5^2 - 1)*((w5^2-w^2)^2+4*w^2*detal1^2));
% d6 = ((p2 - p1)*a*r*w*w*w6^2*(w6^2 - w^2))/((p2 - p1)*9.81*(i6^2 - 1)*((w6^2-w^2)^2+4*w^2*detal1^2));
% d7 = ((p2 - p1)*a*r*w*w*w7^2*(w7^2 - w^2))/((p2 - p1)*9.81*(i7^2 - 1)*((w7^2-w^2)^2+4*w^2*detal1^2));
% d8 = ((p2 - p1)*a*r*w*w*w8^2*(w8^2 - w^2))/((p2 - p1)*9.81*(i8^2 - 1)*((w8^2-w^2)^2+4*w^2*detal1^2));
% d9 = ((p2 - p1)*a*r*w*w*w9^2*(w9^2 - w^2))/((p2 - p1)*9.81*(i9^2 - 1)*((w9^2-w^2)^2+4*w^2*detal1^2));
% d10 = ((p2 - p1)*a*r*w*w*w10^2*(w10^2 - w^2))/((p2 - p1)*9.81*(i10^2 - 1)*((w10^2-w^2)^2+4*w^2*detal1^2));
d1 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i1^2 - 1))*(w1^2*(w1^2 - w^2)*cos(w*t)/((w1^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w1^2*sin(w*t)/((w1^2-w^2)^2+4*detal1^2*w^2));
d2 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i2^2 - 1))*(w2^2*(w2^2 - w^2)*cos(w*t)/((w2^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w2^2*sin(w*t)/((w2^2-w^2)^2+4*detal1^2*w^2)); 
d3 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i3^2 - 1))*(w3^2*(w3^2 - w^2)*cos(w*t)/((w3^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w3^2*sin(w*t)/((w3^2-w^2)^2+4*detal1^2*w^2));
d4 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i4^2 - 1))*(w4^2*(w4^2 - w^2)*cos(w*t)/((w4^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w4^2*sin(w*t)/((w4^2-w^2)^2+4*detal1^2*w^2));
d5 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i5^2 - 1))*(w5^2*(w5^2 - w^2)*cos(w*t)/((w5^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w5^2*sin(w*t)/((w5^2-w^2)^2+4*detal1^2*w^2));
d6 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i6^2 - 1))*(w6^2*(w6^2 - w^2)*cos(w*t)/((w6^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w6^2*sin(w*t)/((w1^6-w^2)^2+4*detal1^2*w^2));
d7 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i7^2 - 1))*(w7^2*(w7^2 - w^2)*cos(w*t)/((w7^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w7^2*sin(w*t)/((w2^7-w^2)^2+4*detal1^2*w^2)); 
d8 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i8^2 - 1))*(w8^2*(w8^2 - w^2)*cos(w*t)/((w8^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w8^2*sin(w*t)/((w3^8-w^2)^2+4*detal1^2*w^2));
d9 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i9^2 - 1))*(w9^2*(w9^2 - w^2)*cos(w*t)/((w9^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w9^2*sin(w*t)/((w4^9-w^2)^2+4*detal1^2*w^2));
d10 = ((p2 - p1)*a*2*r*w*w)/((p2 - p1)*9.81*(i10^2 - 1))*(w10^2*(w10^2 - w^2)*cos(w*t)/((w10^2-w^2)^2+4*w^2*detal1^2)+2*detal1*w*w10^2*sin(w*t)/((w10^2-w^2)^2+4*detal1^2*w^2));
d = d1 + d2 + d3 + d4 + d5 + d6 + d7 + d8 + d9 + d10;  
d_max = max(d);
% tu = d*cos(w*t)
% plot(t,tu)
f(j) = abs(d_max);
j = j+1;
end
x111 = 0.9:0.01:1.1;
plot(f)
plot(x111,f)
F = f'