signals - Comparing FFT of Function to Analytical FT Solution in Matlab -

-

मैं फ़ंक्शन के विश्लेषणात्मक फूरियर ट्रांस्फ़ॉर्म, एक्सपी (- (- ( W ^ 2) / 4) / sqrt (2), आवृत्ति रेंज -3 से 3 पर।

मैंने निम्न मैटलैब कोड लिखा है और इसे सफलतापूर्वक बिना कई बार फिर से चालू किया है। 

  fs = 100; % नमूना आवृत्ति डीटी = 1 / एफएस; टी = 0: डीटी: 10-डीटी; % समय वेक्टर एल = लंबाई (टी); नमूना बिंदुओं की% संख्या% N = 2 ^ अगली पोज़ 2 (एल); %ज़रूरी? Y = एपीपी (- (टी। ^ 2)); वाई = डीटी * ifftshift (पेट (fft (y))); Freq = (-L / 2: एल / 2-1) * एफएस / एल; % Freq vector F = (exp (- (freq। ^ 2) / 4)) / sqrt (2); % विश्लेषणात्मक समाधान% Y_valid_pts = वाई (डब्ल्यू & gt; = - 3 & amp; डब्ल्यू & lt; = 3); % Freq के लिए तुलना = -3 से 3% npts = लंबाई (Y_valid_pts); % W = linspace (-3,3, npts); % फ़े = (एक्सपी (- (डब्ल्यू। ^ 2) / 4)) / एसक्यूआरटी (2); त्रुटि = आदर्श (वाई - एफ) त्रुटि के लिए% L2 मानक; साजिश (freq, वाई, 'आर'); साजिश (freq, एफ, 'बी'); Xlabel ('आवृत्ति, वा'); कथा ( 'संख्यात्मक', 'विश्लेषणात्मक'); रोक कर रखना;   
 आप अभी देख सकते हैं, मैं बस दो भूखंडों को समान दिखने की कोशिश कर रहा हूं। आखिरकार, मैं दो चीजों को करने का एक रास्ता खोजना चाहता हूं: 1) न्यूनतम नमूनाकरण दर, 2) नमूनों की न्यूनतम संख्या पाएं, एक त्रुटि (दो समाधानों के बीच के अंतर के एल 2 आदर्श के रूप में परिभाषित) तक पहुंचने के लिए 10 ^ -4।  
 मुझे लगता है कि यह बहुत आसान है, लेकिन मुझे ये भी लगता है कि ये दोनों ग्राफ़ भी दृष्टि से सहमत नहीं हैं। अगर कोई मुझे बताता है कि मैं गलत कहां रहा हूं और मैं ऊपर दो अंक (न्यूनतम नमूना आवृत्ति और न्यूनतम नमूनों की संख्या) से कैसे निपट सकता हूं, मैं बहुत सराहना करता हूं।  
 धन्यवाद  < नोट करने वाली पहली बात यह है कि फूरियर फ़ंक्शन के लिए जोड़ी बदलती है  exp (-t ^ 2)   
 +/- अनंत सीमा पर, जैसा कि से प्राप्त किया जा सकता है वास्तव में:  
 गाऊसी एफटी: एपीपी (-टी ^ 2) & lt; == & gt; एपीपी (- (2 ^ 2) / 4) * एसक्यूआरटी (पीआई)   
 अंत में, जैसा कि आप फ़ंक्शन बना रहे हैं  एक्सपी (-टी ^ 2) , आप  t  की सीमा को सकारात्मक मूल्यों तक सीमित कर रहे हैं (संपूर्ण +/- अनंत श्रेणी लेने के बजाय)। रिश्ते को पकड़ने के लिए, आपको इस तरह से  exp (-t ^ 2)  उत्पन्न करना होगा:  
  t = 0: dt: 10-dt ; % समय वेक्टर टी = टी - 0.5 * अधिकतम (टी); % के आसपास% t = 0 y = exp (- (ट। ^ 2));   
 फिर, वेरिएंट  w  रेडियन में कोणीय आवृत्ति प्रदर्शित करता है जो कि सामान्यीकृत आवृत्ति  freq  से संबंधित है:  < प्री>  w = 2 * pi * freq;   
 इस प्रकार,  
  F = (exp (- ((2 * pi * freq)। ^ 2) / 4)) * sqrt (pi) ; % विश्लेषणात्मक समाधान

