|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
以下程序:我想要的F应该是1*8的数组,结果显示是1*1的,我不知道哪里出错了,高手请指教一下,多谢<BR>A=[1 0.3 0.5 1];B=[0.6 0.5 0.4 0.7];C=[];<BR>for i=(1:length(A))<BR> C=[C,abs(B.^2/(1-B.^2))];<BR> l=ceil(max(C));<BR> n=length(A)*(l+1)-1;<BR>end;l,n<BR>product=1,y2=[],<BR>syms b;<BR>for i=1:length(A) <BR> y1=A(i).*((B(i)-b).^l);<BR> for j=1:length(A)<BR> if (j~=i)<BR> product=product.*((1-b.*B(j)).^(l+1)); <BR> y2=[y2,y1.*product];<BR> end<BR> end;y2<BR> end<BR>sum=0; F=[]; <BR>for m=0:n <BR> for i=1:length(A), <BR> sum=sum+diff(y2(i),b,m);<BR> end<BR>end<BR> for b=0;<BR> F=[F,(1/prod(m).*sum)]; <BR> end;F |
|