matlab离散求和问题

zhangyong

麻烦各位高手帮帮忙，急用，非常感谢大家帮忙赐教！
r =(1+x^2+y^2-2*x*sin(a)*cos(b)-28y*sin(a)sin(b))^1/2;

      N
p1=∑    (wm*cos(a)*cos(k*rm))/rm^2;
    m=1
      N
p2=∑  ( wm*cos(a)*sin(k*rm))/rm^2;
     m =1
p=(p1^2+p2^2)^1/2;
其中a和b是p的函数，X，Y和wm为已知的一组数（大约几百个,由ansys模拟得出），对应的x,y有对应的wm,m为下标，k为常数。 求a，b和P 的三维图。我主要是不知道已知一组数据怎样在上面的式子求和，都是离散化得。


  x,y为已知，a,b未知，a是(-pi,pi)；b(0,2*pi).k=100*pi,求a,b与p的三维图。数据如下：
NODE      X                 Y                      wm           
       1 -0.15439    -0.15439     0.56647      
       2 -0.48057E-01-0.32269      1.2901         
       3  0.10272    -0.27409      1.0960         
       4  0.18846    -0.15733E-01 0.63840E-01   
       5  0.11820     0.25523     -1.0228         
       6 -0.36382E-01 0.32817     -1.3129         
       7 -0.16703     0.14764    -0.59191        
       8 -0.16703    -0.14764     0.59191        
       9 -0.36382E-01-0.32817      1.3129         
      10  0.11820    -0.25523      1.0228      
      11  0.18846     0.15733E-01-0.63840E-01   
      12  0.10272     0.27409     -1.0960         
      13 -0.48057E-01 0.32269     -1.2901         
      14 -0.15439     0.15439    -0.56647      
      15 -0.32269    -0.48057E-01  1.2901         
      16 -0.10391    -0.10391      2.8761         
      17  0.22366    -0.87927E-01  2.4377         
      18  0.39964    -0.46894E-02 0.13975        
      19  0.25523     0.82236E-01 -2.2749        
      20 -0.79198E-01 0.10593     -2.9238