আমি কীভাবে তরঙ্গসী রূপান্তরটির বনাম প্রস্থের প্লট করতে পারি?

Question

আমি মরলেট অবিচ্ছিন্ন ওয়েভলেট রূপান্তর চালাচ্ছি। আমি wscalogramসিগন্যাল পেয়েছি এবং এখন আমি নীচের ছবির মতো ফ্রিক-ম্যাগনিটিউটি প্লট করতে চাই, তবে এটি কীভাবে করব তা আমি জানি না:

scal2freqস্কেলগুলি সিউডো ফ্রিকোয়েন্সিতে রূপান্তর করতে আমি ম্যাটল্যাব ফাংশন ব্যবহার করেছি । এছাড়াও আমার সিগন্যালে আমার কিছু ফ্রিকোয়েন্সি রয়েছে যার একটি বড় স্যাঁতসেঁতে অনুপাত (4%) থাকে তাই তারা প্লটে ভালভাবে দৃশ্যমান হয় না। আমি কীভাবে এই অত্যধিক স্যাঁতসেঁতে মোডকে বাড়াতে পারি?

আমি ম্যাটল্যাব ব্যবহার করছি, এখানে আমার কোডগুলি রয়েছে:

% Import the text4.txt to matlab workspace. and save it under name "data"
t=linspace(0,30,301);
Fs=ceil(inv(t(2)-t(1)));
x=data(:,4); % use x=data(:,3),x=data(:,5) too. first column is time,second is refrence
wname = 'morl';
scales = 1:1:256;
coefs = cwt(x,scales,wname,'lvlabs');
freq = scal2frq(scales,wname,1/Fs);
surf(t,freq,abs(coefs));shading('interp');
axis tight; xlabel('Seconds'); ylabel('Pseudo-Frequency (Hz)');
axis([0 30 0 1 0 60])
xlabel('Time'); ylabel('Frequency');
figure;
sc=wscalogram('image',coefs,'scales',freq,'ydata',x);
hold on 

আমার বর্তমান প্লটটি এরকম দেখাচ্ছে:

Comments

আপনার তরঙ্গলিপি ট্রান্সফর্ম ডেটা কোন ফর্ম গ্রহণ করে?

— জিম ক্লে

হাই @ জিমক্লে! আমি আপনার প্রশ্নটি বুঝতে পেরেছি কিনা তা নিশ্চিত নই। তবে যদি আপনি এটির মতো দেখতে বোঝাতে চান তবে এটি দ্বিতীয় প্রশ্নটির মতো যা আমি প্রশ্নটিতে বজায় রেখেছি। এবং আমি 1 ম লিঙ্কের মতো প্লট চাই। (২ য় [মাইপ্লট]) cubeupload.com/im/bSSlMI.jpg 1 ম [পছন্দসই] cubeupload.com/im/KbhGMI.jpg এখানে আমার সিগন্যাল < filewap.com/dl/Lwf7n96fAt/test4.txt.html > এখানে আমার এমফাইল < ফাইলস্যাভডড / ডিএল / জিএমএসএলবিএফডিবি / এমফিল.টিএসটিএইচটিএমএল > শুভেচ্ছা।

— বৈদ্যুতিনবিদ

আপনার প্রশ্নটি আসলে কী তা আমি নিশ্চিত নই ...

— তারিন জিয়াআই

হাই @ ব্যবহারকারী 4619; প্রথম চক্রান্ত দেখুন। আপনি ওয়েভলেট স্ক্লোগ্রামের [যা মাতলাব দ্বারা প্রাপ্ত] এর ডান দিকে একটি প্লট [ফ্রিক-ম্যাগনিটিউড] দেখতে পাচ্ছেন। দ্বিতীয় ছবিটি আমার ওয়েভলেট স্কেলোগ্রাম। তবে আমি কীভাবে প্রথম প্লটের ডান পাশের মতো [ফ্রিক-ম্যাগনিটিউড] প্লট করতে জানি না। কোন সাহায্য plz?

— বৈদ্যুতিনবিদ

