24小时热门版块排行榜

返回列表

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

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiaomeim

新虫 (初入文坛)

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

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

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

@laosam280 发自小木虫Android客户端

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有21人回复
导师想让我从独立一作变成了共一第一已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
申请2026年博士已经有6人回复

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

已阅回复此楼关注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)