[求助]想求e指数函数的N次导数 - MATLAB爱好者论坛-LabFans.com

fanxing39 · 2008-06-04, 15:09

syms s p N r
g=exp(-r*(1- s*(1-p+p*s)/((1+p)^(N-1)+s*(1-p+p*s)*(1-(1+p)^(N-1)))));
L1=zeros(100,1); j=1;
g=diff(g,'s');
r=1;s=0; N=200;p=10^(1/200)-1;
f=subs(g);
n=1:100;
L1=f./factorial(n);
plot([1:100],L1);

解释： factorial(n) 就是求阶乘的，你的f 原来只需要算一次，所以放到你的for 循环里，就多计算了99次，增加了它的计算量，并且阶乘也是你自己编程算，这个没有直接调用函数算起来快
在计算一个程序前，如果以前的程序结果你不需要调用，建议加 “clear” （清除内存中的变量）

主题工具
显示可打印版本邮寄本页给好友
显示模式
切换到平板模式切换到混合模式树形模式

相似的主题
主题	主题作者	版面	回复	最后发表
[求助]3维图的X，Y 轴可以随时间变化吗	莱茵河流浪	MATLAB论坛	0	2008-05-18 04:38
请问这个方程可以用MATLAB解没有?	thb296	MATLAB论坛	6	2008-04-28 19:02
求助怎样用matlab计算曲面或球体积分	tangyongliang	MATLAB论坛	0	2008-04-05 22:32
【求助】紧急求助关于simulink的专家	chenlinysu	MATLAB论坛	2	2007-07-30 16:21
【求助】紧急求助	chenlinysu	MATLAB论坛	1	2007-06-25 18:54

2008-06-04, 15:09	#5
fanxing39 高级会员注册日期: 2007-12-02 年龄: 44 帖子: 303 声望力: 31	syms s p N r g=exp(-r(1- s(1-p+ps)/((1+p)^(N-1)+s(1-p+ps)(1-(1+p)^(N-1))))); L1=zeros(100,1); j=1; g=diff(g,'s'); r=1;s=0; N=200;p=10^(1/200)-1; f=subs(g); n=1:100; L1=f./factorial(n); plot([1:100],L1); 解释： factorial(n) 就是求阶乘的，你的f 原来只需要算一次，所以放到你的for 循环里，就多计算了99次，增加了它的计算量，并且阶乘也是你自己编程算，这个没有直接调用函数算起来快在计算一个程序前，如果以前的程序结果你不需要调用，建议加 “clear” （清除内存中的变量）