আর এর সাথে বায়ু সম্পর্কিত তথ্য বিশ্লেষণ করছে

Question

হাই, আমি একটি বায়ু টারবাইন থেকে শক্তি অনুমানের জন্য বায়ু ডেটা অ্যানালাইসিং করছি।
আমি 10 বছরের বায়ু সংক্রান্ত ডেটা নিয়েছি এবং একটি হিস্টোগ্রাম গ্রাফ করেছি;
আমার দ্বিতীয় পর্যায়ে ছিল ডেটাতে একটি ওয়েইবুল বিতরণ।
আমি lmomওয়েবেল আকৃতিটি গণনা করতে প্যাকেজটির সাথে আর ব্যবহার করেছি এবং এটি ব্যবহার করা কোডটি হ'ল:

>library(lmom)    
wind.moments<-samlmu(as.numeric(pp$WS))      
moments<-pelwei(wind.moments)     
x.wei<-rweibull(n=length(pp$WS), shape=moments["delta"], scale=moments["beta"])    
hist(as.numeric(pp$WS), freq=FALSE)    
lines(density(x.wei), col="red", lwd=4)    

দেখে মনে হচ্ছে ডেটা এবং ঘনত্বের কার্যের মধ্যে কিছুটা পিছিয়ে আছে; আপনি আমাকে এই ব্যাপারে সাহায্য করতে পারেন? আর একটি প্রশ্ন আপনি ঘনত্ব ফাংশন থেকে অ্যানুয়াল শক্তি গণনা করতে আমাকে সাহায্য করতে পারেন?

ধন্যবাদ

Comments

ছবিটি সম্পর্কে, কিছু চিত্র হোস্টিংয়ে পোস্ট করুন এবং একটি লিঙ্ক রাখুন - আমি এটিকে একটি আটকানো ছবিতে রূপান্তর করব।

+1, আকর্ষণীয় প্রশ্ন, আপনি খুব শীঘ্রই আপনার যথেষ্ট খ্যাতি অর্জন করতে পারেন :)

— এমপিটকাস

গ্রাফ থেকে বিচার করে, সমস্যাটি পিছিয়ে নেই। আপনি যা চক্রান্ত করেছেন তা মোটামুটি মাপসই ness সুতরাং দেখে মনে হচ্ছে যে ওয়েইবুল বিতরণ আপনার ডেটার জন্য উপযুক্ত নয়। আমি দেখতে পাচ্ছি যে শূন্যের কাছাকাছি একটি গুচ্ছ রয়েছে, আপনার ডেটাতে কি শূন্য মান আছে? সেক্ষেত্রে আপনার আলাদাভাবে শূন্য মানগুলি মডেল করতে হবে। সুতরাং প্রথম পরামর্শটি হ'ল শূন্য অমূল্যের জন্য ওয়েবুলের চেষ্টা করা। এছাড়াও কেন ওয়েইবুল, একই কারণের কোনও নির্দিষ্ট রেফারেন্স রয়েছে?

— এমপিটকাস

মনে রাখবেন যে 'ল্যাগ' একটি শব্দ যা মূলত সময়ের সাথে সাথে বিশ্লেষণে ব্যবহৃত হয়, একের পর এক ঘটতে থাকে। এটি কোনও ল্যাগ নয় - সম্ভবত এটি আরও সঠিকভাবে একটি শিফট বলা হয় - বা একটি অফসেট - তবে শিফট সম্ভবত বিতরণগুলির জন্য সাধারণ, তারা শিফ্ট এবং স্কেল করে।

— স্পেসডম্যান

as.numeric (x) ফ্যাক্টর সহ ব্যবহার সম্পর্কে সতর্কতা অবলম্বন করুন; আপনি ফ্যাক্টরের জন্য সঠিক সংখ্যার মান পেয়েছেন তা নিশ্চিত করার জন্য আপনি প্রকৃতপক্ষে as.numeric (as.character (x)) ব্যবহার করতে চান।

— অ্যান্ডি ক্লিফটন

Answer 1

lmomফাংশন , অবস্থান , স্কেল এবং আকারের পরামিতি সহ তিনটি প্যারামিটার ওয়েইবুল বিতরণ pelweiফিট করে । দ্বি- প্যারামিটার ওয়েইবুল বিতরণের জন্য এলোমেলো সংখ্যা তৈরি করে । আপনাকে অবস্থানের প্যারামিটারটি বিয়োগ করতে হবে । এটি আরও ভাল ফিট করা উচিত , তবে এটি আপনার নির্দিষ্ট ডেটাতে ভাল ফিট দেবে বলে মনে হয় না ।rweibullmoments["zeta"]

আমি লক্ষ্য করি http://www.reuk.co.uk/Wind-Speed-Dist वितरण- Weibull.htm বলছে "বিশ্বের বেশিরভাগ অংশে উইন্ডুল ডিস্ট্রিবিউশন ব্যবহার করে বাতাসের গতি মডেল করা যায়।" সম্ভবত আপনি কেবল দুর্ভাগ্য এবং বিশ্বের এমন একটি অংশে বাস করছেন যেখানে তারা পারবেন না!

Comments

আপনি কি আমার দ্বিতীয় প্রশ্নটি সম্পর্কে আমাকে সহায়তা করতে পারেন: ঘনত্বের ফাংশনের ভিত্তিতে বায়ু শক্তির গণনা করছেন?

— eliavs

আমি পদার্থবিজ্ঞানী নই এবং আমি প্রয়োজনীয় সমীকরণগুলি জানি না, তবে আমি ধারণা করি এটি ঘনত্বের সাথে একটি সংখ্যাসমূহকে জড়িত করবে। আর এর integrate()কাজটি তার জন্য কার্যকর হতে পারে।

