candytbs 发表于 2006-4-28 20:57

【帮忙】谁能编个求f/g（f,g均为含x函数）的N阶导数的Matlab程序

谁能编个求f/g（f,g均为含x函数）的N阶导数的Matlab程序，给我借鉴下啊 谢谢了

zhanglei_1 发表于 2006-4-28 22:55

<P>无知者无畏<BR> <BR>function N_diff=n_diff(p,n)   %"p"待求导函数,"n"求导次数<BR>p1=p;<BR>i=1;<BR>while(i&lt;=n)<BR>   p1=diff(p1);<BR>   i=i+1;<BR> end<BR> N_diff=p1;<BR>    </P>

candytbs 发表于 2006-4-28 23:07

<P>谢谢<BR>不过好像不是我要的哦<BR>我的是<STRONG><FONT color=#000000>求f/g（<FONT color=#ff0000>f,g均为含x函数</FONT>）的N阶导数</FONT></STRONG></P>

lang473036 发表于 2006-4-28 23:08

<STRONG><FONT color=#2b2bd5>syms x;f=3*x+2;g=2*x.^2+x+3;diff(f/g,10)%求f/g的10次导数<BR></FONT> <BR><FONT color=#33cc70>ans =<BR> <BR>-10886400/(2*x^2+x+3)^10*(4*x+1)^9+174182400/(2*x^2+x+3)^9*(4*x+1)^7-914457600/(2*x^2+x+3)^8*(4*x+1)^5+1741824000/(2*x^2+x+3)^7*(4*x+1)^3-870912000/(2*x^2+x+3)^6*(4*x+1)+3628800*(3*x+2)/(2*x^2+x+3)^11*(4*x+1)^10-65318400*(3*x+2)/(2*x^2+x+3)^10*(4*x+1)^8+406425600*(3*x+2)/(2*x^2+x+3)^9*(4*x+1)^6-1016064000*(3*x+2)/(2*x^2+x+3)^8*(4*x+1)^4+870912000*(3*x+2)/(2*x^2+x+3)^7*(4*x+1)^2-116121600*(3*x+2)/(2*x^2+x+3)^6</FONT><BR><FONT color=#f70909>以上是一个具体的事例！！</FONT></STRONG>

candytbs 发表于 2006-4-28 23:09

刚自己写了个比较烂，有错误！哪位帮我改下<br><br>function ans=P1(f,g,n)<br>if n==1<br>    A=diff(g)*f-g*diff(f)<br>    B=f^2<br>else<br>    for i=2:n;<br>      f=A<br>      g=B<br>      A=diff(g)*f-g*diff(f)<br>      B=f^2            <br>    end;<br>end;<br>H=A/B<br><br><br>刚问老师，说是符号运算，不是数值运算！这个程序错的<br>郁闷，没头绪了！有能力的快帮忙啊
[此贴子已经被作者于2006-4-29 10:43:07编辑过]

candytbs 发表于 2006-4-28 23:13

4楼的同样感谢，虽然方法不好

lang473036 发表于 2006-4-28 23:18

呵呵！！<BR>我只能写一些肤浅的东东阿！！<BR>像你们写的function我只懂一点点阿<BR>我是一个初学matlab的小鬼阿<BR>

zhanglei_1 发表于 2006-4-29 09:35

<P>我也是菜鸟,不过相信努力能造就成绩,不过到高手的过程,需要时间和经验!</P>

candytbs 发表于 2006-4-29 14:08

<P>大虾们快帮忙啊</P>

edwardflr 发表于 2006-4-29 18:22

学习

happy 发表于 2006-4-29 21:46

回复：（candytbs）【帮忙】谁能编个求f/g（f,g均为含...

必须写成子程序调用的形式吗？
页: [1]
查看完整版本: 【帮忙】谁能编个求f/g（f,g均为含x函数）的N阶导数的Matlab程序