24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者xiaomeim将赠送您 29 个金币

xiaomeim

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 31.3
散金: 10
帖子: 17
在线: 4.9小时
虫号: 24571323
注册: 2020-11-13
专业: 生态学

[求助] 请各位大神伸出援手帮我解个题（有偿）已有1人参与

想求出f effective(x)，请问有没有大神可以帮忙算一算这个公式？能够用R或C语言算出来最好，手算也行，提供思路也可以（有偿）。感谢感谢感谢！！！

@laosam280 发自小木虫Android客户端

回复此楼

» 猜你喜欢

1楼 2022-06-02 22:35:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

du660283

禁虫 (小有名气)

本帖内容被屏蔽

2楼2022-06-04 05:53:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyh700600

金虫 (小有名气)

应助: 47 (小学生)
金币: 9368.2
红花: 10
帖子: 266
在线: 137.1小时
虫号: 56711
注册: 2005-01-16
专业: 金融工程

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

  
%  用 Matlab 2009a 求解

syms x y

r = sqrt(x*x+y*y)

f_wald = sqrt(3/(2*pi*r^3))*exp( -3*(r-2)^2/(8*r) )

f_eff_1 = f_wald/(2*pi*r)

g = matlabFunction(f_eff_1)

f_eff = @(x) quadgk( @(y) g(x,y) ,-inf,inf)

xx = -2:0.01:2;  %  x坐标

fxx = arrayfun( @(x) f_eff(x) ,xx);  %  f_effective(x)

figure

plot(xx,fxx)