IT - Matlab Software (Phần 2) part 7

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ả
Nội dung trích xuất từ tài liệu:
IT - Matlab Software (Phần 2) part 7Phụ lục-Lệnh và hàm 208tæì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 209ELLIPJ(U,M) âæåüc tênh chênh xaïc âãún EPS. U vaì Mphaíi laì caïc ma tráûn cuìng cåî hoàûc mäüt trong hailaì âaûi læåüng vä hæång. Âãø chàõc chàõn khäng nháömmäâun K våïi tham säú M - chuïng quan hãû theo caïchsau: M = K^2ELLIPKEELLIPKE Têch phán elliptic âáöy âuí [K,E] = ELLIPKE(M) traí vãö giaï trë cuía têch phánelliptic âáöy âuí loaûi 1 vaì 2 trãn M. Nhæ bäø sunghiãûn thåìi, M giåïi haûn 0 < M < 1 Âäü chênh xaïc cuía ELLIPKE(M) laì EPS. Âãø chàõcchàõn khäng nháöm mäâun K våïi tham säú M - chuïng quanhãû theo caïch sau: M = K^2ERFHaìm sai säú y = erf(x) y = 2/sqrt(pi) nhán våïi têch phán tæì 0 âãún x cuíaexp(-t^2) dtERFCHaìm sai säú buì y = erfc(x) y = 2/sqrt(pi) nhán våïi têch phán tæì 0 âãún x cuíaexp(-t^2) dt = 1 - erf(x)ERFCXHaìm sai säú buì coï chia tè lãû y = erfcx(x) y = exp(x^2) * erfc(x) ~ (1/sqrt(pi)) * 1/x våïi xlåïnERFINVHaì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 GCDkhaï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*DGAMMAINC GAMMAINC Haìm gamma chæa hoaìn thaình Y = GAMMAINC(X,A) æåïc læåüng haìm gamma chæa hoaìnthaình taûi táút caí caïc pháön tæí cuía X. X phaíithæû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áûnhàòng cuìng kêch thæåïc våïi caïc âäúi säú khaïcgammainc(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ìngLCMBäü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ì BLEGENDRELEGENDRE Caïc haìm Legendre liãn âåïi P = LEGENDRE(N,X) tênh caïc haìm Legendre liãn âåïibáû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 211LOG2Taï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ãö matráû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êchthæåïc. TH tênh theo âån vë radian [TH,R,Z] = CART2POL(X,Y,Z) Biãún âäøi dæî liãûu læutrong 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êchthæåïcPOL2CART 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 tronghãû 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æutrong 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êchthæåïcSPH2CARTSPH2CART 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æutrong 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ìngkêch thæåïc. AZ vaì EL phaíi theo âån vë radianMåí vaì âoïng tãûpFOPENMåí tãûpFID = 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ãúucáö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ûomåï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âãûmNgáö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äùipermission, 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ûilæåü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 - 2004Phụ lục-Lệnh và hàm 213traí 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ãûucho 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ãûpvaì 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ãûptrong thæ muûc laìm viãûc thç FOPEN tçm xuäúng âæåìngdáùn tçm kiãúm cuía MATLAB [FID, MESSAGE] = FOPEN( filename ,permission,machineformat) måí tãûp chè âënh våïi chãú âäü chè âënhvaì âoüc dæî liãûu bàòng caïch duìng FREAD hoàûc ghidæî liãûu bàòng caïch duìng FWRITE våïi FORMAT âaî chobåíi machineformat. machineformat laì mäüt trong caïcchuäù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ìisäú liãûu 64 bit [FILE ...