@ ইলেক্ট্রিকম্যান একটি স্পেকট্রামগ্রাম স্ক্লোগ্রামের মতো নয়। আপনি যদি কোনও স্ক্লোগ্রামে যাচ্ছেন তবে কোনও ফ্রিকোয়েন্সি ভিএস সময় নেই, কেবলমাত্র ভিএস সময় স্কেল করুন। তাহলে আপনি ঠিক কোনটি করছেন?

— তারিন জিয়াই

Answer 1

আপনার কোডটি কিছুটা সংশোধন করা হচ্ছে, তবে কোনও বড় পরিবর্তন হয়নি এবং আমি যেভাবেই সঠিক ফলাফল পেয়েছি। এই টেমপ্লেট কোডটি এখানে ব্যবহার করুন এবং আপনার কোনও সমস্যা দেখা উচিত নয়। আমি সঠিক ফলাফল পেতে।

clear all; 
t=linspace(0,30,301); 
Fs = (inv(t(2)-t(1))); 
x=randn(100,1);  
wname = 'morl'; 
scales = 1:1:256; 
chefs = cwt(x,scales,wname,'lvlabs');

freq = scal2frq(scales,wname,1/Fs);

figure; 
coefsSquared = abs(coefs).^2; 
imagesc(coefsSquared); 
grid off;

%Pick one of the columsn to plot: 
figure; 
plot(coefsSquared(:,47))

আপনি যে চিত্রটি রেখেছিলেন তার ডানদিকে চিত্রটি পেতে, এটি কেবল সময়-স্লাইস বর্ণালী। অন্য কথায়, আপনার ফ্রিকোয়েন্সি-সময় ম্যাট্রিক্সের একটি কলাম।

উদাহরণস্বরূপ, আমরা ফ্রিকোয়েন্সি বনাম সময়ের মধ্যে নিম্নলিখিত স্কোয়ার্ড-সহগ ম্যাট্রিক্স পেতে পারি:

আসুন এখন ধরা যাক আমরা কলাম 47-তে টাইম-স্লাইসের বর্ণালীটি দেখতে চাই। তারপরে এটি এখানে দেওয়া হল:

এখানে, তারপরে আপনি কলাম 47 এর সাথে সম্পর্কিত টাইম-স্লাইসের জন্য সহগের, বনাম, ফ্রিকোয়েন্সিটির শক্তি দেখতে পাবেন।

Comments

কোড পরিবর্তন করার জন্য ধন্যবাদ। তবে আমি এখনও আমার প্রশ্নের উত্তর পাইনি question প্রশ্নটি হ'ল আমি কীভাবে ফ্রিক-দৈর্ঘ্য বা স্কেল-परिमाण [কোনও পার্থক্য] প্লট করতে পারি না। আমি প্রথম ছবির ডান পাশের মতো প্লট চাই [[এটি প্রশ্নে আছে]। অগ্রিম ধন্যবাদ

— ইলেকট্রিকম্যান 19

@ ইলেক্ট্রিকম্যান ডান হাতের প্লটটি আপনার ফ্রিকোয়েন্সি ভিএস টাইম ম্যাট্রিক্সের কেবল একটি কলাম।

— তারিন জিয়াই

- @ ব্যবহারকারী 4619 এটি ফ্রিক বনাম সময় নয়, সময়সীমা 0-20 এবং ফ্রিকোয়েন্সি রেঞ্জ 0-2.5 তবে ডান হাতের প্লটটি ফ্রিক বনাম প্রস্থের। আমি নিশ্চিত নই যে এর ঠিক কয়ফস বা মাত্রা তবে স্পষ্টতই এটি কোফ বা কোফে to 2 এর সমানুপাতিক। tnx

— বৈদ্যুতিনবিদ

@ ইলেক্ট্রিকম্যান দয়া করে সম্পাদনাগুলি দেখুন। এটা কি এখন পরিষ্কার?

— তারিন জিয়াই

