| 查看: 1018 | 回复: 1 | ||
[求助]
矩阵里面取数值相除消掉最高次项后赋值另一个矩阵怎么可以没有余项
|
|
A*s^4 + B*s^3 + C*s^2 D*s + E / F*s^5 + G*s^4+ H*s^3 + I*s^2 + J*s + 1 = =1/(c5*s + 1/(r5 + 1/(c4*s + 1/(r4 +1/(1/(r3 + 1/(c2*s + 1/(r2 + 1/(c1*s + 1/r1)))) + c3*s))))) 1 = --------------------------------------------------------- 1 c5s + -------------------------------------------------- 1 r5 + --------------------------------------------- 1 c4 s +-------------------------------------- 1 r4 + --------------------------------- 1 c3 s + -------------------------- 1 r3 +--------------------- 1 c2 s +-------------- 1 r2 +--------- 1 c1 s+ -- r1 用的辗转相除的办法，可如附图算出结果。  并且自己学习写了一段，然后发现 输入N和D = format longg N = [1709000 5505000 1632000 47000 144.99] D = [1 36393000 14591000 779000 8043.2 4.506] >> fostertocauerdiode5(N,D) 结果从r5那个值开始，后面的项就消不掉了，可是还是要继续算，好纠结。。。 后来发现从第一次想要减去最高次项的式子里，因为matlab里面除法得到的商再乘以除数，被被除数减完还有剩余的最高此项Z(1,:)=U(1,:)-M(1,:).*w(1); 但是我尝试了从第二次项赋值，还是得不到正确的第二次辗转相除的结果，求问我这个赋值有什么问题？谢谢！ 现在只需要先写五阶的式子，如果有大神可以帮忙仅仅五阶的式子，也十分感激！ 赋值部分如下： （from line 45） w(1)=U(1,6)/M(1,5); w(1) U(1,6)=0; M(1,6)=0; % M(1,5)=0; U (2,:) = M (1,:); Z(1,:)=U(1,:)-M(1,:).*w(1); % Z(1,6)=0; U (2,:) = M (1,:); M (2,:) = Z (1,:); w(2)= U (2,5)/M (2,5); w(2) U(2,5)=0; M(2,5)=0; Z(2,:)=U(2,:)-M(2,:).*w(1); 整体code如下， 当然如果只帮忙看看5阶的代码也感谢！ :) function g = fostertocauerdiode5b(N,D) %clc x=size(N);y=size(D); L=size(N)-size(D); if (L(1,2)>0) D((y(1,2)+1):x(1,2))=0; end if (L(1,2)<0) N((x(1,2)+1):y(1,2))=0; end clear x y L; syms s; L=size(D); P=L(1,2); j=1;q=1; for i=1:P if N(i)~=0 A(j)=N(i)*s^(i-1); j=j+1; end if D(i)~=0 B(q)=D(i)*s^(i-1); q=q+1; end end M(1,:)=A(1,:); % The matrix built by numerator after the first calculation............... U(1,:)=B(1,:); %The matrix built by numerator after the first calculation............... w(1)=U(1,6)/M(1,5); W(1)=w(1,1); %..................... W(1) x=size(M);y=size(U); L=size(M)-size(U); if (L(1,2)>0) U((y(1,2)+1):x(1,2))=0; end if (L(1,2)<0) M((x(1,2)+1):y(1,2))=0; end clear x y L; Z(1,:)=U(1,:)-M(1,:)*w(1);% %z(1)=sum(U(1,:)')'-sum(M(1,:)')'*W(1); k=2; while (P-1)~=0 % while sum(Z(k-1,:)')~=0 % sum(A) treats the columns of it as vectors, returning a row vector of the sums of each column.Then it shows the sum of each row % L=size(Z); U(k,:)=M(k-1,:); %clear M; M(k,:)=Z((k-1),:); % clear L; w(k)=U(k,(P-1))/M(k,(j-1)); W(k)=w(k); % x=size(M);y=size(U); % L=size(M)-size(U); % if (L(1,2)>0) % U(:,(y(1,2)+1):x(1,2))=0; % end % if (L(1,2)<0) % M(:,(x(1,2)+1):y(1,2))=0; % end % clear x y; % % clear x y L; Z(k,:)=Z((k-1),:)-U(k,:)*W(k); W(k) if P == j j=j-1; end if P ~= j P=P-1; end k=k+1; end end 谢谢！ |
回复此楼» 猜你喜欢
关于Filecode分析方法
已经有6人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
|
w(1)=U(1,6)/M(1,5); % 第一次计算电容值是正确的 w(1) % 输出电容值是对的 % U(1,6)=0; 开始尝试想赋值，但是感觉第六列没什么帮助 % M(1,6)=0; Z(1,6)= [s s s s s s];% 完全是为了想知道Z里面是怎么回事 Z(1,  =U(1,-M(1,.*w(1); % 得到分子分离出电容部分之后剩下的部分M(1,5)% 完全是为了想知道M是怎么回事 M(1,5).*w(1);% 完全是为了想知道乘完以后是怎么回事 Z(1,5)%乘完之后减。。。得到第一次的剩余，结果出现错误在这里了 运行如下： ans = (751*s)/24165 ans = (14499*s^4)/100 ans = (40216*s^4)/5 - (2253*s^5)/500 ans = (14499*s^4)/(100*((40216*s^4)/5 - (2253*s^5)/500)) 倒数第二个结果就出现错误了 但是不知道怎么改 我的s设的是syms，但是我不是很清楚错误出在哪里，请大神指出 不胜感激！ |
2楼2014-07-08 21:05:28
