※ 引述《hsuanboy (压力)》之铭言： : f(x)=0.000000000001602*x^5+0.000000003780060*x^4-0.000000034210106*x^3 : -0.000166783274301*x^2+0.000174904487267*x+1.103243573568212 : th1 : ∫ f(x)dx : 0 1 : ------------- = --------- ==> 可求出th1 : 90 100 : ∫ f(x)dx : 0 : th2 : ∫ f(x)dx : th1 1 : ------------- = --------- ==> 可求出th2 : 90 100 : ∫ f(x)dx : 0 : th3 : ∫ f(x)dx : th2 1 : ------------- = --------- ==> 可求出th3 : 90 100 : ∫ f(x)dx : 0 : 以上是我想表达的形式 第一次的积分下限初始值是0 每次求出积分上限值之後 : 会在带入到下一行的积分下限去做运算 之後又会再求出积分上限值 : 依此类推 做100次 : 下面是我自己写的程式 程式有跑出来 但是最後得到角度很怪 : 再做积分时 我是用径度去做就算 算出上限的值都为劲度 : for回圈的最後有再换成角度的形式 角度是以g表示 : 程式如下: : syms x theatai : th(1)=0; : f=0.000000000001602*x^5+0.000000003780060*x^4-0.000000034210106*x^3-0.000166783274301*x^2+0.000174904487267*x+1.103243573568212; : area=int(f,0,pi/2); : integralf=@(s,t) int(f,x,s,t); : for i=1:100 : sol = solve(integralf( th(i),theatai )-area*1/100); : d = eval(sol); : e = d(find( 0<d & d<pi/2 )); : th(i+1) = e; : g(i) = e*(180/pi); : end 针对你的函数是多项式 matlab有几个函数是针对多项式写的 在定义多项式时只要把系数降羃排成列向量 之後可以用polyval计算多项式的值 polyint可以算反导函数的系数 计算会比使用符号运算快 另外solve也可改成fsolve或fzero 避免解出多解还要挑的情况 以下是照你的程式改写给你参考 th=zeros(1,100); a=[0.000000000001602 0.000000003780060 -0.000000034210106 -0.000166783274301 0.000174904487267 1.103243573568212]'; ai=polyint(a); area=polyval(ai,pi/2)-polyval(ai,0); integralf=@(s,t) polyval(ai,t)-polyval(ai,s)-area/100; for i=1:100 sol = fzero(@(t) integralf(th(i),t),th(i)); th(i+1) = sol; g(i) = sol*(180/pi); end 角度差距是先递减再递增 因为你的函数图形是先递增再递减 可以plot(diff(g))来看 : 奇怪的问题点在 : 实际用手去计算每个角度之间的间距值应该是会越来越大才对 : 但是我算完後 每个角度的间距都是差0.9度 程式哪边有写错呢 --

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