Comments

Post a Comment

Popular posts from this blog

python - Receiving "KeyError" after decoding json result from url -

-
I am new to Python I am trying to parse JSON result from a URL. Basically, I was using the following: response = urllib.request.urlopen (url) json_obj = json.load (response) It should be a stroke "str 'not' bytes' in the lines of a given" JSON object, so after searching on the StackoverView Flo, I decode the response in this way: F = urllib.request.urlopen (Url) charset = f.info (). Get_param ('charset', 'utf8') data = f.read () decoded = json.loads (data.decode (charset)) If I print "decode" I is as follows: { 'link': { 'summary data': 'https: // localhost / piwebapi / streams / p0_7qHaW4UHU-RlCaz8tpasAAQAAAAU0hJTExNQU42NDIwXFNJTlVTT0lE / summary' 'value': 'https: // localhost / Piwebapi / streams / P0_7qHaW4UHU-RlCaz8tpasAAQAAAAU0hJTExNQU42NDIwXFNJTlVTT0lE / price ',' InterpolatedData ':' https: // localhost / Piwebapi / streams / P0_7qHaW4UHU-RlCaz8tpasAAQAAAAU0hJTE...
Read more

Out of index C/C++ -

-
सी / सी ++ में एक सरणी के लिए एक विशिष्ट घोषणा है: प्रकार का नाम [तत्व]; जहां एक प्रकार वैध है (जैसे कि इंट, फ्लोट ...), नाम एक मान्य पहचानकर्ता और तत्व फ़ील्ड है (जो हमेशा चौकोर ब्रैकेट्स [] में संलग्न होता है), सरणी की लंबाई को निर्दिष्ट करता है तत्वों की संख्या के संदर्भ में। इसलिए मैं एक सरणी की घोषणा करता हूँ जिसमें 2 तत्व हैं int a [2]; एक [3] = 4; क्यों यह अपवाद नहीं फेंकते? सीमा से बाहर की जाँच कुछ ऐसा है जो आप संभवतः जावा जैसे कुछ उच्च स्तर की भाषा से उपयोग किया जाता है। हालांकि सी / सी ++ में यह डिफ़ॉल्ट रूप से नहीं किया जाता है। यह आपको सरणी की सीमाओं की जांच करने के लिए एक छोटा प्रदर्शन हिट देता है और इसलिए सी तर्क यह है कि आप इसे स्वयं मैन्युअल रूप से कर सकते हैं, जब तक आप इसे सर्वश्रेष्ठ संभव प्रदर्शन की पेशकश कर सकते हैं। सी + + एसटीएल कंटेनर जैसे वेक्टर आमतौर पर एक पर () ऑपरेशन को बाध्य जांच करने के लिए समर्थन देता है और चूंकि आप [] अधिभार लोड कर सकते हैं - आप सरणी के लिए बाउंड-चेक सक्षम कर सकते हैं -स्टाइल एक्सेस। यदि सरणी एक कच्चे पॉइंटर ह...
Read more

java - ImportError: No module named py4j.java_gateway -

-
I'm trying to call the Python using the Java program py4j . I've been installing the plug-in Eclipse and test name Piidvi project. I'm trying to execute the following part of the code found on py4j webpage: Import from py4j.java_gateway to JavaGateway, java_import gateway = JavaGateway () jvm = gateway.jvm java_import (jvm, '' Org.eclipse Kkorkrisorsej. * ") Vrkspes_rut = Jvankresourkesplginkgetvrkspas (). GetRoot (Gateway .help (workspace_root, '* Projects *') project_names = [project.getName () (for projects workpace_root.getProjects))] print (Projekt_nam) But I There is an error in import. I have checked that the P4JJ is present in the Jar Eclipse plugin directory. Can anyone help please? I had to install the py4j application
Read more