matlab-矩阵的奇异值分解算法

%A的第三函数，sa(x)的几次求导，与w,B都无关 function[y]=disanhanshu(x,a); if(a==-1) y=sinint(x); end if(a==0) y=sin(x)/x; end if(a==1) y =cos(x)/x-sin(x)/x^2; end if(a==2) y=-sin(x)/x-2*cos(x)/x^2+2*sin(x)/x^3; end if(a==3) y=-cos(x)/x+3*sin(x)/x^2+6*cos(x)/x^3-6*sin(x)/x^4; end if(a==4) y=sin(x)/x+4*cos(x)/x^2-12*sin(x)/x^3-24*cos(x)/x^4+24*sin(x)/x^5; end if(a==5) y=cos(x)/x-5*sin(x)/x^2-20*cos(x)/x^3+60*sin(x)/x^4+120*cos(x)/x^5-120*sin(x)/x^6; end if(a==6) y=-sin(x)/x-6*cos(x)/x^2+30*sin(x)/x^3+120*cos(x)/x^4-360*sin(x)/x^5-720*cos(x)/x^6+720*sin(x)/x^7; end if(a==7) y=-cos(x)/x+7*sin(x)/x^2+42*cos(x)/x^3-210*sin(x)/x^4-840*cos(x)/x^5+2520*sin(x)/x^6+5040*cos(x)/x^7-5040*sin(x)/x^8; end if(a==8) y=sin(x)/x+8*cos(x)/x^2-56*sin(x)/x^3-336*cos(x)/x^4+1680*sin(x)/x^5+6720*cos(x)/x^6-20160*sin(x)/x^7-40320*cos(x)/x^8+40320*sin(x)/x^9; end if(a==9) y=cos(x)/x-9*sin(x)/x^2-72*cos(x)/x^3+504*sin(x)/x^4+3024*cos(x)/x^5-15120*sin(x)/x^6-60480*cos(x)/x^7+181440*sin(x)/x^8+362880*cos(x)/x^9-362880*sin(x)/x^10; end if(a==10) y=-sin(x)/x-10*cos(x)/x^2+90*sin(x)/x^3+720*cos(x)/x^4-5040*sin(x)/x^5-30240*cos(x)/x^6+151200*sin(x)/x^7+604800*cos(x)/x^8-1814400*sin(x)/x^9-3628800*cos(x)/x^10+3628800*sin(x)/x^11; end if(a==11) y=-cos(x)/x+11*sin(x)/x^2+110*cos(x)/x^3-990*sin(x)/x^4-7920*cos(x)/x^5+55440*sin(x)/x^6+332640*cos(x)/x^7-1663200*sin(x)/x^8-6652800*cos(x)/x^9+19958400*sin(x)/x^10+39916800*cos(x)/x^11-39916800*sin(x)/x^12; end if(a==12) y=sin(x)/x+12*cos(x)/x^2-132*sin(x)/x^3-1320*cos(x)/x^4+11880*sin(x)/x^5+95040*cos(x)/x^6-665280*sin(x)/x^7-3991680*cos(x)/x^8+19958400*sin(x)/x^9+79833600*cos(x)/x^10-239500800*sin(x)/x^11-479001600*cos(x)/x^12+479001600*sin(x)/x^13; end if(a==13) y=-154440*sin(x)/x^6+17160*cos(x)/x^5-1235520*cos(x)/x^7+8648640*sin(x)/x^8-259459200*sin(x)/x^10+3113510400*sin(x)/x^12-1037836800*cos(x)/x^11-156*cos(x)/x^3+1716*sin(x)/x^4+51891840*cos(x)/x^9+cos(x)/x-13*sin(x)/x^2+6227020800*cos(x)/x^13-6227020800*sin(x)/x^14; end %求sa(x)的几次求导