MATLAB简介3MATLABm—file程式与多项式函数.pptVIP

而将两个多项式相乘可以得到一新的多项式 e(x) =? a(x) b(x) ????????????????????????? 如果是两个多项式相除，即 ?????????????????????： ????????????????????????? 上述二个运算式不能直接运算，须要另外定义函数conv做乘法运算以及函数deconv做除法运算。当二多项式相乘，在数学上等于二个阵列做卷积(convolution)运算（因为我们是以阵列来代表一个多项式的各阶系数）， 因此可利用conv函数做乘法运算，其语法为conv(a,b)，其中a, b代表二个多项式的阵列。而二多项式相除就相 当于反卷积(de-convolution) 运算，因此有 deconv 函数，其语法稍有不同 [q,r]=deconv(a,b)，其中q,r分别代表整 除多项式及余数多项式。如果直接使用q=deconv(a,b)则只求出整除多项式。 以下就介绍相关范例，来说明二个多项式的乘除运算： a=[1 2 3 4]; b=[1 4 9 16]; e=conv(a,b) e = 1 6 20 50 75 84 64 g=e+[0 0 0 c] g = 1 6 20 52 81 96 84 (c = [ 2 6 12 20]) [f,r]=deconv(e,b) f = 1 2 3 4 r = 0 0 0 0 0 0 0 % 因为是整除所以余数多项式的各系数皆为零 [h,r]=deconv(g,a) h = 1 4 9 18 r = 0 0 0 0 2 6 12 % 余数多项式为 2*x^2 + 6*x + 12 多项式的根 一个多项式视其阶数而定，它的根可以有一个到数个，可能为实数也可能是复数。要求一高阶多项式的根往 往须借助数值方法，所幸MATLAB已将这些数值方法写成一函数roots(p)，我们只要输入多项式的各阶系数（ 以 p 代表）即可求解到对应的根 p=[1 3 2]; r=roots(p) r = -2 -1 p=[1 -12 0 25 116]; % 注意二阶项系数为零须要输入，否则多项式的阶数就不对 r=roots(p) % 有实数根及复数根 r = 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i * What is an m-file? An m-file is a file with extension .m It can be used to store commands Matlab reads and executes the commands in an m-file if you type the name of the file in the command window Use an ordinary text editor and save the file as filename.m Choose ’new file’ under the File Menu Different ways of creating m-files In the first case the function inv(A) is used to find the inverse. This is then multiplied by b. The Matlab code for this operation is In the second case left- division is performed straight away with the command X= 9.2500 4.2500 2.7500 X= 9.2500 4.2500 2.7500 Example Example solvetime.m A = rand(1000,1000); % Creates a random % matrix A b = rand(1000,1); % Creates a random % vector b det(A) % Calculates the determinant of A tic, % Starts the time-watch x = inv(A)*b; % Solves the system toc