untitled9.m

来自「迭代实现非线性函数（电离层延迟模型）的参数拟合问题」· M 代码 · 共 97 行

t=linspace(100,86000,100);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%变量初始值设置
fm=linspace(-50,50,100);
x01=[      %shi jizhi
1
0.1
0.05
0.005
0.002
0.0000005
  
5000
100
20
1

1
-3
-0.5
-0.01073072

];
A11=x01(1,1);%%%%%%%%%%%%%更新x0的值
B1=x01(2,1);
alpha01=x01(3,1);
alpha11=x01(4,1);
alpha21=x01(5,1);
alpha31=x01(6,1);
beta01=x01(7,1);
beta11=x01(8,1);
beta21=x01(9,1);
beta31=x01(10,1);
gama01=x01(11,1);
gama11=x01(12,1);
gama21=x01(13,1);
gama31=x01(14,1);

for i=1:100
   f01(1,i)=It(A11,B1,alpha01,alpha11,alpha21,alpha31,beta01,beta11,beta21,beta31,gama01,gama11,gama21,gama31,t(1,i),30);%%%%%%%%%%计算的初值
end

plot(t,f01,'r')
hold on  

   
for k=j:j
    
 A10=x00(1,k);%%%%%%%%%%%%%更新x0的值
B0=x00(2,k);
alpha00=x00(3,k);
alpha10=x00(4,k);
alpha20=x00(5,k);
alpha30=x00(6,k);
beta00=x00(7,k);
beta10=x00(8,k);
beta20=x00(9,k);
beta30=x00(10,k);
gama00=x00(11,k);
gama10=x00(12,k);
gama20=x00(13,k);
gama30=x00(14,k);
    
    
 
for i=1:100;
f02(1,i)=It(A10,B0,alpha00,alpha10,alpha20,alpha30,beta00,beta10,beta20,beta30,gama00,gama10,gama20,gama30,t(1,i),30);%%%%%%%%%%计算的初值
%f(1,i)=f0(i,1);
end

%plot(t,f01,'r',t,f02,'b')
plot(t,f02)
axis([0,86400,0,10])
hold on
end