24小时热门版块排行榜

>论坛更新日志 (3409)
>虫友互识 (656)
>导师招生 (184)
>文献求助 (167)
>论文投稿 (100)
>硕博家园 (89)
>考博 (89)
>基金申请 (64)
>教师之家 (63)
>博后之家 (44)
>公派出国 (40)
>考研 (39)
>论文道贺祈福 (32)
>找工作 (29)
>休闲灌水 (28)
>外文书籍求助 (16)

返回列表

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语句定义的向量函数，我实在不懂。

回复此楼

» 猜你喜欢

石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有235人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有1人回复
投稿chemical physical letters不送审？已经有2人回复
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生已经有13人回复

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

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