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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

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

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

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

谢谢 不过好像不是我要的哦 我的是求f/g（f,g均为含x函数）的N阶导数

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

syms x;f=3*x+2;g=2*x.^2+x+3;diff(f/g,10) %求f/g的10次导数 ans = -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 以上是一个具体的事例！！

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

刚自己写了个比较烂，有错误！哪位帮我改下 function ans=P1(f,g,n) if n==1 A=diff(g)*f-g*diff(f) B=f^2 else for i=2:n; f=A g=B A=diff(g)*f-g*diff(f) B=f^2 end; end; H=A/B 刚问老师，说是符号运算，不是数值运算！这个程序错的 郁闷，没头绪了！有能力的快帮忙啊

[此贴子已经被作者于2006-4-29 10:43:07编辑过]

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

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

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

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

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

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

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

大虾们快帮忙啊

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

学习

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

必须写成子程序调用的形式吗？

账号		自动登录	找回密码
密码			我要加入

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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