valarray

alarray定义了一组在两个相同长度和相同类型的valarray类对象之间的数字计算

比如这里的 最小二乘法拟合直线的程序 就会用到alarray:

#include <stdlib.h>

#include <iostream>

#include <valarray>

#include <math.h>
using namespace std;

int main(int argc, char *argv[]) 

{

     cout<<"输入变量个数：";

     int num =0,tmp

     cin>> num;

     valarray<double> data_x(num);

     valarray<double> data_y(num);

     tmp=num;

    for(int x=0;x<num;x++)

   {

         cout<< "输入第"<<x+1<<"项的X:";

         cin>>data_x[x];

     }
for(int y=0;y<num;y++)
{
cout<<"输入第"<<y+1<<"项的Y:";

         cin>> data_y[y];
}

     //系数

     double A = 0.0;

     double B = 0.0;

     double C = 0.0;

     double D = 0.0;

for(int i =0;i<num;i++)

     A = log(data_y[i])+A;

B = data_x.sum ();

C = (data_x * log(data_y)).sum();

D = (data_x * data_x).sum();

double i,j,min=data_y[0],max=data_y[0],e=0;

double k = (A * B - data_x.size () * C) / (B * B  - data_x.size () * D);

double b = (B * C - A * D) / (B * B - data_x.size () * D);

for(tmp=0;tmp<num;tmp++)
{
i=k*data_x[tmp]+b;
j=i-data_y[tmp];

if(j>e)
{
e=j;
}
if(data_y[tmp]>max)
{
max=data_y[tmp];

}

if(data_y[tmp]<min)

{
     min=data_y[tmp];
}
}

     //输出

  
    cout<<"Yi各项求和得A="<<A<<endl;
    cout<<"Xi各项求和得B="<<B<<endl;
    cout<<"Yi*Xi各项求和得C="<<C<<endl;
     cout<<"Xi各项平方再和得D="<<D<<endl;

     cout<<"斜率k:"<<k<<endl;

     cout<<"截距b:"<<b<<endl; 

     cout<<"最佳拟合直线为:\nY="<<exp(b)<<"*Vp^"<<k<<"\n"<<endl;// 最终结果表示为exp(b)*Vp^k

     system("PAUSE");

     return 0;

}