作者superlaba (吼又濕又冷的天氣)
看板CYCU_CE92B
標題Re: 放著code
時間Fri Jul 18 10:24:55 2008
特徵值特徵向量 a0 a1 a2 a3 矩陣B
沒有在workspace 上顯示成 1*length(dt)
他只有在 command windows 上顯示 1*length(dt)
因為如果我有很多值的話 在workspace 的array editor 上找比較方便
如果在workspace 也太佔空間了
而且在plot(dt,eig(B))的時候就失敗了
下面是正確的
clf
clear
clc
dt=1:1:10;
for I=1:length(dt)
A0(I)=1-dt(I)^2;
A1(I)=2*dt(I);
A2(I)=2*dt(I)^3;
A3(I)=dt(I)^2;
B=[A0(I),A1(I);A2(I),A3(I)];
C(I)=max(eig(B))
end
plot(dt,C)
grid on
hold on
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.121.145.65
※ 編輯: superlaba 來自: 140.121.145.65 (07/18 10:30)
※ 編輯: superlaba 來自: 140.121.145.65 (07/18 10:43)
推 c218420:上篇dt跟[B]維度不同 本就不能直接plot丫 122.116.162.13 07/18 12:21
→ superlaba:所以要先給dt維度 才能plot下去 140.121.145.65 07/18 12:26
→ superlaba:不過我想在plot時候dt取成 10^dt 140.121.145.65 07/18 12:27
→ superlaba:他就會失敗耶 140.121.145.65 07/18 12:27
推 c218420:寫在迴圈裡唄 直接在plot改維度會error 122.116.162.13 07/18 13:00
推 c218420: M(I)=10^dt(I) 122.116.162.13 07/18 13:00
推 c218420:plot(M,C);grid on;xlabel('M');ylabel('C 122.116.162.13 07/18 13:01
→ superlaba:嗯 140.121.145.65 07/18 13:06