IT - Matlab Software (Phần 2) part 7

Chia sẻ: Fewgnmerihnweil Bgmrtlihnmeilbni | Ngày: | Loại File: PDF | Số trang:8

Thêm vào BST

Báo xấu

92
lượt xem 6
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Tham khảo tài liệu 'it - matlab software (phần 2) part 7', công nghệ thông tin, kỹ thuật lập trình phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: IT - Matlab Software (Phần 2) part 7

Phụ lục-Lệnh và hàm 208 tæìng âån vë, nghéa laì ALPHA = alpha:1:alpha+n-1, vaì caïc giaï trë phaíi thoaí maîn 0
Phụ lục-Lệnh và hàm 209 ELLIPJ(U,M) âæåüc tênh chênh xaïc âãún EPS. U vaì M phaíi laì caïc ma tráûn cuìng cåî hoàûc mäüt trong hai laì âaûi læåüng vä hæång. Âãø chàõc chàõn khäng nháöm mäâun K våïi tham säú M - chuïng quan hãû theo caïch sau: M = K^2 ELLIPKE ELLIPKE Têch phán elliptic âáöy âuí [K,E] = ELLIPKE(M) traí vãö giaï trë cuía têch phán elliptic âáöy âuí loaûi 1 vaì 2 trãn M. Nhæ bäø sung hiãûn thåìi, M giåïi haûn 0 < M < 1 Âäü chênh xaïc cuía ELLIPKE(M) laì EPS. Âãø chàõc chàõn khäng nháöm mäâun K våïi tham säú M - chuïng quan hãû theo caïch sau: M = K^2 ERF Haìm sai säú y = erf(x) y = 2/sqrt(pi) nhán våïi têch phán tæì 0 âãún x cuía exp(-t^2) dt ERFC Haìm sai säú buì y = erfc(x) y = 2/sqrt(pi) nhán våïi têch phán tæì 0 âãún x cuía exp(-t^2) dt = 1 - erf(x) ERFCX Haìm sai säú buì coï chia tè lãû y = erfcx(x) y = exp(x^2) * erfc(x) ~ (1/sqrt(pi)) * 1/x våïi x låïn ERFINV Haìm ngæåüc cuía haìm sai säú x = erfinv(y) thoía maîn y = erf(x), -1
Phụ lục-Lệnh và hàm 210 G = GCD(A,B) laì æåïc chung låïn nháút cuía caïc säú nguyãn A vaì B G = GCD(0,0) = 0 theo quy æåïc; táút caí caïc GCD khaïc âãöu nguyãn dæång [G,C,D] = GCD(A,B) cuîng traí vãö C vaì D våïi G = A*C + B*D GAMMAINC GAMMAINC Haìm gamma chæa hoaìn thaình Y = GAMMAINC(X,A) æåïc læåüng haìm gamma chæa hoaìn thaình taûi táút caí caïc pháön tæí cuía X. X phaíi thæûc. X vaì A phaíi cuìng kêch thæåïc, ngoaûi træì haìm coï caïc âäúi säú vä hæåïng nhæ caïc ma tráûn hàòng cuìng kêch thæåïc våïi caïc âäúi säú khaïc gammainc(x,a) = (têch phán tæì 0 âãún x cuía t^(a-1) exp(-t) dt)/gamma(a) Læu yï: gammainc(x,a) tiãún vãö 1 khi x tiãún vãö vä cuìng LCM Bäüi chung nhoí nháút LCM(A,B) laì bäüi chung nhoí nháút cuía caïc säú nguyãn dæång A vaì B LEGENDRE LEGENDRE Caïc haìm Legendre liãn âåïi P = LEGENDRE(N,X) tênh caïc haìm Legendre liãn âåïi báûc N vaì theo thæï tæû M = 0, 1, ..., N, tênh trãn X. N laì mäüt säú nguyãn nhoí hån 257. X laì vectå coï caïc pháön tæí X(j) thæûc thoía maîn abs(X(j))
Phụ lục-Lệnh và hàm 211 LOG2 Taïch caïc säú thæûc. Theo chuáøn IEEE laì logb() y = LOG2(x) tênh lä-ga-rit cå säú 2 cuía x [f,e] = LOG2(x) våïi ma tráûn thæûc x, traí vãö ma tráûn thæûc, thæåìng trong miãön 0.5
Phụ lục-Lệnh và hàm 212 [M,N] = SIZE(X), thç Y cuîng phaíi cuìng kêch thæåïc. TH tênh theo âån vë radian [TH,R,Z] = CART2POL(X,Y,Z) Biãún âäøi dæî liãûu læu trong hãû toüa âäü Âãö-caïc sang hãû toüa âäü truû. Nãúu [M,N] = SIZE(X), thç Y vaì Z phaíi cuìng kêch thæåïc POL2CART POL2CART Biãún âäøi hãû toüa âäü cæûc sang hãû toüa âäü Âãö-caïc [X,Y] = POL2CART(TH,R) Biãún âäøi dæî liãûu læu trong hãû toüa âäü cæûc sang hãû toüa âäü Âãö-caïc. Nãúu [M,N] = SIZE(TH), thç R cuîng phaíi cuìng kêch thæåïc. TH phaíi theo âån vë radian . [X,Y,Z] = POL2CART(TH,R,Z) Biãún âäøi dæî liãûu læu trong hãû toüa âäü truû sang hãû toüa âäü Âãö-caïc. Nãúu [M,N] = SIZE(TH), thç R vaì Z phaíi cuìng kêch thæåïc SPH2CART SPH2CART Biãún âäøi hãû toüa âäü cáöu sang hãû toüa âäü Âãö-caïc [X,Y,Z] = SPH2CART(AZ,EL,R) Biãún âäøi dæî liãûu læu trong hãû toüa âäü cáöu sang hãû toüa âäü Âãö-caïc. Nãúu [M,N] = SIZE(AZ), thç EL vaì R cuîng phaíi cuìng kêch thæåïc. AZ vaì EL phaíi theo âån vë radian Måí vaì âoïng tãûp FOPEN Måí tãûp FID = FOPEN('filename',permission) måí tãûp coï tãn chè âënh laì filename våïi chãú âäü cho pheïp chè âënh laì permission. Permission laì mäüt trong caïc chuäùi: 'r'_ âoüc 'w'_ ghi (taûo måïi nãúu cáön) 'a'_ näúi thãm vaìo cuäúi tãûp (taûo måïi nãúu cáön) 'r+' âoüc vaì ghi (khäng taûo måïi) 'w+' càõt tãûp hoàûc taûo måïi âãø âoüc vaì ghi 'a+' âoüc vaì ghi thãm vaìo cuäúi tãûp (taûo måïi nãúu cáön) 'W' ghi maì khäng tæû âäüng xoïa vuìng âãûm 'A' näúi thãm maì khäng tæû âäüng xoïa vuìng âãûm Ngáöm âënh caïc tãûp âæåüc måí trong chãú âäü nhë phán. Âãø måí tãûp vàn baín thç thãm 't' vaìo chuäùi permission, nhæ 'rt' vaì 'wt+' FID = FOPEN(' filename') giaí sæí permission laì 'r'. Nãúu måí thaình cäng thç danh hiãûu FID láúy mäüt âaûi læåüng vä hæåïng nguyãn cuía MATLAB, danh hiãûu tãûp âæåüc duìng nhæ âäúi säú âáöu tiãn cho caïc phuûc vuû vaìo ra tãûp. Nãúu viãûc måí tãûp khäng thaình cäng thç Phan Thanh Tao - 2004
Phụ lục-Lệnh và hàm 213 traí vãö -1 cho FID. Ba danh hiãûu tãûp âæåüc tæû âäüng duìng maì khäng cáön måí. Chuïng laì fid=0 (nháûp chuáøn), fid=1 (xuáút chuáøn), vaì fid=2 (thiãút bë läùi chuáøn) [FID, MESSAGE] = FOPEN(' filename ',permission) traí vãö mäüt läùi hãû thäúng nãúu måí tháút baûi. FOPEN('all') traí vãö mäüt vectå doìng, caïc danh hiãûu cho táút caí caïc tãûp âang måí båíi ngæåìi duìng. (Nhæng khäng phaíi 0, 1, vaì 2) [FILENAME, PERMISSION] = FOPEN(FID) traí vãö tãn tãûp vaì chãú âäü våïi danh hiãûu tãûp âaî cho. Nãúu tãûp âæåüc måí trong chãú âäü 'r' vaì khäng tçm tháúy tãûp trong thæ muûc laìm viãûc thç FOPEN tçm xuäúng âæåìng dáùn tçm kiãúm cuía MATLAB [FID, MESSAGE] = FOPEN(' filename ',permission, machineformat) måí tãûp chè âënh våïi chãú âäü chè âënh vaì âoüc dæî liãûu bàòng caïch duìng FREAD hoàûc ghi dæî liãûu bàòng caïch duìng FWRITE våïi FORMAT âaî cho båíi machineformat. machineformat laì mäüt trong caïc chuäùi sau: 'native' hoàûc 'n' - FORMAT maïy cuûc bäü- ngáöm âënh 'ieee-le' hoàûc 'l' - chuáøn IEEE 'ieee-be' hoàûc 'b' - chuáøn IEEE 'vaxd' hoàûc 'd' - chuáøn VAX D 'vaxg' hoàûc 'g' - chuáøn VAX G 'cray'_ hoàûc 'c' - chuáøn Cray 'ieee-le.l64' hoàûc 'a' - chuáøn IEEE vaì âäü daìi säú liãûu 64 bit [FILENAME,PERMISSION,MACHINEFORMAT] = FOPEN(FID) traí vãö tãn tãûp, chãú âäü måí, vaì daûng maïy våïi daûnh hiãûu tãûp. Chãú âäü 'W' vaì 'A' âæåüc thiãút kãú âãø duìng våïi caïc âéa tæì vaì khäng tæû âäüng xoïa vuìng âãûm hiãûn thåìi sau caïc thao taïc xuáút. Vê duû, måí 1/4" bàng tæì trãn traûm SPARC âãø ghi maì khäng tæû âäüng xoïa vuìng âãûm: fid = fopen('/dev/rst0','W') FCLOSE Âoïng tãûp FCLOSE(FID) âoïng tãûp våïi danh hiãûu FID, laì mäüt säú nguyãn nháûn âæåüc tæì FOPEN() træåïc âáy FCLOSE('all') âoïng táút caí caï tãûp âang måí træì 0, 1 vaì 2 FCLOSE() traí vãö s 0 nãúu thaình cäng, -1 nãúu ngæåüc laûi Phan Thanh Tao - 2004
Phụ lục-Lệnh và hàm 214 Vaìo/Ra tãûp khäng daûng thæïc FREAD Âoüc tãûp nhë phán [A, COUNT] = FREAD(FID,SIZE,PRECISION) âoüc dæî liãûu nhë phán tæì tãûp chè âënh vaì ghi vaìo ma tráûn A. Âäúi säú xuáút tuìy choün COUNT traí vãö säú pháön tæí âoüc thaình cäng. FID laì danh hiãûu nguyãn cuía tãûp âæåüc måí tæì FOPEN. Âäúi säú tuìy choün SIZE; nãúu khäng coï thç âoüc toaìn bäü tãûp; nãúu coï thç caïc âãö muûc håüp lyï laì: N âoüc N pháön tæí vaìo mäüt vectå cäüt. inf âoüc vaìo cuäúi tãûp. [M,N] âoüc caïc pháön tæí âãø láúp âáöy ma tráûn cåî MxN, theo thæï tæû cäüt N coï thãø laì inf, nhæng khäng laì M Âäúi säú PRECISION âiãöu khiãøn säú bit âoüc cho mäùi giaï trë vaì thäng dëch caïc bit naìy nhæ caïc giaï trë kyï tæû, säú nguyãn hoàûc säú thæûc. Caïc chuäùi sau coï thãøø duìng, hoàûc caïc phiãn baín MATLAB, hoàûc C hoàûc Fortran tæång âæång. Nãúu khäng chè âënh thç âæåüc ngáöm âënh laì 'uchar'. Caïc pháön tæí cuía ma tráûn kãút quaí luän læu trong daûng säú thæûc daìi cuía MATLAB MATLAB C hoàûc Fortran Mä taí 'char' 'char' kyï tæû , 8 bit 'schar' 'signed char' kyï tæû coï dáúu, 8 bit 'short' 'short' säú nguyãn, 16 bit 'int' 'int' säú nguyãn, 16 hoàûc 32 bit 'long' 'long' säú nguyãn, 32 hoàûc 64 bit* 'float' 'float' säú thæûc , 32 bit 'double' 'double' säú thæûc daìi , 64 bit 'uchar' 'unsigned char' kyï tæû khäng dáúu, 8 bit 'ushort' 'unsigned short' säú nguyãn khäng dáúu, 16 bit 'uint' 'unsigned int' säú nguyãn khäng dáúu, 16 hoàûc 32 bit 'ulong' 'unsigned long' säú nguyãn khäng dáúu, 32 bit 'char' 'char*1' kyï tæû , 8 bit 'float32' 'real*4' säú thæûc 32 bit 'float64' 'real*8' säú thæûc 64 bit 'int8' 'integer*1' säú nguyãn, 8 bit Phan Thanh Tao - 2004
Phụ lục-Lệnh và hàm 215 'int16' 'integer*2' säú nguyãn, 16 bit 'int32' 'integer*4' säú nguyãn, 32 bit 'integer*8' säú nguyãn, 64 bits** 'intN' säú nguyãn coï dáúu, N bit räüng 'uintN' säú nguyãn khäng dáúu, N bit räüng N biãøu hiãûn giaï trë báút kyì giæîa 1 vaì 32 * 64 bits trãn DEC alpha ** chè duìng trãn DEC alpha Vê duû, fid = fopen('FREAD.m','r'); F = FREAD(fid); s = setstr(F') måí tãûp chæïa âãö muûc HELP naìy, räöi âoüc vaì hiãøn thë toaìn bäü tãûp, ngáöm âënh duìng size = inf vaì precision = 'uchar'. Kãút quaí length(F) laì säú kyï tæû trong tãûp [A, COUNT] = FREAD(FID,SIZE,PRECISION,SKIP) âæa vaìo âäúi säú tuìy choün SKIP âãø chè âënh säú byte boí qua sau mäùi láön âoüc. Caïch naìy thæåìng duìng âãø trêch dæî liãûu trong caïc træåìng khäng liãn tuûc tæì caïc baín ghi coï âäü daìi cäú âënh FWRITE Ghi dæî liãûu vaìo mäüt tãûp nhë phán COUNT = FWRITE(FID,A,PRECISION) ghi caïc pháön tæí cuía ma tráûn A vaìo tãûp chè âënh, dëch caïc giaï trë MATLAB sang âäü chênh xaïc chè âënh laì precision. Dæî liãûu âæåüc ghi theo thæï tæû cäüt. COUNT laì säú pháön tæí âæåüc ghi thaình cäng. FID laì danh hiãûu säú nguyãn cuía tãûp nháûn tæì FOPEN, hoàûc 1 chi xuáút chuáøn hoàûc 2 cho läùi chuáøn. PRECISION âiãöu khiãøn daûng vaì kêch thæåïc cuía kãút quaí. Xem danh saïch PRECISION cho pheïp trong FREAD COUNT = FWRITE(FID,A,PRECISION,SKIP) âæa vaìo âäúi säú tuìy choün SKIP âãø chè âënh säú byte boí qua træåïc mäùi láön ghi. Caïch naìy thæåìng duìng âãø cheìn dæî liãûu trong caïc træåìng khäng liãn tuûc tæì caïc baín ghi coï âäü daìi cäú âënh Vê duû fid = fopen('magic5.bin','wb') FWRITE(fid,magic(5),'integer *4') Taûo måïi mäüt tãûp nhë phán 100-byte, chæïa 25 pháön tæí cuía ma phæång báûc 5, læu nhæ caïc säú nguyãn 4- byte Nháûp/xuáút tãûp coï daûng thæïc FSCANF Âoüc dæî liãûu coï daûng thæïc trong mäüt tãûp Phan Thanh Tao - 2004