曲面积分数值积分方法
34视图(30天)
显示旧的评论
你好。我必须执行这个积分:
这另一个基本上是一样的:
地点:
这另一个基本上是一样的:
地点:
- 一个是一个已知的复杂的矢量,是一个真正的导数的FFT实验室mesuerement吗
是已知的- 年代是一个半径为R的圆周表面(已知)
自
一个
是一个复杂的数值向量,我试图用cumtrapz和积分函数来计算这个积分。错了两个方法给我结果和结果取决于使用的方法不同。
此外,我看到最终结果的方法2积分函数是相同的第一个积分(pa_a)和第二个(p)。
你能帮我找到正确的方式来执行这个计算吗?
% %数据
%的输入
y =数据;
%计算加速度的测量速度
一个= diff (y);
一个= fft (A);
%的测量数据
Fs = 44100;%测量采样frquency
f_start = 2;%初始测量频率(赫兹)
f_end = 5000;%的最终测量频率(赫兹)
%的参数
rho_0 = 1.225;%公斤/立方米
ra = 1;% m
rc = 0;% m
p_0 = 1;%的标准压力
c = 340;%的声音速度米/秒
R = 0.254;%半径
Sc =π* R ^ 2;%膜表面(平方米)
f = linspace (f_start f_end、长度(diff (fft (y))));%频率向量
k = 2 *π* f / c;%波数
% %积分计算
%方法1:梯形集成
domain1 = linspace(0, 2 *π,长度(A));
domain2 = linspace (0, R,长度(A));
pa_a_cumtrapz = rho_0 /(2 *π)* cumtrapz (domain1, cumtrapz (domain2、A / abs (ra-rc)));
p_cumtrapz = rho_0 /(2 *π)* cumtrapz (domain1, cumtrapz (domain2、A / abs (ra-rc)。* exp(我* k * abs (ra-rc))));
% - - - - - - - - - - - - - - - - - - - - - - - -
%方法2:二重积分(我不能使用integral2因为我有一个向量,而不是一个矩阵)
%第一内部积分计算
被积函数= @ (om) / abs (ra-rc);
temp =积分(被积函数0 R,“ArrayValued”,真正的);
%第二积分计算
integrand2 = @ (om)温度;
pa_a_integral = rho_0 /(2 *π)。*积分(integrand2 0 2 *π,“ArrayValued”,真正的);
p_integral = rho_0 /(2 *π)*积分(integrand2 0 2 *π,“ArrayValued”我,真的)。* exp (1 * k * abs (ra-rc));
% %集成结果的情节
图;
重对数(abs (pa_a_cumtrapz));
持有在
重对数(abs (p_cumtrapz) * 0.3)% * 0.3是一个转向检查两个光谱之间的区别
传奇(“pa_a cumtrapz”,“p cumtrapz低转移”)
xlim ([f_start f_end]);
网格小
% - - - - - - - - - - - - - - - - - - - - - - - -
图;
重对数(abs (pa_a_integral));
持有在
重对数(abs (p_integral) * 0.3)% * 0.3是一个转向检查两个光谱之间的区别
传奇(“pa_a积分”,“p积分低转移”)
xlim ([f_start f_end]);
网格小
0评论
答案(0)
翻译的
也seleccionar uno de来说《国家/语言: