Category: Matlab
Category Archives: Matlab
How to find GNSS satellite elevation angle from .obs and .nav data?
Trying to find GNSS/ NavIC satellite elevation angle from obs and nav dataTrying to find GNSS/ NavIC satellite elevation angle from obs and nav data Trying to find GNSS/ NavIC satellite elevation angle from obs and nav data gnss, elevation angle, calculation, .obs, .nav MATLAB Answers — New Questions
How can I use fitdist function
Hello,
I tried to use the function fitdist. The only result is an error:
Undefined function ‘fitdist’ for input arguments of type ‘double’.
Error in FitaNormalDistributiontoDataExample (line 13)
pd = fitdist(x,’Normal’)
This was the result of running the code from help file: C:Users…DocumentsMATLABExamplesFitaNormalDistributiontoDataExampleFitaNormalDistributiontoDataExample.m
%% Fit a Normal Distribution to Data
% Copyright 2015 The MathWorks, Inc.
%%
% Load the sample data. Create a vector containing the patients’ weight data.
load hospital
x = hospital.Weight;
%%
% Create a normal distribution object by fitting it to the data.
pd = fitdist(x,’Normal’)
%%
% Plot the pdf of the distribution.
x_values = 50:1:250;
y = pdf(pd,x_values);
plot(x_values,y,’LineWidth’,2)
What is wrong here?Hello,
I tried to use the function fitdist. The only result is an error:
Undefined function ‘fitdist’ for input arguments of type ‘double’.
Error in FitaNormalDistributiontoDataExample (line 13)
pd = fitdist(x,’Normal’)
This was the result of running the code from help file: C:Users…DocumentsMATLABExamplesFitaNormalDistributiontoDataExampleFitaNormalDistributiontoDataExample.m
%% Fit a Normal Distribution to Data
% Copyright 2015 The MathWorks, Inc.
%%
% Load the sample data. Create a vector containing the patients’ weight data.
load hospital
x = hospital.Weight;
%%
% Create a normal distribution object by fitting it to the data.
pd = fitdist(x,’Normal’)
%%
% Plot the pdf of the distribution.
x_values = 50:1:250;
y = pdf(pd,x_values);
plot(x_values,y,’LineWidth’,2)
What is wrong here? Hello,
I tried to use the function fitdist. The only result is an error:
Undefined function ‘fitdist’ for input arguments of type ‘double’.
Error in FitaNormalDistributiontoDataExample (line 13)
pd = fitdist(x,’Normal’)
This was the result of running the code from help file: C:Users…DocumentsMATLABExamplesFitaNormalDistributiontoDataExampleFitaNormalDistributiontoDataExample.m
%% Fit a Normal Distribution to Data
% Copyright 2015 The MathWorks, Inc.
%%
% Load the sample data. Create a vector containing the patients’ weight data.
load hospital
x = hospital.Weight;
%%
% Create a normal distribution object by fitting it to the data.
pd = fitdist(x,’Normal’)
%%
% Plot the pdf of the distribution.
x_values = 50:1:250;
y = pdf(pd,x_values);
plot(x_values,y,’LineWidth’,2)
What is wrong here? fitdist error MATLAB Answers — New Questions
Problem Saving UIFigures when running in parallel
Attempts to save UIFigure inside a parfor loop using getframe or exportgraphics dont seem to work. Returning the error that these functions cant be used in nodisplay mode. I rely on a 3rd party function ft_topoplot_er to plot certain topographs on my UIFigure and save the resulting video, but these calls take some time done sequentially, so I wished to parallelise them using parfor yet I’m unable to do so.Attempts to save UIFigure inside a parfor loop using getframe or exportgraphics dont seem to work. Returning the error that these functions cant be used in nodisplay mode. I rely on a 3rd party function ft_topoplot_er to plot certain topographs on my UIFigure and save the resulting video, but these calls take some time done sequentially, so I wished to parallelise them using parfor yet I’m unable to do so. Attempts to save UIFigure inside a parfor loop using getframe or exportgraphics dont seem to work. Returning the error that these functions cant be used in nodisplay mode. I rely on a 3rd party function ft_topoplot_er to plot certain topographs on my UIFigure and save the resulting video, but these calls take some time done sequentially, so I wished to parallelise them using parfor yet I’m unable to do so. parfor, uifigure, videowriter MATLAB Answers — New Questions
How many colors with their name present in the image.
I am looking for a way to find the number and name of certain colors present in the image so that I can use them to count for each pixel. I am doing it but actually some other colors or can say mixed colors are also present in the image which when I count are left in the counting and it counts less and as I increase the size of image pixels the uncounted number increase. I am giving the color ranges for the visible 13 colors in the image but their are some other colors present too which I am missing. How can I find them? Thanks in advane for the help, I am stuck with the problem since long time. My image isI am looking for a way to find the number and name of certain colors present in the image so that I can use them to count for each pixel. I am doing it but actually some other colors or can say mixed colors are also present in the image which when I count are left in the counting and it counts less and as I increase the size of image pixels the uncounted number increase. I am giving the color ranges for the visible 13 colors in the image but their are some other colors present too which I am missing. How can I find them? Thanks in advane for the help, I am stuck with the problem since long time. My image is I am looking for a way to find the number and name of certain colors present in the image so that I can use them to count for each pixel. I am doing it but actually some other colors or can say mixed colors are also present in the image which when I count are left in the counting and it counts less and as I increase the size of image pixels the uncounted number increase. I am giving the color ranges for the visible 13 colors in the image but their are some other colors present too which I am missing. How can I find them? Thanks in advane for the help, I am stuck with the problem since long time. My image is all colors number and name, color detection, all unique colors name and count in the image MATLAB Answers — New Questions
How do i connect matlab ui with simulink
i created matlab gui and the simulink model i want it to be combine how can i do those thinki created matlab gui and the simulink model i want it to be combine how can i do those think i created matlab gui and the simulink model i want it to be combine how can i do those think matlab ui, simulink MATLAB Answers — New Questions
Slow initialization at startup during the “Ready” step
I have installed Matlab R2024b on a Kubuntu 24.10 system. I installed only Matlab without an Add-ons and had no trouble during installation. Starting the program results in a long initialization time. I tracked the initialization process with the -timing flag. This results in the following output:
MATLAB Startup Performance Metrics (In Seconds)
total item gap description
=========================================================
0.47 0.47 0.00 MATLAB script
0.66 0.20 -0.00 main
0.71 0.05 0.00 Session Initialize
1.99 0.01 1.27 Toolbox cache load Start
1.99 0.44 0.84 Session Initialize
2.33 0.34 1.99 cachepath
2.50 0.11 2.38 LM Startup
2.60 0.10 0.00 splash
3.03 0.20 0.23 Constant Initialization
3.33 0.73 0.00 Engine Startup
3.52 0.20 0.00 InitSunVM
5.01 1.43 0.06 PostVMInit
5.01 1.68 0.00 mljInit
5.72 0.71 0.00 StartDesktop
5.72 2.39 0.00 Java initialization
5.72 0.00 0.00 psParser
5.72 0.00 0.00 Toolbox cache join
5.74 0.01 0.00 Package Registry initialization
6.17 0.21 0.22 matlabpath
7.08 0.19 0.72 matlabrc
7.09 0.00 0.00 Startup File
9.27 3.55 5.72 Init Desktop
217.84 212.12 5.72 Ready
=========================================================
Items shown account for 99.8% of total startup time
As visible, the last step consumes most of the time. Running the program bare with the -nojvm flag results in a quick start (3.92 s) . I already tried to increase the Java Memory heap space, as suggested in https://de.mathworks.com/matlabcentral/answers/1956139-slow-java-initialization-during-startup and tried several solutions from https://de.mathworks.com/matlabcentral/answers/92566-why-does-matlab-get-stuck-in-the-initializing-or-busy-state-or-take-a-long-time-to-start without improving the starting time.
Thanks in advance for any ideas how to resolve this issue.I have installed Matlab R2024b on a Kubuntu 24.10 system. I installed only Matlab without an Add-ons and had no trouble during installation. Starting the program results in a long initialization time. I tracked the initialization process with the -timing flag. This results in the following output:
MATLAB Startup Performance Metrics (In Seconds)
total item gap description
=========================================================
0.47 0.47 0.00 MATLAB script
0.66 0.20 -0.00 main
0.71 0.05 0.00 Session Initialize
1.99 0.01 1.27 Toolbox cache load Start
1.99 0.44 0.84 Session Initialize
2.33 0.34 1.99 cachepath
2.50 0.11 2.38 LM Startup
2.60 0.10 0.00 splash
3.03 0.20 0.23 Constant Initialization
3.33 0.73 0.00 Engine Startup
3.52 0.20 0.00 InitSunVM
5.01 1.43 0.06 PostVMInit
5.01 1.68 0.00 mljInit
5.72 0.71 0.00 StartDesktop
5.72 2.39 0.00 Java initialization
5.72 0.00 0.00 psParser
5.72 0.00 0.00 Toolbox cache join
5.74 0.01 0.00 Package Registry initialization
6.17 0.21 0.22 matlabpath
7.08 0.19 0.72 matlabrc
7.09 0.00 0.00 Startup File
9.27 3.55 5.72 Init Desktop
217.84 212.12 5.72 Ready
=========================================================
Items shown account for 99.8% of total startup time
As visible, the last step consumes most of the time. Running the program bare with the -nojvm flag results in a quick start (3.92 s) . I already tried to increase the Java Memory heap space, as suggested in https://de.mathworks.com/matlabcentral/answers/1956139-slow-java-initialization-during-startup and tried several solutions from https://de.mathworks.com/matlabcentral/answers/92566-why-does-matlab-get-stuck-in-the-initializing-or-busy-state-or-take-a-long-time-to-start without improving the starting time.
Thanks in advance for any ideas how to resolve this issue. I have installed Matlab R2024b on a Kubuntu 24.10 system. I installed only Matlab without an Add-ons and had no trouble during installation. Starting the program results in a long initialization time. I tracked the initialization process with the -timing flag. This results in the following output:
MATLAB Startup Performance Metrics (In Seconds)
total item gap description
=========================================================
0.47 0.47 0.00 MATLAB script
0.66 0.20 -0.00 main
0.71 0.05 0.00 Session Initialize
1.99 0.01 1.27 Toolbox cache load Start
1.99 0.44 0.84 Session Initialize
2.33 0.34 1.99 cachepath
2.50 0.11 2.38 LM Startup
2.60 0.10 0.00 splash
3.03 0.20 0.23 Constant Initialization
3.33 0.73 0.00 Engine Startup
3.52 0.20 0.00 InitSunVM
5.01 1.43 0.06 PostVMInit
5.01 1.68 0.00 mljInit
5.72 0.71 0.00 StartDesktop
5.72 2.39 0.00 Java initialization
5.72 0.00 0.00 psParser
5.72 0.00 0.00 Toolbox cache join
5.74 0.01 0.00 Package Registry initialization
6.17 0.21 0.22 matlabpath
7.08 0.19 0.72 matlabrc
7.09 0.00 0.00 Startup File
9.27 3.55 5.72 Init Desktop
217.84 212.12 5.72 Ready
=========================================================
Items shown account for 99.8% of total startup time
As visible, the last step consumes most of the time. Running the program bare with the -nojvm flag results in a quick start (3.92 s) . I already tried to increase the Java Memory heap space, as suggested in https://de.mathworks.com/matlabcentral/answers/1956139-slow-java-initialization-during-startup and tried several solutions from https://de.mathworks.com/matlabcentral/answers/92566-why-does-matlab-get-stuck-in-the-initializing-or-busy-state-or-take-a-long-time-to-start without improving the starting time.
Thanks in advance for any ideas how to resolve this issue. initializing, installation, linux, startup time MATLAB Answers — New Questions
Hardware Setup for C2000″ – Unable to Proceed After Selecting Process Family
Hey!
I am encountering an issue while performing the hardware setup for C2000 in MATLAB. Below are the details of the problem:
Setup Step:I am at the step where I need to select the process family for C2000 hardware setup.
Observed Issue: *After selecting the process family (e.g., C2000), the setup page is maximized, and I cannot proceed further. *The "Next" or "Proceed" button is either not visible or unresponsive.
Expected Behavior:After selecting the process family, I should be able to click "Next" or "Proceed" to move to the next step.
Steps Tried:
Resized the window manually.
Tried restarting MATLAB and redoing the hardware setup.
Checked for updates or patches related to MATLAB or C2000.
Can anyone guide me on how to resolve this issue? Any suggestions or troubleshooting steps would be greatly appreciated.
Thank you.Hey!
I am encountering an issue while performing the hardware setup for C2000 in MATLAB. Below are the details of the problem:
Setup Step:I am at the step where I need to select the process family for C2000 hardware setup.
Observed Issue: *After selecting the process family (e.g., C2000), the setup page is maximized, and I cannot proceed further. *The "Next" or "Proceed" button is either not visible or unresponsive.
Expected Behavior:After selecting the process family, I should be able to click "Next" or "Proceed" to move to the next step.
Steps Tried:
Resized the window manually.
Tried restarting MATLAB and redoing the hardware setup.
Checked for updates or patches related to MATLAB or C2000.
Can anyone guide me on how to resolve this issue? Any suggestions or troubleshooting steps would be greatly appreciated.
Thank you. Hey!
I am encountering an issue while performing the hardware setup for C2000 in MATLAB. Below are the details of the problem:
Setup Step:I am at the step where I need to select the process family for C2000 hardware setup.
Observed Issue: *After selecting the process family (e.g., C2000), the setup page is maximized, and I cannot proceed further. *The "Next" or "Proceed" button is either not visible or unresponsive.
Expected Behavior:After selecting the process family, I should be able to click "Next" or "Proceed" to move to the next step.
Steps Tried:
Resized the window manually.
Tried restarting MATLAB and redoing the hardware setup.
Checked for updates or patches related to MATLAB or C2000.
Can anyone guide me on how to resolve this issue? Any suggestions or troubleshooting steps would be greatly appreciated.
Thank you. c2000, hardware setup, f28p55x, microcontroller MATLAB Answers — New Questions
How to retain a program that has be trained by Radial Basis Function?
Hello every one
i made and trained a program using the nntool, i used the Radial Basis Function as a training algorithim. the results were good but i want them with better accuracy, is there any way to change the settings of the fucntion to give better answers, also can i get some details of the trained program such as: no. of iterations, error plot and etc.
thank youHello every one
i made and trained a program using the nntool, i used the Radial Basis Function as a training algorithim. the results were good but i want them with better accuracy, is there any way to change the settings of the fucntion to give better answers, also can i get some details of the trained program such as: no. of iterations, error plot and etc.
thank you Hello every one
i made and trained a program using the nntool, i used the Radial Basis Function as a training algorithim. the results were good but i want them with better accuracy, is there any way to change the settings of the fucntion to give better answers, also can i get some details of the trained program such as: no. of iterations, error plot and etc.
thank you radial basis function, ann MATLAB Answers — New Questions
Trying to use audiorecorder for impulse response testing in a for loop
Hi everyone,
I’m trying to test the impulse response of a microphone that i have attached to a rotary device connected to an arduino board. Basically trying to automate getting the impulse response at different angles. Thing is, I can’t get the audiorecorder function to work properly within the for loop. It only gives me the impulse response for the final iteration. The rest has no values.
Goal is to play a white noise, record the data from mic using audiorecorder, do my ffts and iffts for impulse response, rotate by a certain angle and repeat the process.
My code is below. Any help would be greatly appreciated because I’m at a complete standstill on what to do here.
close all
clear all
% Parameters
fs = 44100; % Sampling frequency
dur = 30; % Duration of recording in seconds
dt = 1/fs; % Time step
% Create time vector for plotting
time = 0:dt:dur-dt;
% Generate white noise for the entire duration
noise=randn(dur*fs,1);
ramp_up = 0:0.01:1;
ramp_up=ramp_up’;
ramp_down= 1:-0.01:0;
ramp_down=ramp_down’;
noise(1:length(ramp_up))=noise(1:length(ramp_up)).*ramp_up;
noise(end-length(ramp_down)+1:end)=noise(end-length(ramp_down)+1:end).*ramp_down;
noise=noise./max(noise);
recObj = audiorecorder(fs, 24, 2, 2);
%Arduino Setup
clear dcm sm shield a;
a = arduino(‘COM5′,’Uno’,’Libraries’,’AdafruitMotorShieldV2′);
shield = addon(a,’AdafruitMotorShieldV2′);
addrs = scanI2CBus(a,0);
sm = stepper(shield,2,200);
sm.RPM = 3;
Angle = 45; %Rotation Angle
n_rot = 180 / Angle; %Number of Rotations
IM_L=1000; % Truncated impulse response length in number of data points
LS_IR=zeros(n_rot,IM_L);
nfft = 96000; % FFT size
df = fs / nfft; % Frequency resolution
faxis = df:df:nfft*df; % Frequency axis
nFreq=length(faxis);
%Initializing data storage
nSamples = fs * dur;
allData_mic = zeros(nSamples, n_rot);
allData_spk = zeros(nSamples,n_rot);
allIR = zeros(nFreq, n_rot);
allIR_trunc=zeros(IM_L, n_rot);
allfrf=zeros(nFreq/2,n_rot);
allfrf_trunc=zeros(IM_L,n_rot);
%for i = 1:1
for i = 1:n_rot+1
%recording
sound(noise, fs);
% Record two-channel microphone input (24-bit, 2 channels)
recordblocking(recObj, dur); % Record for specified duration
% Get recorded data (two channels)
recordedData = getaudiodata(recObj);
spkinput = recordedData(:, 1); % First channel: Speaker
microphoneInput = recordedData(:, 2); % Second channel: Microphone
microphoneInput = microphoneInput./0.0462;%Left Ear (Fixed)
%microphoneInput = microphoneInput./0.0630;%Right Ear
%% FFT
nav = floor(length(spkinput) / nfft) * 2 – 1; % Number of averages
% Initialize accumulators for power spectral densities and cross-spectral density
G1 = zeros(nfft, 1); % Power Spectral Density of Speaker
G2 = zeros(nfft, 1); % Power Spectral Density of Microphone
G12 = zeros(nfft, 1); % Cross-Spectral Density between Microphone 1 and Microphone 2
G21 = zeros(nfft, 1);
% Loop through the signal with 50% overlap and apply FFT
for i = 1:nav
p1 = 1 + nfft * (i – 1) * 0.5; % Start index for the current block
p2 = p1 + nfft-1; % End index for the current block
% Apply Hann window to the blocks
win=hann(length(spkinput(p1:p2))); % Create a Hann window for the current segment
% Applying windows to blocks
spk_win = win.*spkinput(p1:p2); % Mic1 input
mic_win = win.*microphoneInput(p1:p2); % Mic2 input
% FFT of the windowed blocks
Xf = fft(spk_win, nfft); % FFT of Mic1
Yf = fft(mic_win, nfft); % FFT of Mic2
G1 = G1 + Xf .* conj(Xf); % PSD of Speaker input
G2 = G2 + Yf .* conj(Yf); % PSD of Microphone recording
G12 = G12 + Xf .* conj(Yf); % CSD between Microphone and Speaker
G21=G21 + Yf.*conj(Xf);
end
G1 = 2*G1(1:nfft/2)./(nfft*nfft*nav);
G2 = 2*G2(1:nfft/2)./(nfft*nfft*nav);
G12 = 2*G12(1:nfft/2)./(nfft*nfft*nav);
G21 = 2*G21(1:nfft/2)./(nfft*nfft*nav);
% Compute Coherence Function
coherence = abs(G21).^2 ./ (G1 .* G2);
% Compute Frequency Response Function (FRF)
frf = G2./G12; % FRF
magdB=20*log10(abs(frf));
%% Impulse Response
IR = ifft(frf,nfft);
IR_trunc=IR(1:IM_L);
%impulse_time = (0:nfft-1) * dt;
% n_samp=length(frf);
% impulse_time= 0:dt:(length(microphoneInput1)-1)*dt;
%store everything to the two vectors made above
allData_spk(:, i) = spkinput;
allData_mic(:, i) = microphoneInput;
allIR (:,i) = IR;
allIR_trunc (:,i)= IR_trunc;
allfrf (:,i) = frf;
filenameOut = sprintf(‘HeadTesting_Baseline_27Jan_%ddeg.mat’, ((i-1)*5));
%Movement———–
% – moves it to the left (counterclockwise) (standing behind it)
% + moves it to the right (clockwise) (stangind behind it)
if i~=n_rot+1
move(sm, +50*Angle);
end
end
clear dcm sm shield a;Hi everyone,
I’m trying to test the impulse response of a microphone that i have attached to a rotary device connected to an arduino board. Basically trying to automate getting the impulse response at different angles. Thing is, I can’t get the audiorecorder function to work properly within the for loop. It only gives me the impulse response for the final iteration. The rest has no values.
Goal is to play a white noise, record the data from mic using audiorecorder, do my ffts and iffts for impulse response, rotate by a certain angle and repeat the process.
My code is below. Any help would be greatly appreciated because I’m at a complete standstill on what to do here.
close all
clear all
% Parameters
fs = 44100; % Sampling frequency
dur = 30; % Duration of recording in seconds
dt = 1/fs; % Time step
% Create time vector for plotting
time = 0:dt:dur-dt;
% Generate white noise for the entire duration
noise=randn(dur*fs,1);
ramp_up = 0:0.01:1;
ramp_up=ramp_up’;
ramp_down= 1:-0.01:0;
ramp_down=ramp_down’;
noise(1:length(ramp_up))=noise(1:length(ramp_up)).*ramp_up;
noise(end-length(ramp_down)+1:end)=noise(end-length(ramp_down)+1:end).*ramp_down;
noise=noise./max(noise);
recObj = audiorecorder(fs, 24, 2, 2);
%Arduino Setup
clear dcm sm shield a;
a = arduino(‘COM5′,’Uno’,’Libraries’,’AdafruitMotorShieldV2′);
shield = addon(a,’AdafruitMotorShieldV2′);
addrs = scanI2CBus(a,0);
sm = stepper(shield,2,200);
sm.RPM = 3;
Angle = 45; %Rotation Angle
n_rot = 180 / Angle; %Number of Rotations
IM_L=1000; % Truncated impulse response length in number of data points
LS_IR=zeros(n_rot,IM_L);
nfft = 96000; % FFT size
df = fs / nfft; % Frequency resolution
faxis = df:df:nfft*df; % Frequency axis
nFreq=length(faxis);
%Initializing data storage
nSamples = fs * dur;
allData_mic = zeros(nSamples, n_rot);
allData_spk = zeros(nSamples,n_rot);
allIR = zeros(nFreq, n_rot);
allIR_trunc=zeros(IM_L, n_rot);
allfrf=zeros(nFreq/2,n_rot);
allfrf_trunc=zeros(IM_L,n_rot);
%for i = 1:1
for i = 1:n_rot+1
%recording
sound(noise, fs);
% Record two-channel microphone input (24-bit, 2 channels)
recordblocking(recObj, dur); % Record for specified duration
% Get recorded data (two channels)
recordedData = getaudiodata(recObj);
spkinput = recordedData(:, 1); % First channel: Speaker
microphoneInput = recordedData(:, 2); % Second channel: Microphone
microphoneInput = microphoneInput./0.0462;%Left Ear (Fixed)
%microphoneInput = microphoneInput./0.0630;%Right Ear
%% FFT
nav = floor(length(spkinput) / nfft) * 2 – 1; % Number of averages
% Initialize accumulators for power spectral densities and cross-spectral density
G1 = zeros(nfft, 1); % Power Spectral Density of Speaker
G2 = zeros(nfft, 1); % Power Spectral Density of Microphone
G12 = zeros(nfft, 1); % Cross-Spectral Density between Microphone 1 and Microphone 2
G21 = zeros(nfft, 1);
% Loop through the signal with 50% overlap and apply FFT
for i = 1:nav
p1 = 1 + nfft * (i – 1) * 0.5; % Start index for the current block
p2 = p1 + nfft-1; % End index for the current block
% Apply Hann window to the blocks
win=hann(length(spkinput(p1:p2))); % Create a Hann window for the current segment
% Applying windows to blocks
spk_win = win.*spkinput(p1:p2); % Mic1 input
mic_win = win.*microphoneInput(p1:p2); % Mic2 input
% FFT of the windowed blocks
Xf = fft(spk_win, nfft); % FFT of Mic1
Yf = fft(mic_win, nfft); % FFT of Mic2
G1 = G1 + Xf .* conj(Xf); % PSD of Speaker input
G2 = G2 + Yf .* conj(Yf); % PSD of Microphone recording
G12 = G12 + Xf .* conj(Yf); % CSD between Microphone and Speaker
G21=G21 + Yf.*conj(Xf);
end
G1 = 2*G1(1:nfft/2)./(nfft*nfft*nav);
G2 = 2*G2(1:nfft/2)./(nfft*nfft*nav);
G12 = 2*G12(1:nfft/2)./(nfft*nfft*nav);
G21 = 2*G21(1:nfft/2)./(nfft*nfft*nav);
% Compute Coherence Function
coherence = abs(G21).^2 ./ (G1 .* G2);
% Compute Frequency Response Function (FRF)
frf = G2./G12; % FRF
magdB=20*log10(abs(frf));
%% Impulse Response
IR = ifft(frf,nfft);
IR_trunc=IR(1:IM_L);
%impulse_time = (0:nfft-1) * dt;
% n_samp=length(frf);
% impulse_time= 0:dt:(length(microphoneInput1)-1)*dt;
%store everything to the two vectors made above
allData_spk(:, i) = spkinput;
allData_mic(:, i) = microphoneInput;
allIR (:,i) = IR;
allIR_trunc (:,i)= IR_trunc;
allfrf (:,i) = frf;
filenameOut = sprintf(‘HeadTesting_Baseline_27Jan_%ddeg.mat’, ((i-1)*5));
%Movement———–
% – moves it to the left (counterclockwise) (standing behind it)
% + moves it to the right (clockwise) (stangind behind it)
if i~=n_rot+1
move(sm, +50*Angle);
end
end
clear dcm sm shield a; Hi everyone,
I’m trying to test the impulse response of a microphone that i have attached to a rotary device connected to an arduino board. Basically trying to automate getting the impulse response at different angles. Thing is, I can’t get the audiorecorder function to work properly within the for loop. It only gives me the impulse response for the final iteration. The rest has no values.
Goal is to play a white noise, record the data from mic using audiorecorder, do my ffts and iffts for impulse response, rotate by a certain angle and repeat the process.
My code is below. Any help would be greatly appreciated because I’m at a complete standstill on what to do here.
close all
clear all
% Parameters
fs = 44100; % Sampling frequency
dur = 30; % Duration of recording in seconds
dt = 1/fs; % Time step
% Create time vector for plotting
time = 0:dt:dur-dt;
% Generate white noise for the entire duration
noise=randn(dur*fs,1);
ramp_up = 0:0.01:1;
ramp_up=ramp_up’;
ramp_down= 1:-0.01:0;
ramp_down=ramp_down’;
noise(1:length(ramp_up))=noise(1:length(ramp_up)).*ramp_up;
noise(end-length(ramp_down)+1:end)=noise(end-length(ramp_down)+1:end).*ramp_down;
noise=noise./max(noise);
recObj = audiorecorder(fs, 24, 2, 2);
%Arduino Setup
clear dcm sm shield a;
a = arduino(‘COM5′,’Uno’,’Libraries’,’AdafruitMotorShieldV2′);
shield = addon(a,’AdafruitMotorShieldV2′);
addrs = scanI2CBus(a,0);
sm = stepper(shield,2,200);
sm.RPM = 3;
Angle = 45; %Rotation Angle
n_rot = 180 / Angle; %Number of Rotations
IM_L=1000; % Truncated impulse response length in number of data points
LS_IR=zeros(n_rot,IM_L);
nfft = 96000; % FFT size
df = fs / nfft; % Frequency resolution
faxis = df:df:nfft*df; % Frequency axis
nFreq=length(faxis);
%Initializing data storage
nSamples = fs * dur;
allData_mic = zeros(nSamples, n_rot);
allData_spk = zeros(nSamples,n_rot);
allIR = zeros(nFreq, n_rot);
allIR_trunc=zeros(IM_L, n_rot);
allfrf=zeros(nFreq/2,n_rot);
allfrf_trunc=zeros(IM_L,n_rot);
%for i = 1:1
for i = 1:n_rot+1
%recording
sound(noise, fs);
% Record two-channel microphone input (24-bit, 2 channels)
recordblocking(recObj, dur); % Record for specified duration
% Get recorded data (two channels)
recordedData = getaudiodata(recObj);
spkinput = recordedData(:, 1); % First channel: Speaker
microphoneInput = recordedData(:, 2); % Second channel: Microphone
microphoneInput = microphoneInput./0.0462;%Left Ear (Fixed)
%microphoneInput = microphoneInput./0.0630;%Right Ear
%% FFT
nav = floor(length(spkinput) / nfft) * 2 – 1; % Number of averages
% Initialize accumulators for power spectral densities and cross-spectral density
G1 = zeros(nfft, 1); % Power Spectral Density of Speaker
G2 = zeros(nfft, 1); % Power Spectral Density of Microphone
G12 = zeros(nfft, 1); % Cross-Spectral Density between Microphone 1 and Microphone 2
G21 = zeros(nfft, 1);
% Loop through the signal with 50% overlap and apply FFT
for i = 1:nav
p1 = 1 + nfft * (i – 1) * 0.5; % Start index for the current block
p2 = p1 + nfft-1; % End index for the current block
% Apply Hann window to the blocks
win=hann(length(spkinput(p1:p2))); % Create a Hann window for the current segment
% Applying windows to blocks
spk_win = win.*spkinput(p1:p2); % Mic1 input
mic_win = win.*microphoneInput(p1:p2); % Mic2 input
% FFT of the windowed blocks
Xf = fft(spk_win, nfft); % FFT of Mic1
Yf = fft(mic_win, nfft); % FFT of Mic2
G1 = G1 + Xf .* conj(Xf); % PSD of Speaker input
G2 = G2 + Yf .* conj(Yf); % PSD of Microphone recording
G12 = G12 + Xf .* conj(Yf); % CSD between Microphone and Speaker
G21=G21 + Yf.*conj(Xf);
end
G1 = 2*G1(1:nfft/2)./(nfft*nfft*nav);
G2 = 2*G2(1:nfft/2)./(nfft*nfft*nav);
G12 = 2*G12(1:nfft/2)./(nfft*nfft*nav);
G21 = 2*G21(1:nfft/2)./(nfft*nfft*nav);
% Compute Coherence Function
coherence = abs(G21).^2 ./ (G1 .* G2);
% Compute Frequency Response Function (FRF)
frf = G2./G12; % FRF
magdB=20*log10(abs(frf));
%% Impulse Response
IR = ifft(frf,nfft);
IR_trunc=IR(1:IM_L);
%impulse_time = (0:nfft-1) * dt;
% n_samp=length(frf);
% impulse_time= 0:dt:(length(microphoneInput1)-1)*dt;
%store everything to the two vectors made above
allData_spk(:, i) = spkinput;
allData_mic(:, i) = microphoneInput;
allIR (:,i) = IR;
allIR_trunc (:,i)= IR_trunc;
allfrf (:,i) = frf;
filenameOut = sprintf(‘HeadTesting_Baseline_27Jan_%ddeg.mat’, ((i-1)*5));
%Movement———–
% – moves it to the left (counterclockwise) (standing behind it)
% + moves it to the right (clockwise) (stangind behind it)
if i~=n_rot+1
move(sm, +50*Angle);
end
end
clear dcm sm shield a; audiorecorder, signal processing, audio, arduino, for loop MATLAB Answers — New Questions
MATLAB GRADER VIEW LEARNER ANSWERS BY LEARNER ID
Dear All,
I am currently marking students’ assignments in MATLAB Grader. The assignments involve several questions for learners to answer. One good feature of MATLAB Grader is the ability to search for answers using the solution ID.
However, it would be helpful to have a search option using the learner ID, which would bring up all their answers for the respective assignment. I’m not sure if this feature already exists.
Any help would be appreciated.
Regards,
FaisalDear All,
I am currently marking students’ assignments in MATLAB Grader. The assignments involve several questions for learners to answer. One good feature of MATLAB Grader is the ability to search for answers using the solution ID.
However, it would be helpful to have a search option using the learner ID, which would bring up all their answers for the respective assignment. I’m not sure if this feature already exists.
Any help would be appreciated.
Regards,
Faisal Dear All,
I am currently marking students’ assignments in MATLAB Grader. The assignments involve several questions for learners to answer. One good feature of MATLAB Grader is the ability to search for answers using the solution ID.
However, it would be helpful to have a search option using the learner ID, which would bring up all their answers for the respective assignment. I’m not sure if this feature already exists.
Any help would be appreciated.
Regards,
Faisal matlab grader, search, userid MATLAB Answers — New Questions
Removing comas from the text file
I have a text file containing some numbers separated with coma.How can I remove those comas.
Example:
data.txt={1,2,3,4………}.Is there any method to remove those comas?
My text file when loaded to matlab is 2517×1001 matrix.
Kindly respondI have a text file containing some numbers separated with coma.How can I remove those comas.
Example:
data.txt={1,2,3,4………}.Is there any method to remove those comas?
My text file when loaded to matlab is 2517×1001 matrix.
Kindly respond I have a text file containing some numbers separated with coma.How can I remove those comas.
Example:
data.txt={1,2,3,4………}.Is there any method to remove those comas?
My text file when loaded to matlab is 2517×1001 matrix.
Kindly respond removing comas MATLAB Answers — New Questions
programing in Simulink a lamp
i need the lamp to display different color when it gets certain value from the matlab function. i add another Matlab Function to program the lamp but i am not sure if that such a thing to exist.
this is simulink. Here is my code that is in the Matlab Function to program the lamp as you see how i connect the lamp to the function. the code below is the code inside the matlab function. i am not even sure if my code is right at all.
function lampStates = caution(distance)
lampState1.Value = ( distance >= 0.0) && ( distance <= 0.5);
lampState2.Value = ( distance == 0.5);
lampState3.Value = ( distance == 1.5);
lampState4.Value = ( distance == 2.5);
lampState5.Value = ( distance >= 2.5);
if ( distance >= 0.0) && ( distance <= 0.5)
lampState1.Color = [254 51 10];
elseif ( distance == 0.5)
lampState2.Color = [237 177 32];
elseif ( distance == 1.5)
lampState3.Color = [100 212 19];
elseif ( distance == 2.5)
lampState4.Color = [237 177 32];
elseif ( distance >= 2.5)
lampState5.Color = [254 51 10];
end
lampStates = [lampState1 lampState2 lampState3 lampState4 lampState5];i need the lamp to display different color when it gets certain value from the matlab function. i add another Matlab Function to program the lamp but i am not sure if that such a thing to exist.
this is simulink. Here is my code that is in the Matlab Function to program the lamp as you see how i connect the lamp to the function. the code below is the code inside the matlab function. i am not even sure if my code is right at all.
function lampStates = caution(distance)
lampState1.Value = ( distance >= 0.0) && ( distance <= 0.5);
lampState2.Value = ( distance == 0.5);
lampState3.Value = ( distance == 1.5);
lampState4.Value = ( distance == 2.5);
lampState5.Value = ( distance >= 2.5);
if ( distance >= 0.0) && ( distance <= 0.5)
lampState1.Color = [254 51 10];
elseif ( distance == 0.5)
lampState2.Color = [237 177 32];
elseif ( distance == 1.5)
lampState3.Color = [100 212 19];
elseif ( distance == 2.5)
lampState4.Color = [237 177 32];
elseif ( distance >= 2.5)
lampState5.Color = [254 51 10];
end
lampStates = [lampState1 lampState2 lampState3 lampState4 lampState5]; i need the lamp to display different color when it gets certain value from the matlab function. i add another Matlab Function to program the lamp but i am not sure if that such a thing to exist.
this is simulink. Here is my code that is in the Matlab Function to program the lamp as you see how i connect the lamp to the function. the code below is the code inside the matlab function. i am not even sure if my code is right at all.
function lampStates = caution(distance)
lampState1.Value = ( distance >= 0.0) && ( distance <= 0.5);
lampState2.Value = ( distance == 0.5);
lampState3.Value = ( distance == 1.5);
lampState4.Value = ( distance == 2.5);
lampState5.Value = ( distance >= 2.5);
if ( distance >= 0.0) && ( distance <= 0.5)
lampState1.Color = [254 51 10];
elseif ( distance == 0.5)
lampState2.Color = [237 177 32];
elseif ( distance == 1.5)
lampState3.Color = [100 212 19];
elseif ( distance == 2.5)
lampState4.Color = [237 177 32];
elseif ( distance >= 2.5)
lampState5.Color = [254 51 10];
end
lampStates = [lampState1 lampState2 lampState3 lampState4 lampState5]; display different color lamp MATLAB Answers — New Questions
LAG VALUE OF A VARIABLE
Hello,
I want to have a lag value of variable X (PANEL of nature), but since Matlab does not accept different matrix size variables, what should I do. Variable is 4100*1 and lag of it is 3936*1.Hello,
I want to have a lag value of variable X (PANEL of nature), but since Matlab does not accept different matrix size variables, what should I do. Variable is 4100*1 and lag of it is 3936*1. Hello,
I want to have a lag value of variable X (PANEL of nature), but since Matlab does not accept different matrix size variables, what should I do. Variable is 4100*1 and lag of it is 3936*1. lag values MATLAB Answers — New Questions
syms error, cannot use syms command
I am having an issue using the syms command:
syms x(t) y(t)
getting the following error:
Error using assignin
Invalid variable name "x(t)" in ASSIGNIN.
Error in syms (line 56)
assignin(‘caller’,varargin{i},sym(varargin{i}));
Error in Untitled (line 1)
syms x(t) y(t)
I am tried the following solution but I can not solve my problem"
which -all syms
C:Program FilesMATLABR2018btoolboxmaplesyms.m
My Matlab ver: R2018b
On my system Maple 2018 was installed at: C:Program FilesMaple 2018I am having an issue using the syms command:
syms x(t) y(t)
getting the following error:
Error using assignin
Invalid variable name "x(t)" in ASSIGNIN.
Error in syms (line 56)
assignin(‘caller’,varargin{i},sym(varargin{i}));
Error in Untitled (line 1)
syms x(t) y(t)
I am tried the following solution but I can not solve my problem"
which -all syms
C:Program FilesMATLABR2018btoolboxmaplesyms.m
My Matlab ver: R2018b
On my system Maple 2018 was installed at: C:Program FilesMaple 2018 I am having an issue using the syms command:
syms x(t) y(t)
getting the following error:
Error using assignin
Invalid variable name "x(t)" in ASSIGNIN.
Error in syms (line 56)
assignin(‘caller’,varargin{i},sym(varargin{i}));
Error in Untitled (line 1)
syms x(t) y(t)
I am tried the following solution but I can not solve my problem"
which -all syms
C:Program FilesMATLABR2018btoolboxmaplesyms.m
My Matlab ver: R2018b
On my system Maple 2018 was installed at: C:Program FilesMaple 2018 syms error MATLAB Answers — New Questions
Invalid operand. Variables of type “sym” cannot be combined with other models.
Hello, how can I resolve the error "Invalid operand. Variables of type "sym" cannot be combined with other models." in this code. Separately if it works for me. Thanks
d=5;%sobreoscilacion
tp=4; %tiempo pico
k=100;
syms e w
ecA=((-e*pi)/sqrt(1-e^2))-log(d)== 0;
e= vpa(solve(ecA,e)); %valor epsilon
ecB=(pi/(w*(sqrt(1-e^2))))-tp==0;
w= vpa(solve(ecB,w)); %valor ohmega
%%%%
s=tf(‘s’);
Q=((k*(w)^2)/(s^2+2*e+s+(w)^2))*exp(-tp*s);
Q_apro=pade(Q,1);Hello, how can I resolve the error "Invalid operand. Variables of type "sym" cannot be combined with other models." in this code. Separately if it works for me. Thanks
d=5;%sobreoscilacion
tp=4; %tiempo pico
k=100;
syms e w
ecA=((-e*pi)/sqrt(1-e^2))-log(d)== 0;
e= vpa(solve(ecA,e)); %valor epsilon
ecB=(pi/(w*(sqrt(1-e^2))))-tp==0;
w= vpa(solve(ecB,w)); %valor ohmega
%%%%
s=tf(‘s’);
Q=((k*(w)^2)/(s^2+2*e+s+(w)^2))*exp(-tp*s);
Q_apro=pade(Q,1); Hello, how can I resolve the error "Invalid operand. Variables of type "sym" cannot be combined with other models." in this code. Separately if it works for me. Thanks
d=5;%sobreoscilacion
tp=4; %tiempo pico
k=100;
syms e w
ecA=((-e*pi)/sqrt(1-e^2))-log(d)== 0;
e= vpa(solve(ecA,e)); %valor epsilon
ecB=(pi/(w*(sqrt(1-e^2))))-tp==0;
w= vpa(solve(ecB,w)); %valor ohmega
%%%%
s=tf(‘s’);
Q=((k*(w)^2)/(s^2+2*e+s+(w)^2))*exp(-tp*s);
Q_apro=pade(Q,1); syms, sym MATLAB Answers — New Questions
figure does not print or print preview
When I try print a figure or even print preview the figure I get the error below. This happens with 2019b version. Thanks.
rror using printopt
Too many input arguments.
Error in matlab.graphics.internal.mlprintjob
Error in printjob (line 22)
pj = matlab.graphics.internal.mlprintjob;
Error in print>LocalCreatePrintJob (line 102)
pj = printjob([handles{:}]);
Error in print (line 38)
[pj, inputargs] = LocalCreatePrintJob(varargin{:});
Error in printdlg (line 63)
print(Fig, ‘-v’);
Error while evaluating Menu Callback.When I try print a figure or even print preview the figure I get the error below. This happens with 2019b version. Thanks.
rror using printopt
Too many input arguments.
Error in matlab.graphics.internal.mlprintjob
Error in printjob (line 22)
pj = matlab.graphics.internal.mlprintjob;
Error in print>LocalCreatePrintJob (line 102)
pj = printjob([handles{:}]);
Error in print (line 38)
[pj, inputargs] = LocalCreatePrintJob(varargin{:});
Error in printdlg (line 63)
print(Fig, ‘-v’);
Error while evaluating Menu Callback. When I try print a figure or even print preview the figure I get the error below. This happens with 2019b version. Thanks.
rror using printopt
Too many input arguments.
Error in matlab.graphics.internal.mlprintjob
Error in printjob (line 22)
pj = matlab.graphics.internal.mlprintjob;
Error in print>LocalCreatePrintJob (line 102)
pj = printjob([handles{:}]);
Error in print (line 38)
[pj, inputargs] = LocalCreatePrintJob(varargin{:});
Error in printdlg (line 63)
print(Fig, ‘-v’);
Error while evaluating Menu Callback. print MATLAB Answers — New Questions
select only rows and not columns
Does anybody knows how to select only rows and not columns? I’m trying to select all the first values of each column, is that possible?Does anybody knows how to select only rows and not columns? I’m trying to select all the first values of each column, is that possible? Does anybody knows how to select only rows and not columns? I’m trying to select all the first values of each column, is that possible? rows, columns, matrix MATLAB Answers — New Questions
SIL on autosar composition, Different Types error
I have an autsoar architecture model(autosar composition).
At first I created a test harness on the architecture model and wanted to make sil test. I got error it is not supported I shall set the SimulationMode of architecture model to normal and the SimulationMode of the autosar compome below to the SIL. I did it
Then I got the error Different types as below shown
I follow another way as mathwork proposal, Testing AUTOSAR Compositions
At the end I got the same error Different types with no any helpful info what going wrong :(I have an autsoar architecture model(autosar composition).
At first I created a test harness on the architecture model and wanted to make sil test. I got error it is not supported I shall set the SimulationMode of architecture model to normal and the SimulationMode of the autosar compome below to the SIL. I did it
Then I got the error Different types as below shown
I follow another way as mathwork proposal, Testing AUTOSAR Compositions
At the end I got the same error Different types with no any helpful info what going wrong 🙁 I have an autsoar architecture model(autosar composition).
At first I created a test harness on the architecture model and wanted to make sil test. I got error it is not supported I shall set the SimulationMode of architecture model to normal and the SimulationMode of the autosar compome below to the SIL. I did it
Then I got the error Different types as below shown
I follow another way as mathwork proposal, Testing AUTOSAR Compositions
At the end I got the same error Different types with no any helpful info what going wrong 🙁 autosar sil test MATLAB Answers — New Questions
Unterminated quote in @-file: n” error when appending pdfs
error when appending pdfserror when appending pdfs error when appending pdfs append pdf error MATLAB Answers — New Questions
How to add extra text to plot legend?
Hi, guys
Let us suppose that I have the following simple code:
clc; clear all; close all;
x = linspace(0, 2*pi, 100);
figure
hold on;grid on; box on;
plot(x,sin(x));
plot(x, cos(x));
plot(x, tan(x));
axis([0 2*pi -4 4]);
legend(‘sin’, ‘cos’, ‘tan’);
My goal is to add some extra text in the legend box as shown in the following image example (image modified with gimp):
Can you help me to codify this task?Hi, guys
Let us suppose that I have the following simple code:
clc; clear all; close all;
x = linspace(0, 2*pi, 100);
figure
hold on;grid on; box on;
plot(x,sin(x));
plot(x, cos(x));
plot(x, tan(x));
axis([0 2*pi -4 4]);
legend(‘sin’, ‘cos’, ‘tan’);
My goal is to add some extra text in the legend box as shown in the following image example (image modified with gimp):
Can you help me to codify this task? Hi, guys
Let us suppose that I have the following simple code:
clc; clear all; close all;
x = linspace(0, 2*pi, 100);
figure
hold on;grid on; box on;
plot(x,sin(x));
plot(x, cos(x));
plot(x, tan(x));
axis([0 2*pi -4 4]);
legend(‘sin’, ‘cos’, ‘tan’);
My goal is to add some extra text in the legend box as shown in the following image example (image modified with gimp):
Can you help me to codify this task? legend, text, plots MATLAB Answers — New Questions
