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 的三维图。我主要是不知道已知一组数据怎样在上面的式子求和，都是离散化得。

yufeng

首先确定ab范围 然后利用xy求出r 带入p1 p2 求出p 离散的 话就一个一个求解

wqsoooooooooo

a b P 都是未知数？？
假设a b都为已知的
r =(1+x^2+y^2-2*x*sin(a)*cos(b)-28y*sin(a)sin(b))^1/2;
sin(a)算出来应该是个矩阵
x^2 y^2 算出来是一个数 这样的话 r出来的就是一个矩阵
r1 是一个矩阵r2是一个矩阵………………
这里也不知道想x y 从几到几啊
下面的三个公式看懂了！！
请楼主说详细一些！！！

最好给出数据！！

zhangyong

非常感谢！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

ChaChing

个人水平不高, 还是不太了解LZ的意思!
建议楼主说详细些! 或待高人路过!

wqsoooooooooo

这样求出的好像是个  空间曲线啊 ！

zhangyong

对，就是空间曲线，你可不可以给我个编程过程，这些离散数据实现求和我不会。不知道怎样代入。谢谢！

zhangyong

谢谢各位，我已解决！

ChaChing

怎解决问题的?
可以的话, 回帖分享你的经验和心得!

zhangyong

p1=0;
for i=1:20
  p1=p1+W(i).*cos(x).*cos(110.21.*pi.*sqrt(1+X(i).^2+Y(i).^2-2.*X(i).*sin(x).*cos(y)-2.*Y(i).*sin(x).*sin(y)))./(1+X(i).^2+Y(i).^2-2.*X(i).*sin(x).*cos(y)-2.*Y(i).*sin(x).*sin(y));end;
p2=0;
for i=1:20
  p2=p2+W(i).*cos(x).*sin(110.21.*pi.*sqrt(1+X(i).^2+Y(i).^2-2.*X(i).*sin(x).*cos(y)-2.*Y(i).*sin(x).*sin(y)))./(1+X(i).^2+Y(i).^2-2.*X(i).*sin(x).*cos(y)-2.*Y(i).*sin(x).*sin(y));end;
syms p3
p3=sqrt(p1.^2+p2.^2);

friendchj

syms p3
这一句可以去掉