matlab离散求和
麻烦各位高手帮帮忙，急用，非常感谢大家帮忙赐教！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 的三维图。我主要是不知道已知一组数据怎样在上面的式子求和，都是离散化得。 首先确定ab范围 然后利用xy求出r 带入p1 p2 求出p 离散的 话就一个一个求解 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 从几到几啊
下面的三个公式看懂了！！
请楼主说详细一些！！！
最好给出数据！！
回复 板凳 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
30.10272 -0.27409 1.0960
40.18846 -0.15733E-01 0.63840E-01
50.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
100.11820 -0.25523 1.0228
110.18846 0.15733E-01-0.63840E-01
120.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-011.2901
16 -0.10391 -0.10391 2.8761
170.22366 -0.87927E-012.4377
180.39964 -0.46894E-02 0.13975
190.25523 0.82236E-01 -2.2749
20 -0.79198E-01 0.10593 -2.9238 个人水平不高, 还是不太了解LZ的意思!
建议楼主说详细些! 或待高人路过!
回复 地板 zhangyong 的帖子
这样求出的好像是个空间曲线啊 ！回复 6楼 wqsoooooooooo 的帖子
对，就是空间曲线，你可不可以给我个编程过程，这些离散数据实现求和我不会。不知道怎样代入。谢谢！ 谢谢各位，我已解决！回复 8楼 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);
回复 10楼 zhangyong 的帖子
syms p3这一句可以去掉
页:
[1]