提高Matlab 编程效率的小技巧.pdfVIP

提高Matlab 编程效率的小技巧 景明利 预分配内存 生成一个大小为n 1的向量，预先分配内存与未分配内存的对比； 未分配内存 function PreAllocmem1 n=4000; tic; for k=1:n a(k)=1; end time=toc; disp([未分配内存下赋值运算时间:,num2str(time),秒!]) 分配内存： function PreAllocmem2 n=4000; b=zeros(n,1); tic; for k=1:n b(k)=1; end time=toc; disp([预分配内存下赋值运算时间:,num2str(time),秒!]) 实验结果对比： 未分配内存下赋值运算时间:0.015792 秒! 预分配内存下赋值运算时间:4.6424e-005 秒! 2.循环中次序的问题 在多层循环中，循环次序尤为重要！ 求出一个n m 矩阵中非零元素的和 function sumat=nonzero_sum(A) [m,n]=size(A); S=0; tic for i=1:m for j=1:n if A(i,j)~=0 S=S+A(i,j); else continue end end end sumat=S; toc 循环次数多的至于内层，循环次数少的至于外层 function sumat2=nonzero_sum2(A) [m,n]=size(A); S=0; tic if n=m for i=1:m for j=1:n if A(i,j)~=0 S=S+A(i,j); else continue end end end else for j=1:n for i=1:m if A(i,j)~=0 S=S+A(i,j); else continue end end end end sumat2=S; toc 两种结果对比： A=randn(40000,1000); sumat=nonzero_sum(A) Elapsed time is 0.842068 seconds. sumat = -8.1031e+003 sumat2=nonzero_sum2(A) Elapsed time is 0.471576 seconds. sumat2 = -8.1031e+003 A=randn(40000,10); sumat=nonzero_sum(A) Elapsed time is 0.021663 seconds. sumat = 1.3607e+003 sumat2=nonzero_sum2(A) Elapsed time is 0.006993 seconds. sumat2 = 1.3607e+003 0.842068/0.471576 ans = 1.7856 0.021663/ 0.006993 ans = 3.0978 3.矩阵与向量相乘 1 2 f (x) y Φx 2 利用最速下降法求解， 2 y Φ 已知 ， ，求解 ？ x x=rand(100,1); phi=randn(100); phi=orth(ph