声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2678|回复: 10

[编程技巧] matlab离散求和

[复制链接]
发表于 2009-4-23 09:52 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
麻烦各位高手帮帮忙,急用,非常感谢大家帮忙赐教!
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 的三维图。我主要是不知道已知一组数据怎样在上面的式子求和,都是离散化得。
回复
分享到:

使用道具 举报

发表于 2009-4-24 11:05 | 显示全部楼层
首先确定ab范围 然后利用xy求出r 带入p1 p2 求出p 离散的 话就一个一个求解
发表于 2009-4-24 11:08 | 显示全部楼层
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 从几到几啊
下面的三个公式看懂了!!
请楼主说详细一些!!!

最好给出数据!!
 楼主| 发表于 2009-4-24 11:40 | 显示全部楼层

回复 板凳 wqsoooooooooo 的帖子

非常感谢!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
发表于 2009-5-3 12:27 | 显示全部楼层
个人水平不高, 还是不太了解LZ的意思!
建议楼主说详细些! 或待高人路过!
发表于 2009-5-4 09:06 | 显示全部楼层

回复 地板 zhangyong 的帖子

这样求出的好像是个  空间曲线啊 !
 楼主| 发表于 2009-5-6 10:09 | 显示全部楼层

回复 6楼 wqsoooooooooo 的帖子

对,就是空间曲线,你可不可以给我个编程过程,这些离散数据实现求和我不会。不知道怎样代入。谢谢!
 楼主| 发表于 2009-6-1 15:26 | 显示全部楼层
谢谢各位,我已解决!
发表于 2009-6-1 16:31 | 显示全部楼层

回复 8楼 zhangyong 的帖子

怎解决问题的?
可以的话, 回帖分享你的经验和心得!
 楼主| 发表于 2009-7-1 11:22 | 显示全部楼层
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);

评分

1

查看全部评分

发表于 2009-7-1 15:32 | 显示全部楼层

回复 10楼 zhangyong 的帖子

syms p3
这一句可以去掉

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 12:28 , Processed in 0.055569 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表