※ 引述《s907330204 (Bubu)》之铭言： : ※ 引述《s907330204 (Bubu)》之铭言： : : 我是要解一个数值积分问题 : : 但是方程式很多项不知道怎麽设 : : 我想到如下 但是是行不通的 : : for i=1:1000 : : y(i)=@(x) besselh(i,1,x); : : end : : y=sum(y); : : quad(y,0,pi/2); : : 问一下板上大大 : : 假如方程式很多项 要累加在一起要怎麽设计 : 我在写清楚点好了 : y=@(x) besselh(1,1,x)+besselh(2,1,x)+besselh(3,1,x)+...+besselh(1000,1,x); : 要怎麽设计 : 谢谢各位大大 真正方程式是这样 除了x以外，其他符号是已知数 Beta_Jl=@(x) exp(i*k*p_Jl*cos(x-Phi_Jl)); Beta_LI=@(x) exp(i*k*p_LI*cos(x-Phi_LI)); p %--这里方程式 SUM besselh(p,1,k*(p_lL))*exp(-i*p*(Phi_lL-x-(pi/2)))-- -p Alpha_lL=@(x) ... besselh(-16,1,k*(p_lL))*exp(-i*-16*(Phi_lL-x-(pi/2)))+... besselh(-15,1,k*(p_lL))*exp(-i*-15*(Phi_lL-x-(pi/2)))+... besselh(-14,1,k*(p_lL))*exp(-i*-14*(Phi_lL-x-(pi/2)))+... besselh(-13,1,k*(p_lL))*exp(-i*-13*(Phi_lL-x-(pi/2)))+... besselh(-12,1,k*(p_lL))*exp(-i*-12*(Phi_lL-x-(pi/2)))+... besselh(-11,1,k*(p_lL))*exp(-i*-11*(Phi_lL-x-(pi/2)))+... besselh(-10,1,k*(p_lL))*exp(-i*-10*(Phi_lL-x-(pi/2)))+... besselh(-9,1,k*(p_lL))*exp(-i*-9*(Phi_lL-x-(pi/2)))+... besselh(-8,1,k*(p_lL))*exp(-i*-8*(Phi_lL-x-(pi/2)))+... besselh(-7,1,k*(p_lL))*exp(-i*-7*(Phi_lL-x-(pi/2)))+... besselh(-6,1,k*(p_lL))*exp(-i*-6*(Phi_lL-x-(pi/2)))+... besselh(-5,1,k*(p_lL))*exp(-i*-5*(Phi_lL-x-(pi/2)))+... besselh(-4,1,k*(p_lL))*exp(-i*-4*(Phi_lL-x-(pi/2)))+... besselh(-3,1,k*(p_lL))*exp(-i*-3*(Phi_lL-x-(pi/2)))+... besselh(-2,1,k*(p_lL))*exp(-i*-2*(Phi_lL-x-(pi/2)))+... besselh(-1,1,k*(p_lL))*exp(-i*-1*(Phi_lL-x-(pi/2)))+... besselh(0,1,k*(p_lL))*exp(-i*0*(Phi_lL-x-(pi/2)))+... besselh(1,1,k*(p_lL))*exp(-i*1*(Phi_lL-x-(pi/2)))+... besselh(2,1,k*(p_lL))*exp(-i*2*(Phi_lL-x-(pi/2)))+... besselh(3,1,k*(p_lL))*exp(-i*3*(Phi_lL-x-(pi/2)))+... besselh(4,1,k*(p_lL))*exp(-i*4*(Phi_lL-x-(pi/2)))+... besselh(5,1,k*(p_lL))*exp(-i*5*(Phi_lL-x-(pi/2)))+... besselh(6,1,k*(p_lL))*exp(-i*6*(Phi_lL-x-(pi/2)))+... besselh(7,1,k*(p_lL))*exp(-i*7*(Phi_lL-x-(pi/2)))+... besselh(8,1,k*(p_lL))*exp(-i*8*(Phi_lL-x-(pi/2)))+... besselh(16,1,k*(p_lL))*exp(-i*16*(Phi_lL-x-(pi/2)))+... besselh(15,1,k*(p_lL))*exp(-i*15*(Phi_lL-x-(pi/2)))+... besselh(14,1,k*(p_lL))*exp(-i*14*(Phi_lL-x-(pi/2)))+... besselh(13,1,k*(p_lL))*exp(-i*13*(Phi_lL-x-(pi/2)))+... besselh(12,1,k*(p_lL))*exp(-i*12*(Phi_lL-x-(pi/2)))+... besselh(11,1,k*(p_lL))*exp(-i*11*(Phi_lL-x-(pi/2)))+... besselh(10,1,k*(p_lL))*exp(-i*10*(Phi_lL-x-(pi/2)))+... besselh(9,1,k*(p_lL))*exp(-i*9*(Phi_lL-x-(pi/2))); %----------------------------------------------------------------- f=@(x) 1/(2*pi).*(Beta_Jl(x).*Alpha_lL(x).*Beta_LI(x)); H_besselh=quad(f,0,2*pi) 以上p不等於16，应该大於16，那是我硬干的写法 所以想说，各位高手，有没有设计方法可以解决很多项累加问题 谢谢各位大大。 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 114.34.20.7