— onestop

আমি জানি আমার প্রোব্লাম সমীকরণটি কি আমি বায়ুর প্রতিটি গতিবেগের কত শতাংশ সময় গণনা করতে চাই

— eliavs

আমি যা করছি তা আপনি এই সংহতকরণে আমাকে সহায়তা করতে পারেন -> আপনাকে ধন্যবাদ

— এলিয়াভস

Answer 2

আমি আপনার প্লটটি http://hawaii.gov/dbedt/ert/winddata/krab0192.txt থেকে ডেটা দিয়ে পুনরায় তৈরি করেছি (আমি 1200 পরিমাপ নিয়েছি)। আমি সাধারণত আপনার কোড ব্যবহার করে ডেটার একটি শালীন ফিট পেয়েছি:

library(lmom)

daten <- read.delim("wind.txt")
wind.avg <- na.omit(as.numeric(daten[,"X12"]))
wind.moments<-samlmu(wind.avg)
moments<-pelwei(wind.moments)
x.wei<-rweibull(n=length(wind.avg), shape=moments["delta"], scale=moments["beta"])
hist(as.numeric(wind.avg), freq=FALSE)
lines(density(x.wei), col="red", lwd=4)

দুঃখিত, আমি সমস্যা হইনি যদি আপনার সমস্যা হতে পারে তবে আমি মনে করি আপনার ডেটাতে ওয়েলবুল ফিট করতে পারা উচিত। যা আমাকে সন্দেহজনক করে তোলে তা হ'ল আপনার ঘনত্বের চক্রান্তের বেল-বক্রতা, আমি জানি না যে এটি কোথা থেকে এসেছে।

আমি উত্পন্ন মুহুর্তগুলি এখানে:

wind.moments

       l_1         l_2         t_3         t_4 
15.17287544  4.80372580  0.14963501  0.06954438

মুহূর্ত

     zeta      beta     delta 
 0.516201 16.454233  1.745413 

বার্ষিক আউটপুটে ডব্লিউটিআর: আমি মনে করি আমি সম্ভাবনা ঘনত্ব ফাংশনের জন্য পৃথক মান উত্পন্ন করতাম, আউটপুট ফাংশনটির সাথে এই মানগুলিকে গুণিত করতাম এবং সমষ্টি করতাম। বিকল্পভাবে, আপনি কেবলমাত্র আপনার কাঁচা ডেটা ব্যবহার করতে পারেন, আউটপুট ফাংশনটির সাথে মানগুলি গুন করতে পারেন, এটির যোগফল তৈরি করতে এবং বার্ষিক গড় গণনা করতে পারেন, আপনাকে উপযুক্ত উপায়ে মৌসুমীকরণের জন্য নিয়ন্ত্রণ করা উচিত (যেমন পুরো বছর ব্যবহার করা নিশ্চিত করুন বা তদনুসারে ওজন অনুযায়ী) ।

এখানে অনিয়ন্ত্রিত আউটপুট ( http://www.articlesbase.com/diy-articles/determining-wind-turbine-annual-power-output-a-simple-forula-based-upon-blade-diype- এর সূত্র ব্যবহার করে ) এবং-আপনার-অবস্থান-এ-গড়-বায়ু-গতি-এ-অবস্থান -1303080 html )

years  <- length(wind.avg)/365
diameter <- 150
Power = (0.01328*diameter^2)*((wind.avg)^3)
(annual.power <- sum(Power)/years)
[1] 791828306

Answer 3

এখানে উইন্ডো টারবাইনগুলির বিষয়ে এসওতে একটি সাম্প্রতিক পোস্ট। এই লিঙ্কটিতে আমার উত্তরটির তিনটি লিঙ্ক রয়েছে যাতে আপনি আগ্রহী হতে পারেন:

/programming/4843194/r-language-sorting-data-into-ranges-averaging-ignore-outliers/4844783#4844783

আমি উপরের এসও উত্তরে কেবলমাত্র ওয়েইবুল লিঙ্কগুলির মধ্যে একটি চেক করেছি। কিছু কারণে লিঙ্কটি নিচে রয়েছে। এখানে কয়েকটি লিঙ্ক রয়েছে যা একই বেসিক তথ্য সরবরাহ করে:

http://www.gso.uri.edu/ozone/

http://www.weru.ksu.edu/new_weru/publications/pdf/Comparison%20of%20the%20Weibull%20model%20with%20measured%20wind%20speed%20distributions%20for%20stochastic%20wind%20genera.pdf

http://www.kfupm.edu.sa/ri/publication/cer/41_JP_Weibull_parameters_for_wind_speed_distribution_in_Saudi_Arabia.pdf

http://journal.dogus.edu.tr/13026739/2008/cilt9/sayi1/M00195.pdf

http://www.eurojournals.com/ejsr_26_1_01.pdf

এছাড়াও, বাতাস থেকে উত্পন্ন বিদ্যুৎ থেকে, seasonতুরতা সুস্পষ্ট।

Answer 4

আমি নিশ্চিত না যে কেউ ইতিমধ্যে এই পয়েন্টটি তৈরি করেছে কিনা তবে পেলওয়েই একটি নির্দিষ্ট বাউন্ডে যোগ করে 2 প্যারামিটার ওয়েবেল ফাংশন হিসাবে কাজ করতে বাধ্য হতে পারে।

কল করার জন্য moments<-pelwei(wind.moments)আপনাকে কেবল কল করা উচিতmoments<-pelwei(wind.moments,bound=0)

জিটা মানটি কী তা আপনি সর্বদা যাচাই করতে পারেন। যদি এটি 0 না হয় এবং আপনি ডুইবুল ব্যবহার করছেন তবে আপনার এটি সম্পর্কে কিছু করা দরকার।