| 查看: 941 | 回复: 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 谢谢！ |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有8人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
|
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