大侠们帮帮忙，我要达到这样的功效，该怎么办？
#include<iostream>#include<math.h>
using namespace std;
const int N=5690,M=1000,K=1000;
int main()
{
double a[N],n[M],m[K];
int i;
double a1,a2,a3,a4;
for( a1=0.01,i=1; a1<1.01;i++)
{
a[i]=a1;
a1=a1+0.01;
}
for( a2=1.1,i=101; a2<50.1;i++)
{
a[i]=a2;
a2=a2+0.1;
}
for( a3=50.5,i=591; a3<200.5;i++)
{
a[i]=a3;
a3=a3+0.5;
}
for( a4=201,i=891; a4<5001;i++)
{
a[i]=a4;
a4=a4+1;
}
cout<<a[5690]<<endl;
double n1;
for( n1=0.05,i=1; n1<50.05;i++)
{
n[i]=n1;
n1=n1+0.05;
}
cout<<n[1000]<<endl;
double m1;
for( m1=0.01,i=1; m1<10.01;i++)
{
m[i]=m1;
m1=m1+0.01;
}
cout<<m[1000]<<endl;
double x1=0.1,y1=0.4,x2=1,y2=0.399,x3=10,y3=0.38,x4=100,y4=0.24,x5=1000,y5=0.1;
const int L=N*M*K;
double F[L];
for(int l=1;l<L;l++)
for(i=1;i<N;i++)
for(j=1;j<M;j++)
for(k=1;k<K;k++)
//F[l]=(1/(ln(e+(x1/a[i])^n[j]))^m[k]-y1)^2/y1+(1/(ln(e+(x2/a[i])^n[j]))^m[k]-y2)^2/y2+....+(1/(ln(e+(x5/a[i])^n[j]))^m[k]-y5)^2/y5;
//minF[l]
cout<<F[l]<<endl;
//return a[i],n[j],m[k];
cout<<"a="<<a[i]<<";"<<endl;
cout<<"n="<<n[j]<<";"<<endl;
cout<<"m="<<m[k]<<";"<<endl;
return 0;
}
黑体之前的，我有试过，是可以运行的，产生3个数组，也就是本例中函数f(x)=1/[ln(e+(x/a)^n]^m 中三个函数的可能值。 现在有5组数据，（x1,y1);(x2,y2),...(x5,y5)
这5组数据不一定在函数曲线上，我现在要做的是，在所有参数组合中，找出标准平方差最小的那对组合，也就是∑（f(xi)-yi)^2/yi ,最小,此时返回参数a[i];n[j];m[k].
小弟最近正为这个搞得焦头烂额，还请大侠们帮帮忙。 万分感谢。
