আমি অজগর ব্যবহার করি না, তবে যদি আমি সঠিকভাবে বুঝতে পারি তবে বাই দ্বারা
এফ( র ) =∫R0Y( x ) dএক্স
তুমি এমন কিছু ভাবছো
F=integrate(y,x)
কোথায়
F=[F1,...,Fn] গ্রিডের উপর অবিচ্ছেদ্য নমুনা দিচ্ছেন একটি ভেক্টর
x।
তবে আপনার নমুনা নেই x এবং yবরং আপনার কাছে নমুনা রয়েছে x^=log(x) এবং y^=log(y)।
অবশ্যই সহজ পদ্ধতির হবে
F=integrate(exp(y^),exp(x^)),
তবে এটি ত্রুটি-প্রবণ হবে, কারণ
y(x) মসৃণ নয়, যদিও
y^(x^) হয়।
এখন ট্র্যাপিজয়েডাল নিয়মটি আপনার ইনপুটটিকে মূলত ধরে নেয়y(x)অংশবিশেষ রৈখিক হয়। সুতরাং সাধারণ জেনারেলাইজেশনটি আপনার এটি ধরে নেওয়া উচিতy^(x^) অংশবিশেষ রৈখিক হয়।
এই ক্ষেত্রে, সংজ্ঞায়িত ΔFk=Fk+1−Fk, তোমার আছে
ΔFk=∫xk+1xky(x)dx=∫x^k+1x^key^(x^)ex^dx^=∫x^k+1x^ky~(x^)dx^
তারপরে, সংজ্ঞায়িত করা হচ্ছে t=(x^−x^k) / Δx^ট, তোমার আছে
Y^কে + টি≈y^ট+ + T ΔY^ট
এবং
Y~( T ) ≈ কইখ টি, সঙ্গে
একটি =ইY^ট+ +এক্স^ট এবং
খ = ΔY^ট+ + Δএক্স^ট।
সুতরাং অবিচ্ছেদ্য হয়
Δএফট≈ ক Δএক্স^∫10ইখ টিঘT = একটি Δএক্স^ইখ- 1খ
মতলব এ দেখতে এরকম কিছু লাগবে
dlogx=diff(logx); dlogy=diff(logy); k=1:length(logx)-1;
b=dlogx+dlogy; a=exp(logx+logy);
dF=a(k).*dlogx.*(exp(b)-1)./b;
F=cumsum([0,dF]);
আশাকরি এটা সাহায্য করবে!
(সম্পাদনা করুন: আমার উত্তরটি মূলত অনেক বেশি সংক্ষিপ্ত জবাবের সাথে একই রকম যে আমি টাইপ করার সময় দামেস্ক স্টিল দিয়েছিল। পার্থক্য কেবলমাত্র সেই ক্ষেত্রেই আমি একটি বিশেষ সমাধান দেওয়ার চেষ্টা করেছি যেখানে "নির্দিষ্ট Y( এক্স )"একটি বিভক্ত-লিনিয়ার Y^(এক্স^) একটি বিযুক্তির উপর বিচক্ষণ এক্স^ জাল, সাথে এফ(এক্স^1) = 0।)