数值方法作业一解答.docVIP

數值方法作業一 解答 (a) X=[-1000:10:1000] (b) X=[0:0.01:1] (c) X=[1:100] (d) X=[ones(1,50) zeros(1,50)] (e) X=[diag( 1 : 100 )] (f) X=1:100; repmat(X,100,1) (g) X = eye(3 , 4)\ (h) rand(10,20)*10-5 (i) X= randn(10, 20) (j) 1. A = 1 : 100 ; X = repmat(A,100,1) ; repmat(X,10,1) 2. X = repmat(A,100,1) ; repmat(X,1,10) (k) A= [1 2;1 2] ; X= repmat(A,10,10) (2) (a) diag(1:2:7) 創造一個對角矩陣，對角元素是從1每隔2到7 其餘為0。 (b) [以從1到7,且間隔為2的元素，置於每個對角線元素下方，其餘元素為0]+[以從1到9,且間隔為2的向量當做對角線,其餘元素均為0的矩陣]+[以從1到7,且間隔為2的向量,置於每個對角線元素右方,其餘元素 均為0的矩陣] (c) a =b所以此時a已經等於b的值3,b又要等於a,因為a已經是3所以b的值還是等於3,so a+b=6。 (d) x是1到5的數字，y是5到1的數字,現在要求的是z = x y， 如果結果正確則顯示為1，錯誤則顯示為0。 z =0 0 0 1 1 (e) 造一個線性方程式為x^2-1，然後把x = 1和2的結果帶入 (f) f=inline(x^2 - 1);f([1 2]) ??? Error using == inlineeval Error in inline expression == x^2 - 1 ??? Error using == mpower Matrix must be square. 與上一題類似，但會出現錯誤，因為X 的平方項表示方式錯誤，而且輸入的矩陣需要是方陣才能執行 (g) 同上 令a為[1 2 3 4 5 6]的矩陣，接著讓其以2 x 3的形式重新排列，第一行為前兩個元素，第二行為皆下來的兩個元素，以此類推 (i)把a矩陣先轉置再重新排為2x3的矩陣 (3) (a) x = [-1 1 -1 1 -1 -1 1 1] ; y = [1 1 -1 -1 1 -1 -1 1] ; plot(x,y,.) ; hold on ; plot(x,y) (b) x=2*rand(1,800)-1; y=2*rand(1,800)-1; plot(x,y,.); (c) x=linspace( -2*pi , 2*pi); y=cos(x) + sin(x); plot(x,y) (d) x=linspace(-2*pi,2*pi); subplot(2,2,1); plot(x,cos(x)); subplot(2,2,2); plot(x,cos(2*x)); subplot(2,2,3); plot(x,cos(3*x)); subplot(2,2,4); plot(x,cos( 4*x)); (4) i=1:20; y=i.^4; sum(y)=722666