构造函数中向量的初始化,该怎么处理

构造函数中向量的初始化
请教：
我编写了一个类:
1、可是在构造函数中给向量赋初值时总出现错误,如何给这个向量初始化?
2、在复数和实数运算中，总出现类型不匹配，为什么？
程序如下：    
//
#include   <iostream>
#include   <cmath>
#include   <complex>
#include   <vector>

using   namespace   std;
typedef   double   dbl;
typedef   complex <double>   cdbl;
class   Nlayer
{
        public:
                Nlayer(unsigned   i_n,
                              vector <dbl> &   i_p,
                              vector <dbl> &   i_h);
                ~Nlayer();
                //计算单个频率下地表视电阻率  
                void   caculate_ps(dbl   w,                           //单个角频率  
                                                  dbl   ps);                       //视电阻率
                                      //         vector <cdbl> &   Zn);   //层阻抗  
                //计算一系列频率下地表视电阻率
                void   caculate_ps(vector <dbl> &   ws,       //w：频率向量  
                                                  vector <dbl> &   pss,     //pss：层电阻率向量  
                                  //           vector <cdbl> &   wZn，   //层阻抗  
                                                unsigned   m);                 //m：频点数
                                                                                               
        private:
                unsigned   int   n;                                           //地层数  
                vector <dbl>   p;                                             //层电阻率  
                vector <dbl>   h;                                             //层厚度  

                static   const   dbl   u0=1.256e-6;               //磁导率  
};
//-------------------------
Nlayer::Nlayer(unsigned   int   i_n,
                              vector <dbl>   i_p,
                              vector <dbl>   i_h)
                              {
                                      n=i_n;
                                      p=i_p;
                                      h=i_h;
                              }          
//-------------------------    
                inline
        void   Nlayer::caculate_ps(dbl   w,                           //单个角频率  
                                                  dbl   ps)                       //视电阻率
//                                                 vector <cdbl> &   Zn);     //层阻抗  
        {
                  cdbl   k(0.0,sqrt(w*u0/p[n-1]));                                                   //k：波数   ，k2：波数的平方  
                  cdbl   Zom(0,w*u0);                                               //
                  cdbl   ZZ(0,sqrt(w*u0/p[n-1]));                                                       //中间变量
                  cdbl   Zn;  
                  Zom=-Zom/k;
                  for(unsigned   int   i=n-2;   i> =0;   ++i)
                  {
                          Zn=Zom*(1+((ZZ-Zom)/(ZZ+Zom))*(exp(-2.0*k*h[i])))/(1-((ZZ-Zom)/(ZZ+Zom))*(exp(-2.0*k*h[i])));
                          k=(0,sqrt(w*u0/p[i]));

                          Zom=-w*u0/k;
                          ZZ=Zn;
                  }
                  ps=abs(Zn)*abs(Zn)/(w*u0);
                  return;
                 
          }        
//-----------------------------            
inline
void         Nlayer::caculate_ps(vector <dbl> &   ws,       //w：频率向量  
                                                  vector <dbl> &   pss,     //pss：层电阻率向量  
  //                                           vector <cdbl> &   wZn，   //层阻抗  
                                                  unsigned   m)               //m：频点数
          {
                  for(unsigned   i=0;   i <m;   ++i)
                  {
                          dbl   w1=ws[i];
                          dbl   pr=0;
                          this-> caculate_ps(w1,pr);
                          pss[i]=pr;
                  }        
          }        
                       
37   D:\C++学习\C++class\MTlayerRes\Nlayer.h   prototype   for   `   Nlayer::Nlayer(unsigned   int,   std::vector <dbl,   std::allocator <dbl>   > ,    
  error   D:\C++学习\C++class\MTlayerRes\Nlayer.h:11   candidates   are:     Nlayer::Nlayer(const   Nlayer&)  
15   D:\C++学习\C++class\MTlayerRes\Nlayer.h                                     Nlayer::Nlayer(unsigned   int,   std::vector <dbl,   std::allocator <dbl>   > &,    
  D:\C++学习\C++class\MTlayerRes\Nlayer.h   In   member   function   `void     Nlayer::caculate_ps(double,   double) ':  
55   D:\C++学习\C++class\MTlayerRes\Nlayer.h   no   match   for   'operator+ '   in     '1   +   std::operator*(const   std::complex <_Tp> &,   const   std::complex <_Tp> &)    
55   D:\C++学习\C++class\MTlayerRes\Nlayer.h   no   match   for   'operator- '   in     '1   -   std::operator*(const   std::complex <_Tp> &,   const   std::complex <_Tp> &)    

[解决办法]
Nlayer(unsigned i_n,vector <dbl> & i_p,vector <dbl> & i_h);
[解决办法]
你是这样声明，后面定义确不一样..（&）
[解决办法]
Nlayer::Nlayer(unsigned int i_n,
vector <dbl> & i_p,
vector <dbl> & i_h)
:n(i_n),p(i_P),h(i_h){}