24小时热门版块排行榜

返回列表

vs570588

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3112.1
散金: 456
红花: 1
帖子: 923
在线: 420.6小时
虫号: 822119
注册: 2009-08-04
性别: GG
专业: 环境工程

[求助] 大家看看这个matlab中自定义函数的意思

function Y=mynf1(beta,t)
x=213.63;
for j=1:size(beta,1)
q(j)=beta(j,1);
K(j)=(beta(j,2))^2;
for i=1:length(t)
Y(j,i)=fsolve(@(S)(K(j)*log(S)+S+x*q(j)*t(i)-(3.8106*K(j)+x*q(j)*19.9167))/130.6542-1,120);
end;
end;
Y=Y';
大家帮忙，看看这个函数，最好能给我解释每一步，我matlab一点基础都没有。for语句定义的向量函数，我实在不懂。

回复此楼

» 猜你喜欢

浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
物理学I论文润色/翻译怎么收费? 已经有62人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有1人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复
apl投稿状态已经有2人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab的函数调用问题已经有3人回复
MATLAB中产生高斯白噪声，涉及到awgn和wgn函数已经有3人回复
请教个matlab中的rand函数已经有3人回复
【求助】大家看看，matlab安装过程的截图，找出错已经有14人回复
大家帮我看看这个拒稿信是什么意思？已经有10人回复
【求助】精通matlab最优化计算，一书中的funval函数是什么意思已经有10人回复
【求助】利用matlab求解微分方程得出解函数中含0的意思已经有3人回复
【求助】大家帮我看看这个MATLAB程序？已经有6人回复

1楼 2011-06-21 16:07:41

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

libralibra

至尊木虫 (著名写手)

骠骑将军

计算强帖: 3
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

★ ★ ★
vs570588(金币+2): 2011-06-23 20:25:02
sunyang1988(金币+3): 谢谢指点 2011-06-25 12:42:15
ben_ladeng: 2011-08-17 00:35:54
vs570588(金币+2): 2011-08-19 21:17:34

应该是根据输入参数解方程的,大概解释下:

CODE:

% 子函数:

% 输入: 

%   beta: 矩阵

%   t:  向量

% 输出:

%   Y:  矩阵

function Y=mynf1(beta,t)

% 常数x

x = 213.63;

% 循环1,次数等于beta矩阵行数

for j=1:size(beta,1)

    % q等于beta第一列

    q(j)=beta(j,1);

    % K等于beta第二列的平方

    K(j)=(beta(j,2))^2;

    % 循环2,次数等于t长度

    for i=1:length(t)

        % Y第j行第i个元素等于解出来的S

        Y(j,i)=fsolve(@(S)(K(j)*log(S)+S+x*q(j)*t(i)-(3.8106*K(j)+x*q(j)*19.9167))/130.6542-1,120);

    end;

end;

% Y转置输出

Y=Y';

end

其实此函数有待优化,应避免操作变长向量,而应该一次赋值,加快处理速度
修改版如下(其实q和K这2个中间变量完全可以不用,只是在解方程函数fsolve中用beta(j,1)和beta(j,2)^2代替即可):

CODE:

% 子函数:

% 输入: 

%   beta: 矩阵

%   t:  向量

% 输出:

%   Y:  矩阵

function Y=mynf1(beta,t)

% 常数x

x = 213.63;

% q和k一次性赋值:q等于beta第一列;K等于beta第二列的平方

q = beta(:,1);

K = beta(:,2).^2;

% 循环1,次数等于beta矩阵行数

for j=1:size(beta,1)

    % 循环2,次数等于t长度

    for i=1:length(t)

        % Y第j行第i个元素等于解出来的S

        Y(j,i)=fsolve(@(S)(K(j)*log(S)+S+x*q(j)*t(i)-(3.8106*K(j)+x*q(j)*19.9167))/130.6542-1,120);

    end;

end;

% Y转置输出

Y=Y';

end