- @ ব্যবহারকারী 4619 এটি কাজ করে। আপনার উত্তরের জন্য ধন্যবাদ. আমি মাইফিলের শেষে একটি সামান্য কোড রেখেছি। এবং উল্লিখিত প্লটের ডান হাতের মতো একই প্লট পেয়েছি। cof2 = ABS (coefs) ^ 2। cof2trans = cof2 '; maxmods = MAX (cof2trans); চিত্রে; প্লট (ম্যাকমোডস, ফ্রিক) অক্ষ ([0 60 0 1]) শুভেচ্ছা

— বৈদ্যুতিনিক

Answer 2

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

F = scal2frq (A, 'wname', DELTA) এ দ্বারা প্রদত্ত স্কেলগুলির সাথে সম্পর্কিত সিউডো-ফ্রিকোয়েন্সিগুলি দেয়, তরঙ্গীকরণ ফাংশন 'wname' (আরও তথ্যের জন্য ওয়েভফান দেখুন) এবং স্যাম্পলিং পিরিয়া DELTA দেয়।

আপনি যেমন পড়তে পারেন, ফ্রিকোয়েন্সিগুলি সিগন্যালের আসল ফ্রিকোয়েন্সি নয়, তবে আনুমানিক। আপনি যদি আসল ফ্রিকোয়েন্সি চান তবে আপনার স্বল্প সময়ের ফুরিয়ার ট্রান্সফর্ম (এসটিএফটি) ব্যবহার করা উচিত।

সুতরাং, যখন কোন সিডাব্লুটি ব্যবহার করা উচিত? একটি সাধারণ নীতি হিসাবে, সংকেতগুলি ক্ষণস্থায়ী হলে (যেমন দ্রুত পরিবর্তন হয়) তরঙ্গপত্রগুলি ভালভাবে কাজ করে ।

সম্পাদনা: হাহা, আমি জানি না আমি এখানে কী করছি। দুঃখিত :-D আপনার প্রশ্নের বিন্দুটি পুরোপুরি মিস করেছে। শুধু শিরোনাম পড়ুন ....

Comments

@ ওয়েভ আপনার তথ্য সঠিক তবে কি? আমি কেবল প্রথম ছবির ডান পাশের মতো একটি প্লট রাখতে চাই। freq বনাম परिमाण। এবং যতক্ষণ না আমি স্কেল 2ফ্রেইক ব্যবহার করতে পারি ততক্ষণ যদি আমি স্কেল বনাম প্রস্থের মাত্রা পাই তবে তা ঠিক।

— ইলেকট্রিকম্যান

"আপনি যদি আসল ফ্রিকোয়েন্সি চান তবে আপনার স্বল্প সময়ের ফুরিয়ার ট্রান্সফর্ম (এসটিএফটি) ব্যবহার করা উচিত।" মরলেট তরঙ্গলেটের যেমন গাউস-উইন্ডো এসটিএফটির ফ্রিকোয়েন্সি রয়েছে তেমন ফ্রিকোয়েন্সি রয়েছে। উভয় রূপান্তরগুলির "পরমাণু" একই রকম।

— endolith

Answer 3

আপনার এখনও সহায়তা প্রয়োজন কিনা তা নিশ্চিত নন ...

আপনি যদি ফ্রিকোয়েন্সি বনাম প্রস্থ চান তবে সহগের প্লট করতে একটি "জাল / সার্ফ" কমান্ড ব্যবহার করুন।

তারপরে প্লটের ওরিয়েন্টেশনটি পরিবর্তন করুন যাতে আপনি যে ফলাফলটি সন্ধান করেন তা পান। এর জন্য "ভিউ" কমান্ডটি দেখুন।

Comments

ধন্যবাদ, আমি বেশ কয়েক মাস আগে এই সমস্যাটি সমাধান করেছি এবং আমার ফলাফল প্রকাশ করেছি। তবে আমি আপনার বিকল্পটি পরের বার চেষ্টা করব b তবে দয়া করে আপনার উত্তরটি মন্তব্য করার জায়গায় নিয়ে যান। শুভকামনা @ asilva732

— বৈদ্যুতিনবিদ