声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 918|回复: 1

[求助]请帮忙看以下程序哪里错了

[复制链接]
发表于 2006-3-31 11:41 | 显示全部楼层 |阅读模式

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

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

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
回复
分享到:

使用道具 举报

发表于 2006-3-31 19:35 | 显示全部楼层

回复:(jifang29)[求助]请帮忙看以下程序哪里错了

<DIV class=quote><B>以下是引用<I>jifang29</I>在2006-3-31 11:41:13的发言:</B><BR>以下程序:我想要的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><FONT color=#ff0000>   for b=0;<BR>    F=[F,(1/prod(m).*sum)]; <BR>   end;F</FONT></DIV>
<P>这是什么啊?是不是忘了什么东西</P>
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-20 21:30 , Processed in 0.054357 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表