[求助]请教数组求导的问题

jifang29 · 发表于 2006-3-30 18:03

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

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

x

请问数组求导可以直接用diff(f,x,n)命令吗？

happy · 发表于 2006-3-30 18:03

可以，doc diff看例子

jifang29 · 发表于 2006-3-30 18:07

doc diff???

happy · 发表于 2006-3-30 18:08

在命令窗口输入doc diff你就可以看到这个函数的帮助了

jifang29 · 发表于 2006-3-31 09:27

A=[1 0.3 0.5 1];B=[0.6 0.5 0.4 0.7];C=[]; for i=(1:length(A)) C=[C,abs(B.^2/(1-B.^2))]; l=ceil(max(C)); n=length(A)*(l+1)-1; end;l,n product=1,y2=[], syms b; for i=1:length(A) y1=A(i).*((B(i)-b).^l); for j=1:length(A) if (j~=i) product=product.*((1-b.*B(j)).^(l+1)); y2=[y2,y1.*product]; end end;y2 end sun=0; F=[]; for i=1:length(A), m=0:n f=diff(y2(i),b,m); sum=sum+f; end;sum for b=0; F=[F,(1/prod(m)*sum)]; end;F 出现以下错误??? ??? Error using ==> sym/diff Do not recognize argument number 3
Error in ==> E:\study\filter\laguerrenew1.m On line 20 ==> f=diff(y2(i),b,m);

happy · 发表于 2006-3-31 20:04

<DIV class=quote>以下是引用jifang29在2006-3-31 9:27:11的发言： 
A=[1 0.3 0.5 1];B=[0.6 0.5 0.4 0.7];C=[]; for i=(1:length(A)) C=[C,abs(B.^2/(1-B.^2))]; l=ceil(max(C)); n=length(A)*(l+1)-1; end;l,n product=1,y2=[], syms b; for i=1:length(A) y1=A(i).*((B(i)-b).^l); for j=1:length(A) if (j~=i) product=product.*((1-b.*B(j)).^(l+1)); y2=[y2,y1.*product]; end end;y2 end sun=0; F=[]; for i=1:length(A), m=0:n f=diff(y2(i),b,m); sum=sum+f; end;sum for b=0; F=[F,(1/prod(m)*sum)]; end;F 出现以下错误??? ??? Error using ==> sym/diff Do not recognize argument number 3
Error in ==> E:\study\filter\laguerrenew1.m On line 20 ==> f=diff(y2(i),b,m);</DIV>
 刚才看前面的帖子，这个错误你已经改过来了，不用再回复了吧

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

[求助]请教数组求导的问题

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

回复：（jifang29）[求助]请教数组求导的问题

回复：（jifang29）[求助]请教数组求导的问题

[求助]请教数组求导的问题

回复：（jifang29）[求助]请教数组求导的问题