Redian新闻
>
请教:带参数的积分—公式比较长!
avatar
请教:带参数的积分—公式比较长!# Computation - 科学计算
p*a
1
参数:Xi,eta,对eta积分,积分上下限eta1,eta2是Xi的函数,被积函数G是Xi,eta
的函数,给定Xi范围,846点(前211点一种被积函数形式、积分上下限,其他点另外一
种被积函数形式&积分上下限),目的是分别求每个Xi点对应的积分
程序如下:其中被积函数G复杂,有好几个过渡的表示
clear all;
clc;
syms eta;
Et = 1; r = 100; % transmitter energy; los distance
C = 3e8; Xi = 1.155:0.001:2;
t = Xi*r/C; eta1 = -0.366; eta2 = 0.366; % light speed, spheroidal
coordinate parameters
Ks = 0.231; Ka = 0.746; Ke = Ks+Ka; g = 0.7; % g is mie scattering parameter
K: coefficient
beta_tx = 45*2*pi/360; beta_rx = 45*2*pi/360; %transmitter re
avatar
s*b
2
晕倒。。。
一个玩具而已,还真拿matlab来做积分阿

eta

【在 p*******a 的大作中提到】
: 参数:Xi,eta,对eta积分,积分上下限eta1,eta2是Xi的函数,被积函数G是Xi,eta
: 的函数,给定Xi范围,846点(前211点一种被积函数形式、积分上下限,其他点另外一
: 种被积函数形式&积分上下限),目的是分别求每个Xi点对应的积分
: 程序如下:其中被积函数G复杂,有好几个过渡的表示
: clear all;
: clc;
: syms eta;
: Et = 1; r = 100; % transmitter energy; los distance
: C = 3e8; Xi = 1.155:0.001:2;
: t = Xi*r/C; eta1 = -0.366; eta2 = 0.366; % light speed, spheroidal

avatar
p*a
3
对时间要求不高,只要能够算出正确的结果就成。

【在 s**b 的大作中提到】
: 晕倒。。。
: 一个玩具而已,还真拿matlab来做积分阿
:
: eta

avatar
s*b
4
运行巨慢无比,我还以为你等不及了。
anyway,程序没错就行,matlab就是太慢,特别是带参数积分多冲击分之类的

【在 p*******a 的大作中提到】
: 对时间要求不高,只要能够算出正确的结果就成。
avatar
p*a
5
对前212个Xi点确实没有错误,只有警告,第213个点开始报错,其实和前面的公式都一
样,就是不知道咋错了,能指点下吗?

【在 s**b 的大作中提到】
: 运行巨慢无比,我还以为你等不及了。
: anyway,程序没错就行,matlab就是太慢,特别是带参数积分多冲击分之类的

avatar
s*b
6
没看懂你为啥要分成两个循环

【在 p*******a 的大作中提到】
: 对前212个Xi点确实没有错误,只有警告,第213个点开始报错,其实和前面的公式都一
: 样,就是不知道咋错了,能指点下吗?

avatar
p*a
7
因为被积函数和积分上下限不同。
我试过跑一个Xi点,对于Xi=1:212,有警告,但能得到值
从Xi=213开始,就又错了。

【在 s**b 的大作中提到】
: 没看懂你为啥要分成两个循环
avatar
s*b
8
你还是贴算法吧,把你的思路写出来,不然要看懂你的code还真的花点时间

【在 p*******a 的大作中提到】
: 因为被积函数和积分上下限不同。
: 我试过跑一个Xi点,对于Xi=1:212,有警告,但能得到值
: 从Xi=213开始,就又错了。

avatar
p*a
9


【在 s**b 的大作中提到】
: 你还是贴算法吧,把你的思路写出来,不然要看懂你的code还真的花点时间
avatar
p*a
10
公式比较绕人,555.
我就是想把那个积分积出来,被积函数中的各项可以归结为Xi、eta,积分上下限的选
择依赖于Xi的取值。
指定了Xi的范围,讨论好积分上下限,直接调用int积分的,对于前212个Xi值,只有警
告,有值,后面的就出错了。
avatar
p*a
11
最后还是算出来了,虽然有警告。
问题处在积分区间没整对,多谢。

【在 p*******a 的大作中提到】
: 公式比较绕人,555.
: 我就是想把那个积分积出来,被积函数中的各项可以归结为Xi、eta,积分上下限的选
: 择依赖于Xi的取值。
: 指定了Xi的范围,讨论好积分上下限,直接调用int积分的,对于前212个Xi值,只有警
: 告,有值,后面的就出错了。

相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。