HW3_solutions精选.pdfVIP

Jan. 27, 2010 Homework 3 Solutions 6.094: Introduction to Matlab Homework 3 Solutions 1-3. The following script solves problems 1-3: % shortProblems % problem 1 disp(Problem 1); % make the matrix and right hand side vector, then solve using \ mat=[3 6 4;1 5 0;0 7 7]; y=[1;2;3]; abc=mat\y; % check that the solution is correct, if it is, this should be close to 0 errorVector=mat*abc-y % problem 2 disp() disp(Problem 2); % make the x vector x=linspace(0,5,1000); approximate=trapz(x,x.*exp(-x/3)); actual=-24*exp(-5/3)+9; disp([Difference between approximate and actual integral is ... num2str(abs(approximate-actual))]); % problem 3 disp() disp(Problem 3) temp=[1 2;3 4]; i=inv(temp) identity=i*temp The script displays the following to the screen: 1 Jan. 27, 2010 Homework 3 Solutions 6.094: Introduction to Matlab 2 Jan. 27, 2010 Homework 3 Solutions 6.094: Introduction to Matlab 4. The polynomial fitting script is pasted below % polynomialFitting % load the data load randomData % make a figure and plot the data figure plot(x,y, k.,markersize,15); hold on % fit polynomials of various degrees and plot them on the figure degree=1:5; colors=jet(5); for n=1:length(degree) % fit the data, centering and scaling it [P,S,MU]=polyfit(x,y,degree(n)); % plot the data plot(x,polyval(P,x,S,MU), color,col