 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 981 | 回復(fù): 1 | |||
misamchen新蟲 (初入文壇)
|
[求助]
矩陣里面取數(shù)值相除消掉最高次項后賦值另一個矩陣怎么可以沒有余項
|
|
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 用的輾轉(zhuǎn)相除的辦法，可如附圖算出結(jié)果。 并且自己學(xué)習(xí)寫了一段，然后發(fā)現(xiàn) 輸入N和D = format longg N = [1709000 5505000 1632000 47000 144.99] D = [1 36393000 14591000 779000 8043.2 4.506] >> fostertocauerdiode5(N,D) 結(jié)果從r5那個值開始，后面的項就消不掉了，可是還是要繼續(xù)算，好糾結(jié)。。。 后來發(fā)現(xiàn)從第一次想要減去最高次項的式子里，因為matlab里面除法得到的商再乘以除數(shù)，被被除數(shù)減完還有剩余的最高此項Z(1,:)=U(1,:)-M(1,:).*w(1); 但是我嘗試了從第二次項賦值，還是得不到正確的第二次輾轉(zhuǎn)相除的結(jié)果，求問我這個賦值有什么問題？謝謝！ 現(xiàn)在只需要先寫五階的式子，如果有大神可以幫忙僅僅五階的式子，也十分感激！ 賦值部分如下： （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 謝謝！ |
回復(fù)此樓» 猜你喜歡
316求調(diào)劑
已經(jīng)有6人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
-
一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！
已經(jīng)有3人回復(fù)
-
310求調(diào)劑
已經(jīng)有3人回復(fù)
-
機械專碩299求調(diào)劑至材料
已經(jīng)有4人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有12人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
misamchen
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 20
- 帖子: 2
- 在線: 33分鐘
- 蟲號: 3296786
- 注冊: 2014-06-28
- 專業(yè): 電力電子學(xué)
|
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)%乘完之后減。。。得到第一次的剩余，結(jié)果出現(xiàn)錯誤在這里了 運行如下： 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)) 倒數(shù)第二個結(jié)果就出現(xiàn)錯誤了 但是不知道怎么改 我的s設(shè)的是syms，但是我不是很清楚錯誤出在哪里，請大神指出 不勝感激！ |
2樓2014-07-08 21:05:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 misamchen 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +5 | 77chaselx 2026-03-20 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|