怎么用matlab生成三对角矩阵？
材料/工具：电脑、matlab软件
1、首先打开电脑上的“matlab”软件，在命令行窗口输入一个向量n，向量n有4个元素。接着使用diag函数生成对角矩阵。
2、diag(n,k)可以把向量放在第k条对角线上，k为正值，表示右上。k为负值，则表示左下。在命令行输入diag(n,2)，即可把向量放在右上的第二条对角线上。
3、在命令行输入diag(n,-2)，即可把向量放在左下的第二条对角线上，运行结果。
4、接下来在命令行输入一个矩阵，此处以随机矩阵为例，输入rand(3)生成3*3的随机矩阵A，结果。
5、下面使用diag函数提取矩阵A的对角线元素，对角线元素个数为3，提取的对角线元素。
6、使用命令diag(A,k)即可获得第k条对角线上的元素组成的向量，k为正值，表示右上。
7、k为负值，则表示左下。输入命令diag(A,1)即可获得右上的第1条对角线上的元素。
8、下面输入命令diag(A,-1)即可获得左下的第1条对角线上的元素。具体运行结果。
9、如果要根据矩阵的对角线元素生成对角矩阵，可使用两个diag函数，输入命令diag(diag(A))，即可使用A的对角线元素生成对角矩阵。
10、blkdiag函数可以根据多个矩阵生成准对角矩阵，此处使用rand(3)和ones(2)，生成准对角矩阵，运行结果。
绛旓細1銆侀鍏堟墦寮鐢佃剳涓婄殑鈥渕atlab鈥濊蒋浠讹紝鍦ㄥ懡浠よ绐楀彛杈撳叆涓涓悜閲弉锛屽悜閲弉鏈4涓厓绱銆傛帴鐫浣跨敤diag鍑芥暟鐢熸垚瀵硅鐭╅樀銆2銆乨iag(n,k)鍙互鎶婂悜閲忔斁鍦ㄧk鏉″瑙掔嚎涓婏紝k涓烘鍊硷紝琛ㄧず鍙充笂銆俴涓鸿礋鍊硷紝鍒欒〃绀哄乏涓嬨傚湪鍛戒护琛岃緭鍏iag(n,2)锛屽嵆鍙妸鍚戦噺鏀惧湪鍙充笂鐨勭浜屾潯瀵硅绾夸笂銆3銆佸湪鍛戒护琛岃緭鍏iag(n,-...
绛旓細diag(1:9, 1) + diag(2:10, -1) + eye(10)
绛旓細濡備綍鐢╩atlab鍒涘缓涓涓猲琛宯鍒楀瑙掔嚎涓10,绗竴涓婂瑙掔嚎涓1绗竴涓嬪瑙掔嚎涓2鐨涓夊瑙掔煩闃碉紵棰樹富鐨勫瑙掔煩闃甸棶棰樺彲浠ヤ娇鐢╠iag锛堬級鍑芥暟鏉ュ疄鐜般傝鍑芥暟鐨勪娇鐢ㄦ牸寮忔槸 D = diag(v,k)涓婅堪鏍煎紡鐨勬剰鎬濇槸锛 灏嗗悜閲 v 鐨勫厓绱犳斁缃湪绗 k 鏉″瑙掔嚎涓娿俴=0 琛ㄧず涓诲瑙掔嚎锛宬>0 浣嶄簬涓诲瑙掔嚎涓婃柟锛宬<0 浣嶄簬...
绛旓細diag(A)+diag(B,1)+diag(C,-1)
绛旓細浣跨敤Lanczos绠楁硶鐨勫叿浣撹緭鍏ユ牸寮忎负涓涓嚱鏁拌皟鐢锛屽lanczos(matrix, vector, nmax)銆傝繖閲岀殑matrix鏄綘闇瑕佷笁瀵硅鍖栫殑鐭╅樀锛屼緥濡俒1 2 3;4 5 6;7 8 9]锛泇ector鏄綘鎻愪緵鐨勫悜閲忥紝濡俒1;1;1]锛沶max鍒欐槸浣犳兂瑕佺殑鐭╅樀缁村害锛屽涓婅堪鐨12銆傝皟鐢ㄨ繖涓嚱鏁板悗锛岀畻娉曚細杩斿洖涓涓12x12鐨勪笁瀵硅鐭╅樀锛岃繖鏄畻娉曡繍琛...
绛旓細N=5; %鏈鍚庣煩闃靛皢鏄3N*3N澶у皬 I0=eye(3);A0=magic(3); %杩欓噷鐢ㄩ瓟鏂圭煩闃垫祴璇曪紝鏍规嵁闇瑕佹敼浣犺嚜宸辨暟鎹 E1=eye(N); %N闃跺崟浣嶇煩闃 E2=diag(ones(1,N-1),1);%涓瀵硅绾夸笂涓瀵硅绾 E2=E2+E2'; %涓诲瑙掔嚎涓婁笅鍚勪竴瀵硅绾 Ag=kron(E1,A0)+kron(E2,I0);...
绛旓細浣跨敤diag(a)鍛戒护鐢熸垚瀵硅鐭闃碉紝a涓烘煇涓悜閲忥紝濡備笅鎵绀猴細浣跨敤diag(a,i)鍛戒护鐢熸垚锛宎涓烘煇涓悜閲忥紝i涓篴鍚戦噺鐩稿涓诲瑙掔嚎鍋忕Щ鐨勫垪鏁帮紙鍚戜笂涓烘锛屽悜涓嬩负璐燂級銆傚綋i=0鏃讹紝鍙互鐩存帴鍐欐垚diag(a)銆傚叿浣撴儏鍐靛涓嬶細锛1锛塧(1:3,1) = 1 锛2锛塧 = repmat([1],3,1)锛3锛塷nes(3,1)鎴栬 锛氬厛鐢ˋ...
绛旓細n=100000;瀵硅绾夸笂鏂,瀵硅绾夸互涓嬪瑙掔嚎涓嬩笅鏂硅繘琛屾瀯閫 A=sparse([1:(n-1),1:n,3:n],[2:n,1:n,1:(n-2)],[1:(n-1),1:n,1:(n-2)],n,n);鍙互閫氳繃full(A)鏌ョ湅n杈冨皬鐨勬儏鍐 eigs(A,1,'sm')
绛旓細杈撳叆杈撳嚭 绋嬪簭鐨勮緭鍏ヨ緭鍑哄弬鏁板涓嬶細A浠ｈ〃浠绘剰涓涓渶瑕佷笁瀵硅鍖栫殑鐭╅樀锛沚鏄换鎰忎竴涓悜閲忥紝涓攂鐨勮鏁颁笌A鐨勫垪鏁扮浉鍚屽洜涓鸿鐢ㄥ埌v = A*q锛沶max鏄綘鎯宠寰楀埌鐨勭煩闃电殑澶у皬锛屼緥濡俷max=12锛屾渶鍚庡緱鍒12*12鐨涓夊瑙掔煩闃点傝緭鍑篢鏄竴涓笁瀵硅鐭╅樀銆傝皟鐢ㄥ疄渚 >> A = rand(4);>> b = rand(4,1);>> ...
绛旓細function [ x ] = Chase_method( A, b )Chase method 杩借刀娉曟眰涓夊瑙掔煩闃电殑瑙 A涓轰笁瀵硅鐭╅樀鐨勭郴鏁帮紝b涓虹瓑寮忓彸绔殑甯告暟椤癸紝杩斿洖鍊紉鍗充负鏈缁堢殑瑙 娉細A灏介噺涓烘柟闃碉紝b涓瀹氳涓哄垪鍚戦噺 姹傝拷璧舵硶鎵闇L鍙奤 T = A;for i = 2 : size(T,1)T(i,i-1) = T(i,i-1)/T(i-1,i-1);...
