I have a steady state heat transfer problem trying to model a thin film heater wedged between two significantly thicker slabs. The slab thicknesses may be as much as three orders of magnitude larger than the heater. I want temperature resolution within the slabs and a temperature of the heater but I don’t need temperature resolution within the film per se. I tried treating the entire thin film as dissipating energy using an internal heat source command of this form [internalHeatSource(thermalmodel,6000,"Face",2)]. My checks on model validity included looking for higher temperature within the heater compared to the two bounding slabs using a command of this form [pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), … "Contour","on",… "ColorMap","jet")]. I also checked energy conservation in the steady state by comparing the known amount of power dissipated within the heater to the energy rate leaving the boundaries of the system using commands of the form [Qn1= evaluateHeatRate(thermalresults, "Edge",1)] for each edge and summing up these energy rates..
This approach was not successful. I did not see significant temperature increases that I would expect for heater power levels that I used. Conservation of energy was not observed. I plotted the mesh results and in some or all cases would get strange looking results. In some or all cases Matlab returned the following feedback: “Warning: Found elements with poor shape quality.” I tried manipulating the mesh size to get finer resolution using commands of the form [mesh= generateMesh(thermalmodel, ‘Hmax’, 1e-6, ‘Hmin’, 1e-7)]. I was not successful.
I feel like the problem is related to the mismatch in mesh resolution that may be required in the heater compared to the larger slabs, but I don’t know if that’s true and I don’t know how to handle it if that is true.
The code is included below.
Any help would be greatly appreciated.
john

%% Create a thermal model for transient analysis
thermalmodel = createpde("thermal","steadystate");

% Define the geometry components
Rfilmbottom = [3,4,[-5.0, 5.0, 5.0, -5.0, 0.0, 0.0, 0.124375,0.124375]/1000]’;
Rfilmhtr = [3,4,[-5.0, 5.0, 5.0, -5.0, 0.124375,0.124375,0.125625,0.125625]/1000]’;
Rfilmtop= [3,4,[-5.0, 5.0, 5.0, -5.0, 0.125625,0.125625, 0.25, 0.25]/1000]’;

% Combine the geometry data
gdm = [ Rfilmbottom Rfilmhtr Rfilmtop ];
ns = char(‘Rfilmbottom’, ‘Rfilmhtr’, ‘Rfilmtop’);
g = decsg(gdm, ‘Rfilmbottom + Rfilmhtr + Rfilmtop’, ns’);

% Create geometry from edges
geometryFromEdges(thermalmodel, g);

% Create a larger figure window
figure(‘Position’, [100, 100, 1200, 800]);

% Figure 1
% Plot with edge labels
pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");

%Figure 2
figure
% Adjust limits and aspect ratio for better visibility
% Overall view
pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");
xlim([-5.0e-3 5.0e-3]);
ylim([0.0e-3 0.5e-3])
daspect([50 1 1]);

% Around heater
%pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");
%xlim([-5.0e-3 5.0e-3]);
%ylim([1.1e-4 1.4e-4])
%daspect([500 1 1]);

% Properties assigned assuming corresponding order as per above

% Substrate
%thermalProperties(thermalmodel,"ThermalConductivity",0.050);

%Film
thermalProperties(thermalmodel,"ThermalConductivity",0.025);

% Boundary conditions

% Edge 1 is bottom edge of film
thermalBC(thermalmodel,"Edge",1,"Temperature",30);

% Edge 2 is top edge of film
thermalBC(thermalmodel,"Edge",2,"Temperature",30);

% Edge 3 is bottom edge of heater film- no BC for continuity
%thermalBC(thermalmodel,"Edge",3,"HeatFlux",0);

% Edge 4 is top edge of heater film- no BC for continuity
%thermalBC(thermalmodel,"Edge",4,"HeatFlux",0);

% Edge 5 is right edge of bottom film- adiabatic
thermalBC(thermalmodel,"Edge",5,"HeatFlux",0);

% Edge 6 is right edge of heater film- adiabatic
thermalBC(thermalmodel,"Edge",6,"HeatFlux",0);

% Edge 7 is right edge of top film- adiabatic
thermalBC(thermalmodel,"Edge",7,"HeatFlux",0);

% Edge 8 is left edge of bottom film- adiabatic
thermalBC(thermalmodel,"Edge",8,"HeatFlux",0);

% Edge 9 is left edge of heater film- adiabatic
thermalBC(thermalmodel,"Edge",9,"HeatFlux",0);

% Edge 10 is left edge of top film- adiabatic
thermalBC(thermalmodel,"Edge",10,"HeatFlux",0);

internalHeatSource(thermalmodel,6000,"Face",2);

mesh= generateMesh(thermalmodel);
%mesh= generateMesh(thermalmodel, ‘Hmax’, 5e-8, ‘Hmin’, 5e-9);
%mesh= generateMesh(thermalmodel, ‘Hmax’, 1e-6, ‘Hmin’, 1e-7);

% Figure 5
figure
pdemesh(thermalmodel)

pdeplot(mesh)

%xlim([-1.0e-3 1.0e-3]);
%ylim([-2.0e-3 1.0e-3])
%daspect([1 1 1]);

xlim([-5.0e-3 5.0e-3]);
ylim([0.0e-3 0.5e-3])
daspect([50 1 1]);

title("Mesh with Refined Elements");

%Warning: Found elements with poor shape quality.

thermalresults = solve(thermalmodel)
T= thermalresults.Temperature;
msh= thermalresults.Mesh;

[qx,qy] = evaluateHeatFlux(thermalresults);

%figure
%pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), …
% "Contour","on",…
% "FlowData",[qx(:,end),qy(:,end)], …
%"ColorMap","jet")
%title("Temperature (C) & Heat flux fields (W/m2)")
%xlabel("Width (m)")
%ylabel("Height (m)")

% Figure 6
figure
pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), …
"Contour","on",…
"ColorMap","jet")
title("Temperature field (C)")
xlabel("Width (m)")
ylabel("Height (m)")

xlim([-1.0e-3 1.0e-3]);
ylim([-2.0e-3 1.0e-3])
daspect([1 1 1]);

% Figure 7
figure
pdeplot(thermalmodel,"FlowData",[qx(:,end), qy(:,end)])
title("Heat flux field (W/m2)")
xlabel("Width (m)")
ylabel("Height (m)")

xlim([-1.0e-3 1.0e-3]);
ylim([-2.0e-3 1.0e-3])
daspect([1 1 1]);

%isp(qx(:,end))
%disp(qy(:,end))

Qn1= evaluateHeatRate(thermalresults, "Edge",1)
Qn2= evaluateHeatRate(thermalresults, "Edge",2)
Qn3= evaluateHeatRate(thermalresults, "Edge",3)
Qn4= evaluateHeatRate(thermalresults, "Edge",4)
Qn5= evaluateHeatRate(thermalresults, "Edge",5)
Qn6= evaluateHeatRate(thermalresults, "Edge",6)
Qn7= evaluateHeatRate(thermalresults, "Edge",7)
Qn8= evaluateHeatRate(thermalresults, "Edge",8)
Qn9= evaluateHeatRate(thermalresults, "Edge",9)
Qn10= evaluateHeatRate(thermalresults, "Edge",10)

Edgetotalwatts= Qn1+Qn2+Qn4+Qn5+Qn6+Qn7+Qn8+Qn9+Qn10I have a steady state heat transfer problem trying to model a thin film heater wedged between two significantly thicker slabs. The slab thicknesses may be as much as three orders of magnitude larger than the heater. I want temperature resolution within the slabs and a temperature of the heater but I don’t need temperature resolution within the film per se. I tried treating the entire thin film as dissipating energy using an internal heat source command of this form [internalHeatSource(thermalmodel,6000,"Face",2)]. My checks on model validity included looking for higher temperature within the heater compared to the two bounding slabs using a command of this form [pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), … "Contour","on",… "ColorMap","jet")]. I also checked energy conservation in the steady state by comparing the known amount of power dissipated within the heater to the energy rate leaving the boundaries of the system using commands of the form [Qn1= evaluateHeatRate(thermalresults, "Edge",1)] for each edge and summing up these energy rates..
This approach was not successful. I did not see significant temperature increases that I would expect for heater power levels that I used. Conservation of energy was not observed. I plotted the mesh results and in some or all cases would get strange looking results. In some or all cases Matlab returned the following feedback: “Warning: Found elements with poor shape quality.” I tried manipulating the mesh size to get finer resolution using commands of the form [mesh= generateMesh(thermalmodel, ‘Hmax’, 1e-6, ‘Hmin’, 1e-7)]. I was not successful.
I feel like the problem is related to the mismatch in mesh resolution that may be required in the heater compared to the larger slabs, but I don’t know if that’s true and I don’t know how to handle it if that is true.
The code is included below.
Any help would be greatly appreciated.
john

%% Create a thermal model for transient analysis
thermalmodel = createpde("thermal","steadystate");

% Define the geometry components
Rfilmbottom = [3,4,[-5.0, 5.0, 5.0, -5.0, 0.0, 0.0, 0.124375,0.124375]/1000]’;
Rfilmhtr = [3,4,[-5.0, 5.0, 5.0, -5.0, 0.124375,0.124375,0.125625,0.125625]/1000]’;
Rfilmtop= [3,4,[-5.0, 5.0, 5.0, -5.0, 0.125625,0.125625, 0.25, 0.25]/1000]’;

% Combine the geometry data
gdm = [ Rfilmbottom Rfilmhtr Rfilmtop ];
ns = char(‘Rfilmbottom’, ‘Rfilmhtr’, ‘Rfilmtop’);
g = decsg(gdm, ‘Rfilmbottom + Rfilmhtr + Rfilmtop’, ns’);

% Create geometry from edges
geometryFromEdges(thermalmodel, g);

% Create a larger figure window
figure(‘Position’, [100, 100, 1200, 800]);

% Figure 1
% Plot with edge labels
pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");

%Figure 2
figure
% Adjust limits and aspect ratio for better visibility
% Overall view
pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");
xlim([-5.0e-3 5.0e-3]);
ylim([0.0e-3 0.5e-3])
daspect([50 1 1]);

% Around heater
%pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");
%xlim([-5.0e-3 5.0e-3]);
%ylim([1.1e-4 1.4e-4])
%daspect([500 1 1]);

% Properties assigned assuming corresponding order as per above

% Substrate
%thermalProperties(thermalmodel,"ThermalConductivity",0.050);

%Film
thermalProperties(thermalmodel,"ThermalConductivity",0.025);

% Boundary conditions

% Edge 1 is bottom edge of film
thermalBC(thermalmodel,"Edge",1,"Temperature",30);

% Edge 2 is top edge of film
thermalBC(thermalmodel,"Edge",2,"Temperature",30);

% Edge 3 is bottom edge of heater film- no BC for continuity
%thermalBC(thermalmodel,"Edge",3,"HeatFlux",0);

% Edge 4 is top edge of heater film- no BC for continuity
%thermalBC(thermalmodel,"Edge",4,"HeatFlux",0);

% Edge 5 is right edge of bottom film- adiabatic
thermalBC(thermalmodel,"Edge",5,"HeatFlux",0);

% Edge 6 is right edge of heater film- adiabatic
thermalBC(thermalmodel,"Edge",6,"HeatFlux",0);

% Edge 7 is right edge of top film- adiabatic
thermalBC(thermalmodel,"Edge",7,"HeatFlux",0);

% Edge 8 is left edge of bottom film- adiabatic
thermalBC(thermalmodel,"Edge",8,"HeatFlux",0);

% Edge 9 is left edge of heater film- adiabatic
thermalBC(thermalmodel,"Edge",9,"HeatFlux",0);

% Edge 10 is left edge of top film- adiabatic
thermalBC(thermalmodel,"Edge",10,"HeatFlux",0);

internalHeatSource(thermalmodel,6000,"Face",2);

mesh= generateMesh(thermalmodel);
%mesh= generateMesh(thermalmodel, ‘Hmax’, 5e-8, ‘Hmin’, 5e-9);
%mesh= generateMesh(thermalmodel, ‘Hmax’, 1e-6, ‘Hmin’, 1e-7);

% Figure 5
figure
pdemesh(thermalmodel)

pdeplot(mesh)

%xlim([-1.0e-3 1.0e-3]);
%ylim([-2.0e-3 1.0e-3])
%daspect([1 1 1]);

xlim([-5.0e-3 5.0e-3]);
ylim([0.0e-3 0.5e-3])
daspect([50 1 1]);

title("Mesh with Refined Elements");

%Warning: Found elements with poor shape quality.

thermalresults = solve(thermalmodel)
T= thermalresults.Temperature;
msh= thermalresults.Mesh;

[qx,qy] = evaluateHeatFlux(thermalresults);

%figure
%pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), …
% "Contour","on",…
% "FlowData",[qx(:,end),qy(:,end)], …
%"ColorMap","jet")
%title("Temperature (C) & Heat flux fields (W/m2)")
%xlabel("Width (m)")
%ylabel("Height (m)")

% Figure 6
figure
pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), …
"Contour","on",…
"ColorMap","jet")
title("Temperature field (C)")
xlabel("Width (m)")
ylabel("Height (m)")

xlim([-1.0e-3 1.0e-3]);
ylim([-2.0e-3 1.0e-3])
daspect([1 1 1]);

% Figure 7
figure
pdeplot(thermalmodel,"FlowData",[qx(:,end), qy(:,end)])
title("Heat flux field (W/m2)")
xlabel("Width (m)")
ylabel("Height (m)")

xlim([-1.0e-3 1.0e-3]);
ylim([-2.0e-3 1.0e-3])
daspect([1 1 1]);

%isp(qx(:,end))
%disp(qy(:,end))

Qn1= evaluateHeatRate(thermalresults, "Edge",1)
Qn2= evaluateHeatRate(thermalresults, "Edge",2)
Qn3= evaluateHeatRate(thermalresults, "Edge",3)
Qn4= evaluateHeatRate(thermalresults, "Edge",4)
Qn5= evaluateHeatRate(thermalresults, "Edge",5)
Qn6= evaluateHeatRate(thermalresults, "Edge",6)
Qn7= evaluateHeatRate(thermalresults, "Edge",7)
Qn8= evaluateHeatRate(thermalresults, "Edge",8)
Qn9= evaluateHeatRate(thermalresults, "Edge",9)
Qn10= evaluateHeatRate(thermalresults, "Edge",10)

Edgetotalwatts= Qn1+Qn2+Qn4+Qn5+Qn6+Qn7+Qn8+Qn9+Qn10 I have a steady state heat transfer problem trying to model a thin film heater wedged between two significantly thicker slabs. The slab thicknesses may be as much as three orders of magnitude larger than the heater. I want temperature resolution within the slabs and a temperature of the heater but I don’t need temperature resolution within the film per se. I tried treating the entire thin film as dissipating energy using an internal heat source command of this form [internalHeatSource(thermalmodel,6000,"Face",2)]. My checks on model validity included looking for higher temperature within the heater compared to the two bounding slabs using a command of this form [pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), … "Contour","on",… "ColorMap","jet")]. I also checked energy conservation in the steady state by comparing the known amount of power dissipated within the heater to the energy rate leaving the boundaries of the system using commands of the form [Qn1= evaluateHeatRate(thermalresults, "Edge",1)] for each edge and summing up these energy rates..
This approach was not successful. I did not see significant temperature increases that I would expect for heater power levels that I used. Conservation of energy was not observed. I plotted the mesh results and in some or all cases would get strange looking results. In some or all cases Matlab returned the following feedback: “Warning: Found elements with poor shape quality.” I tried manipulating the mesh size to get finer resolution using commands of the form [mesh= generateMesh(thermalmodel, ‘Hmax’, 1e-6, ‘Hmin’, 1e-7)]. I was not successful.
I feel like the problem is related to the mismatch in mesh resolution that may be required in the heater compared to the larger slabs, but I don’t know if that’s true and I don’t know how to handle it if that is true.
The code is included below.
Any help would be greatly appreciated.
john

%% Create a thermal model for transient analysis
thermalmodel = createpde("thermal","steadystate");

% Define the geometry components
Rfilmbottom = [3,4,[-5.0, 5.0, 5.0, -5.0, 0.0, 0.0, 0.124375,0.124375]/1000]’;
Rfilmhtr = [3,4,[-5.0, 5.0, 5.0, -5.0, 0.124375,0.124375,0.125625,0.125625]/1000]’;
Rfilmtop= [3,4,[-5.0, 5.0, 5.0, -5.0, 0.125625,0.125625, 0.25, 0.25]/1000]’;

% Combine the geometry data
gdm = [ Rfilmbottom Rfilmhtr Rfilmtop ];
ns = char(‘Rfilmbottom’, ‘Rfilmhtr’, ‘Rfilmtop’);
g = decsg(gdm, ‘Rfilmbottom + Rfilmhtr + Rfilmtop’, ns’);

% Create geometry from edges
geometryFromEdges(thermalmodel, g);

% Create a larger figure window
figure(‘Position’, [100, 100, 1200, 800]);

% Figure 1
% Plot with edge labels
pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");

%Figure 2
figure
% Adjust limits and aspect ratio for better visibility
% Overall view
pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");
xlim([-5.0e-3 5.0e-3]);
ylim([0.0e-3 0.5e-3])
daspect([50 1 1]);

% Around heater
%pdegplot(thermalmodel, "EdgeLabels", "on", "FaceLabels", "on");
%xlim([-5.0e-3 5.0e-3]);
%ylim([1.1e-4 1.4e-4])
%daspect([500 1 1]);

% Properties assigned assuming corresponding order as per above

% Substrate
%thermalProperties(thermalmodel,"ThermalConductivity",0.050);

%Film
thermalProperties(thermalmodel,"ThermalConductivity",0.025);

% Boundary conditions

% Edge 1 is bottom edge of film
thermalBC(thermalmodel,"Edge",1,"Temperature",30);

% Edge 2 is top edge of film
thermalBC(thermalmodel,"Edge",2,"Temperature",30);

% Edge 3 is bottom edge of heater film- no BC for continuity
%thermalBC(thermalmodel,"Edge",3,"HeatFlux",0);

% Edge 4 is top edge of heater film- no BC for continuity
%thermalBC(thermalmodel,"Edge",4,"HeatFlux",0);

% Edge 5 is right edge of bottom film- adiabatic
thermalBC(thermalmodel,"Edge",5,"HeatFlux",0);

% Edge 6 is right edge of heater film- adiabatic
thermalBC(thermalmodel,"Edge",6,"HeatFlux",0);

% Edge 7 is right edge of top film- adiabatic
thermalBC(thermalmodel,"Edge",7,"HeatFlux",0);

% Edge 8 is left edge of bottom film- adiabatic
thermalBC(thermalmodel,"Edge",8,"HeatFlux",0);

% Edge 9 is left edge of heater film- adiabatic
thermalBC(thermalmodel,"Edge",9,"HeatFlux",0);

% Edge 10 is left edge of top film- adiabatic
thermalBC(thermalmodel,"Edge",10,"HeatFlux",0);

internalHeatSource(thermalmodel,6000,"Face",2);

mesh= generateMesh(thermalmodel);
%mesh= generateMesh(thermalmodel, ‘Hmax’, 5e-8, ‘Hmin’, 5e-9);
%mesh= generateMesh(thermalmodel, ‘Hmax’, 1e-6, ‘Hmin’, 1e-7);

% Figure 5
figure
pdemesh(thermalmodel)

pdeplot(mesh)

%xlim([-1.0e-3 1.0e-3]);
%ylim([-2.0e-3 1.0e-3])
%daspect([1 1 1]);

xlim([-5.0e-3 5.0e-3]);
ylim([0.0e-3 0.5e-3])
daspect([50 1 1]);

title("Mesh with Refined Elements");

%Warning: Found elements with poor shape quality.

thermalresults = solve(thermalmodel)
T= thermalresults.Temperature;
msh= thermalresults.Mesh;

[qx,qy] = evaluateHeatFlux(thermalresults);

%figure
%pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), …
% "Contour","on",…
% "FlowData",[qx(:,end),qy(:,end)], …
%"ColorMap","jet")
%title("Temperature (C) & Heat flux fields (W/m2)")
%xlabel("Width (m)")
%ylabel("Height (m)")

% Figure 6
figure
pdeplot(thermalmodel,"XYData",thermalresults.Temperature(:,end), …
"Contour","on",…
"ColorMap","jet")
title("Temperature field (C)")
xlabel("Width (m)")
ylabel("Height (m)")

xlim([-1.0e-3 1.0e-3]);
ylim([-2.0e-3 1.0e-3])
daspect([1 1 1]);

% Figure 7
figure
pdeplot(thermalmodel,"FlowData",[qx(:,end), qy(:,end)])
title("Heat flux field (W/m2)")
xlabel("Width (m)")
ylabel("Height (m)")

xlim([-1.0e-3 1.0e-3]);
ylim([-2.0e-3 1.0e-3])
daspect([1 1 1]);

%isp(qx(:,end))
%disp(qy(:,end))

Qn1= evaluateHeatRate(thermalresults, "Edge",1)
Qn2= evaluateHeatRate(thermalresults, "Edge",2)
Qn3= evaluateHeatRate(thermalresults, "Edge",3)
Qn4= evaluateHeatRate(thermalresults, "Edge",4)
Qn5= evaluateHeatRate(thermalresults, "Edge",5)
Qn6= evaluateHeatRate(thermalresults, "Edge",6)
Qn7= evaluateHeatRate(thermalresults, "Edge",7)
Qn8= evaluateHeatRate(thermalresults, "Edge",8)
Qn9= evaluateHeatRate(thermalresults, "Edge",9)
Qn10= evaluateHeatRate(thermalresults, "Edge",10)

Edgetotalwatts= Qn1+Qn2+Qn4+Qn5+Qn6+Qn7+Qn8+Qn9+Qn10 heat transfer, thin heater, mesh problems MATLAB Answers — New Questions

​

I would like to prototype software-defined radio (SDR) systems using MATLAB & Simulink using the Xilinx Virtex-6 and Spartan-6 FPGA boards. Is there an available Support Package for Xilinx FPGA Radio?I would like to prototype software-defined radio (SDR) systems using MATLAB & Simulink using the Xilinx Virtex-6 and Spartan-6 FPGA boards. Is there an available Support Package for Xilinx FPGA Radio? I would like to prototype software-defined radio (SDR) systems using MATLAB & Simulink using the Xilinx Virtex-6 and Spartan-6 FPGA boards. Is there an available Support Package for Xilinx FPGA Radio? xilinx, fpga, support, package, communications, system, toolbox, zynq, sdr MATLAB Answers — New Questions

​

Hello everyone, I’m very new to MATLAB. Could you help me with this question please?
How to write codes that calculate e^x by serializing the x value entered from the keyboard into a series equal to the number of terms (N) entered from the keyboard?Hello everyone, I’m very new to MATLAB. Could you help me with this question please?
How to write codes that calculate e^x by serializing the x value entered from the keyboard into a series equal to the number of terms (N) entered from the keyboard? Hello everyone, I’m very new to MATLAB. Could you help me with this question please?
How to write codes that calculate e^x by serializing the x value entered from the keyboard into a series equal to the number of terms (N) entered from the keyboard? #e^x #maclaurin #taylorseries MATLAB Answers — New Questions

​

I have a parameter in a model reference that I would like to tune while the simulation is running in accelerator mode. How do I do this?I have a parameter in a model reference that I would like to tune while the simulation is running in accelerator mode. How do I do this? I have a parameter in a model reference that I would like to tune while the simulation is running in accelerator mode. How do I do this?  MATLAB Answers — New Questions

​

Hi there,
I have been getting different outputs from Matlab to JASP and SPSS, both of which agree that Matlab is wrong, which I think is due to the intercept being included as an interaction term in the model spec.

I have run a 2b*6w*2w Mixed ANOVA in Matlab. (Between-2=Consciousness, Within-6 = Conguency, Within-2 =Target) However I am getting the intercept as an interaction term in my ranova table:
Matlab Output

The factors Congruency and Target should be reported only as a main effects, not as an interaction with the intercept.

I’m including a JASP output so to compare the difference between the two softwares.

Could you take a look at my code, I think I’ve done something wrong in the Wilkinson notation of the model specification, and tell me how I stop it using the intercept as an interaction term.

OverallDataTable=readtable(‘Data.xlsx’)

WithinDesign=table(categorical({‘Related Hand’,’Related Foot’,’Unrelated Hand’,’Unrelated Foot’, …
‘Neutral Hand’,’Neutral Foot’,’Related Hand’,’Related Foot’,’Unrelated Hand’,’Unrelated Foot’, …
‘Neutral Hand’,’Neutral Foot’}’),categorical({‘BP’,’BP’,’BP’,’BP’,’BP’,’BP’,’O’,’O’,’O’,’O’,’O’,’O’}’), …
‘VariableNames’,{‘Congruency’,’Target’});

OverallMixedANOVA=fitrm(OverallDataTable,’BP_RH_Con-O_NF_Sub~Consciousness’,WithinDesign=WithinDesign);

OverallMixedANOVATable=ranova(OverallMixedANOVA,"WithinModel",’Congruency*Target’)

Thanks for any and all advice

DanHi there,
I have been getting different outputs from Matlab to JASP and SPSS, both of which agree that Matlab is wrong, which I think is due to the intercept being included as an interaction term in the model spec.

I have run a 2b*6w*2w Mixed ANOVA in Matlab. (Between-2=Consciousness, Within-6 = Conguency, Within-2 =Target) However I am getting the intercept as an interaction term in my ranova table:
Matlab Output

The factors Congruency and Target should be reported only as a main effects, not as an interaction with the intercept.

I’m including a JASP output so to compare the difference between the two softwares.

Could you take a look at my code, I think I’ve done something wrong in the Wilkinson notation of the model specification, and tell me how I stop it using the intercept as an interaction term.

OverallDataTable=readtable(‘Data.xlsx’)

WithinDesign=table(categorical({‘Related Hand’,’Related Foot’,’Unrelated Hand’,’Unrelated Foot’, …
‘Neutral Hand’,’Neutral Foot’,’Related Hand’,’Related Foot’,’Unrelated Hand’,’Unrelated Foot’, …
‘Neutral Hand’,’Neutral Foot’}’),categorical({‘BP’,’BP’,’BP’,’BP’,’BP’,’BP’,’O’,’O’,’O’,’O’,’O’,’O’}’), …
‘VariableNames’,{‘Congruency’,’Target’});

OverallMixedANOVA=fitrm(OverallDataTable,’BP_RH_Con-O_NF_Sub~Consciousness’,WithinDesign=WithinDesign);

OverallMixedANOVATable=ranova(OverallMixedANOVA,"WithinModel",’Congruency*Target’)

Thanks for any and all advice

Dan Hi there,
I have been getting different outputs from Matlab to JASP and SPSS, both of which agree that Matlab is wrong, which I think is due to the intercept being included as an interaction term in the model spec.

I have run a 2b*6w*2w Mixed ANOVA in Matlab. (Between-2=Consciousness, Within-6 = Conguency, Within-2 =Target) However I am getting the intercept as an interaction term in my ranova table:
Matlab Output

The factors Congruency and Target should be reported only as a main effects, not as an interaction with the intercept.

I’m including a JASP output so to compare the difference between the two softwares.

Could you take a look at my code, I think I’ve done something wrong in the Wilkinson notation of the model specification, and tell me how I stop it using the intercept as an interaction term.

OverallDataTable=readtable(‘Data.xlsx’)

WithinDesign=table(categorical({‘Related Hand’,’Related Foot’,’Unrelated Hand’,’Unrelated Foot’, …
‘Neutral Hand’,’Neutral Foot’,’Related Hand’,’Related Foot’,’Unrelated Hand’,’Unrelated Foot’, …
‘Neutral Hand’,’Neutral Foot’}’),categorical({‘BP’,’BP’,’BP’,’BP’,’BP’,’BP’,’O’,’O’,’O’,’O’,’O’,’O’}’), …
‘VariableNames’,{‘Congruency’,’Target’});

OverallMixedANOVA=fitrm(OverallDataTable,’BP_RH_Con-O_NF_Sub~Consciousness’,WithinDesign=WithinDesign);

OverallMixedANOVATable=ranova(OverallMixedANOVA,"WithinModel",’Congruency*Target’)

Thanks for any and all advice

Dan fitrm, statistics, mixed anova, ranova, wilkinson notation MATLAB Answers — New Questions

​

Hi,
The full error is:
"Error: File: builtin.m Line: 1 Column: 24
Invalid text character. Check for unsupported symbol, invisible character, or pasting of non-ASCII characters."
I get it whenever I try to do anything, including close Matlab. The file name changes slightly depending on what I do: for example, when I try to open a new file, I get the following error:
"Error in matlab.unittest.internal.ui.toolstrip.getFileInfoForToolstrip (line 8)
isClassBasedTest = false;",
preceded by the above text with "false.m" as the file name.

I think the problem started when I tried to set the path? I got the following error:
"Unable to run P-code file. The file header is damaged.
Error in pathtool (line 13)
isDesktop = desktop(‘-inuse’);"
But I’m not sure. Any help that you can give me would be appreciated!Hi,
The full error is:
"Error: File: builtin.m Line: 1 Column: 24
Invalid text character. Check for unsupported symbol, invisible character, or pasting of non-ASCII characters."
I get it whenever I try to do anything, including close Matlab. The file name changes slightly depending on what I do: for example, when I try to open a new file, I get the following error:
"Error in matlab.unittest.internal.ui.toolstrip.getFileInfoForToolstrip (line 8)
isClassBasedTest = false;",
preceded by the above text with "false.m" as the file name.

I think the problem started when I tried to set the path? I got the following error:
"Unable to run P-code file. The file header is damaged.
Error in pathtool (line 13)
isDesktop = desktop(‘-inuse’);"
But I’m not sure. Any help that you can give me would be appreciated! Hi,
The full error is:
"Error: File: builtin.m Line: 1 Column: 24
Invalid text character. Check for unsupported symbol, invisible character, or pasting of non-ASCII characters."
I get it whenever I try to do anything, including close Matlab. The file name changes slightly depending on what I do: for example, when I try to open a new file, I get the following error:
"Error in matlab.unittest.internal.ui.toolstrip.getFileInfoForToolstrip (line 8)
isClassBasedTest = false;",
preceded by the above text with "false.m" as the file name.

I think the problem started when I tried to set the path? I got the following error:
"Unable to run P-code file. The file header is damaged.
Error in pathtool (line 13)
isDesktop = desktop(‘-inuse’);"
But I’m not sure. Any help that you can give me would be appreciated! matlab, set path MATLAB Answers — New Questions

​

I have a pretty large script that analyzes a somewhat large data file (~5gbs). My script works fine when I first run it but then when I go to run a codeblock or even just try to get the output from a single simple variable in the command window AFTER I loaded everything into my workspace, MATLAB will stall for a minute or even more before starting to run whatever command I gave it. I’ve been monitoring my PC resources and it doesn’t seem like I am running out of RAM or anything (working with 64gbs). I even went through and cleared many larger variables that were not needed for later parts of the script and the problem persists. I do not recieve any errors, it is just very slow to do simple things. Once it starts executing the command, it runs at the expected speed (I’ve verfied with some manual progress bars I coded in).
The data that I load is from a single .MAT file which has a structure in it with all of my data. I’ve also ran this script on 3 other PCs and had the same issue.I have a pretty large script that analyzes a somewhat large data file (~5gbs). My script works fine when I first run it but then when I go to run a codeblock or even just try to get the output from a single simple variable in the command window AFTER I loaded everything into my workspace, MATLAB will stall for a minute or even more before starting to run whatever command I gave it. I’ve been monitoring my PC resources and it doesn’t seem like I am running out of RAM or anything (working with 64gbs). I even went through and cleared many larger variables that were not needed for later parts of the script and the problem persists. I do not recieve any errors, it is just very slow to do simple things. Once it starts executing the command, it runs at the expected speed (I’ve verfied with some manual progress bars I coded in).
The data that I load is from a single .MAT file which has a structure in it with all of my data. I’ve also ran this script on 3 other PCs and had the same issue. I have a pretty large script that analyzes a somewhat large data file (~5gbs). My script works fine when I first run it but then when I go to run a codeblock or even just try to get the output from a single simple variable in the command window AFTER I loaded everything into my workspace, MATLAB will stall for a minute or even more before starting to run whatever command I gave it. I’ve been monitoring my PC resources and it doesn’t seem like I am running out of RAM or anything (working with 64gbs). I even went through and cleared many larger variables that were not needed for later parts of the script and the problem persists. I do not recieve any errors, it is just very slow to do simple things. Once it starts executing the command, it runs at the expected speed (I’ve verfied with some manual progress bars I coded in).
The data that I load is from a single .MAT file which has a structure in it with all of my data. I’ve also ran this script on 3 other PCs and had the same issue. structures MATLAB Answers — New Questions

​

I’m working on a MATLAB plot where I want to achieve the following functionality:
Change the color of a line when it’s clicked.
Enable panning by dragging while clicking off the line.
Enable zooming using the scroll wheel.
Here’s the simple code I’m using:

% Plotting a simple line
h = plot([0,1],[0,1],’LineWidth’,3);

% Setting the ButtonDownFcn to change line color (placeholder code)
h.ButtonDownFcn = @(~,~) disp(h);

In this example, clicking on the line displays its handle h in the Command Window, which is a placeholder for my actual code to change the line color.

The Problem:

Assigning a ButtonDownFcn to the line object seems to override MATLAB’s built-in pan and zoom functionalities. After setting the ButtonDownFcn, I’m unable to pan by clicking and dragging off the line, and the scroll wheel zoom no longer works. It appears that the custom callback interferes with the default interactive behaviors, even when I’m not interacting directly with the line.

My Questions:

Why does setting the ButtonDownFcn on the line object disable panning and zooming, even when interacting off the line?
Is there a way to have both the custom click behavior (changing the line color when clicked) and retain the default pan and zoom functionalities when interacting elsewhere on the plot?I’m working on a MATLAB plot where I want to achieve the following functionality:
Change the color of a line when it’s clicked.
Enable panning by dragging while clicking off the line.
Enable zooming using the scroll wheel.
Here’s the simple code I’m using:

% Plotting a simple line
h = plot([0,1],[0,1],’LineWidth’,3);

% Setting the ButtonDownFcn to change line color (placeholder code)
h.ButtonDownFcn = @(~,~) disp(h);

In this example, clicking on the line displays its handle h in the Command Window, which is a placeholder for my actual code to change the line color.

The Problem:

Assigning a ButtonDownFcn to the line object seems to override MATLAB’s built-in pan and zoom functionalities. After setting the ButtonDownFcn, I’m unable to pan by clicking and dragging off the line, and the scroll wheel zoom no longer works. It appears that the custom callback interferes with the default interactive behaviors, even when I’m not interacting directly with the line.

My Questions:

Why does setting the ButtonDownFcn on the line object disable panning and zooming, even when interacting off the line?
Is there a way to have both the custom click behavior (changing the line color when clicked) and retain the default pan and zoom functionalities when interacting elsewhere on the plot? I’m working on a MATLAB plot where I want to achieve the following functionality:
Change the color of a line when it’s clicked.
Enable panning by dragging while clicking off the line.
Enable zooming using the scroll wheel.
Here’s the simple code I’m using:

% Plotting a simple line
h = plot([0,1],[0,1],’LineWidth’,3);

% Setting the ButtonDownFcn to change line color (placeholder code)
h.ButtonDownFcn = @(~,~) disp(h);

In this example, clicking on the line displays its handle h in the Command Window, which is a placeholder for my actual code to change the line color.

The Problem:

Assigning a ButtonDownFcn to the line object seems to override MATLAB’s built-in pan and zoom functionalities. After setting the ButtonDownFcn, I’m unable to pan by clicking and dragging off the line, and the scroll wheel zoom no longer works. It appears that the custom callback interferes with the default interactive behaviors, even when I’m not interacting directly with the line.

My Questions:

Why does setting the ButtonDownFcn on the line object disable panning and zooming, even when interacting off the line?
Is there a way to have both the custom click behavior (changing the line color when clicked) and retain the default pan and zoom functionalities when interacting elsewhere on the plot? callback, plot MATLAB Answers — New Questions

​

function [y, signalT, notch_intensity] = signalfunccycle(Gridlength, startfreq, endfreq, notchstart1, notchend1, datalength)
close all;
clear all;
PIx2 = 2 * pi; % used to calculate the frequency values of the output signal
fstart = 0;
Gridlength = 200e-3;
startfreq = 40e3;
endfreq = 200e3;
stimeStep = 5e-8;
datalength = round(Gridlength / stimeStep);

notchstart1 = 116000;
notchend1 = 119000;
dataN = datalength;
Stimelength = Gridlength;
fend = 10^6;
Fs = dataN / Gridlength;
NFFT = dataN; % Number of points used in the FFT analysis
fstep = (Fs – fstart) / NFFT; % Point step

T = 1 / Fs;
Stime = linspace(0, Stimelength, dataN);
t = (0:dataN-1) * T; % Length of time
signal = zeros(1, length(t));
f1 = notchstart1;
f2 = notchend1;
f3 = notchstart1;
f4 = notchend1;

for freq = startfreq:0.5*10^3:endfreq
if (f1 <= freq && freq <= f2) || (f3 <= freq && freq <= f4)
continue;
end
signal = sin(PIx2 * freq * t + pi) + signal;
end

%% Time-domain signal graphs
figure(1);
plot(Stime, signal);
grid on;
title(‘signal’);
xlabel(‘time’);
ylabel(‘amplitude’);

%% FFT
fMax = NFFT / 2 + 1; % The maximum frequency to plot
signalFFT0 = fft(signal, NFFT); % FFT of the signal
signalFFT_phase = angle(fft(signal, NFFT)); % Phase of the FFT result
signalFFT = abs(fft(signal, NFFT)); % Absolute value (modulus) of the FFT result
signalFFTShow = 2 * signalFFT / dataN; % Actual amplitude values
signalFFTShow(1) = signalFFTShow(1) / 2; % DC component
f = Fs / 2 * linspace(0, 1, fMax); % Actual frequencies

% Plot FFT of the signal
figure(2);
plot(f, signalFFTShow(1:fMax), ‘r’, ‘LineWidth’, 3); % Red color, line width set to 3
grid on; % Display grid lines
title(‘FFT Signal’);
xlabel(‘Frequency’);
ylabel(‘Amplitude’);

%% Phase function modulation
signalstart = startfreq;
signalend = endfreq;
Istart = round(signalstart / fstep);
Iend = round(signalend / fstep);
Itotal = Iend – Istart;

N1start1 = round(f1);
N1end1 = round(f2);
N1start2 = round(f3);
N1end2 = round(f4);
signalphase(1:NFFT) = 0;

for I = 1:NFFT
if (I <= N1end1 && I >= N1start1) || (I <= N1end2 && I >= N1start2)
signalphase(I) = 0;
else
signalphase(I) = pi / 2 * (I – Istart) + 1 * pi / Itotal / 2 * (I – Istart)^2;
end
end

signalF = 1/2 * signalFFTShow .* exp(1i * (signalphase + signalFFT_phase));
signalT = ifft(NFFT * signalF);
t1 = (0:NFFT-1) * T;

%% Phase spectrogram
x1 = 1:NFFT;
figure(3);
plot(x1, signalphase)
xlabel(‘frequency’);
ylabel(‘Phase’);
title(‘Phase spectrum’);

%% Signal waveform after modulation
signaltime = (0:NFFT-1) * T;
figure(4);
plot(signaltime, real(signalT));
xlabel(‘time’);
ylabel(‘amplitude’);
signalT_real = real(signalT);
savefile = ‘signalT_real.txt’;
x = Stime’;
y = signalT_real’;
save(savefile, ‘x’, ‘y’, ‘-ASCII’)

%% STFT Calculation
NN = 10;
win_size = round(datalength / NN);
overlap = 0;
fft_size = win_size * 2 + 1;
stft_result = [];

for i = 1:win_size*(1-overlap):length(signalT)-win_size
window = signalT(i:i+win_size-1);
fft_result = fft(window, fft_size+1);
stft_result = [stft_result; abs(fft_result(1:fft_size))];
end

time_axis = linspace(0, Gridlength, size(stft_result, 1));
freq_axis = linspace(fstart, fend, fft_size);

% Extract intensity for notch frequencies (116 to 119 kHz)
notch_indices = find(freq_axis >= 116000 & freq_axis <= 119000);
notch_intensity = mean(stft_result(:, notch_indices), 2);

% Plot the STFT spectrogram
figure(5);
imagesc(time_axis * 1000, freq_axis / 1000, stft_result.’);
axis([Gridlength * 950 Gridlength * 1050 110e3 / 1000 125e3 / 1000]);
colorbar;
ylabel(‘Frequency (kHz)’);
xlabel(‘Time (ms)’);

% % Plot the intensity over time for the notch frequencies
% figure(6);
% plot(time_axis * 1000, notch_intensity, ‘LineWidth’, 2);
% xlabel(‘Time (ms)’);
% ylabel(‘Intensity’);
% title(‘Intensity at Notch Frequencies (116 to 119 kHz)’);
% grid on;
endfunction [y, signalT, notch_intensity] = signalfunccycle(Gridlength, startfreq, endfreq, notchstart1, notchend1, datalength)
close all;
clear all;
PIx2 = 2 * pi; % used to calculate the frequency values of the output signal
fstart = 0;
Gridlength = 200e-3;
startfreq = 40e3;
endfreq = 200e3;
stimeStep = 5e-8;
datalength = round(Gridlength / stimeStep);

notchstart1 = 116000;
notchend1 = 119000;
dataN = datalength;
Stimelength = Gridlength;
fend = 10^6;
Fs = dataN / Gridlength;
NFFT = dataN; % Number of points used in the FFT analysis
fstep = (Fs – fstart) / NFFT; % Point step

T = 1 / Fs;
Stime = linspace(0, Stimelength, dataN);
t = (0:dataN-1) * T; % Length of time
signal = zeros(1, length(t));
f1 = notchstart1;
f2 = notchend1;
f3 = notchstart1;
f4 = notchend1;

for freq = startfreq:0.5*10^3:endfreq
if (f1 <= freq && freq <= f2) || (f3 <= freq && freq <= f4)
continue;
end
signal = sin(PIx2 * freq * t + pi) + signal;
end

%% Time-domain signal graphs
figure(1);
plot(Stime, signal);
grid on;
title(‘signal’);
xlabel(‘time’);
ylabel(‘amplitude’);

%% FFT
fMax = NFFT / 2 + 1; % The maximum frequency to plot
signalFFT0 = fft(signal, NFFT); % FFT of the signal
signalFFT_phase = angle(fft(signal, NFFT)); % Phase of the FFT result
signalFFT = abs(fft(signal, NFFT)); % Absolute value (modulus) of the FFT result
signalFFTShow = 2 * signalFFT / dataN; % Actual amplitude values
signalFFTShow(1) = signalFFTShow(1) / 2; % DC component
f = Fs / 2 * linspace(0, 1, fMax); % Actual frequencies

% Plot FFT of the signal
figure(2);
plot(f, signalFFTShow(1:fMax), ‘r’, ‘LineWidth’, 3); % Red color, line width set to 3
grid on; % Display grid lines
title(‘FFT Signal’);
xlabel(‘Frequency’);
ylabel(‘Amplitude’);

%% Phase function modulation
signalstart = startfreq;
signalend = endfreq;
Istart = round(signalstart / fstep);
Iend = round(signalend / fstep);
Itotal = Iend – Istart;

N1start1 = round(f1);
N1end1 = round(f2);
N1start2 = round(f3);
N1end2 = round(f4);
signalphase(1:NFFT) = 0;

for I = 1:NFFT
if (I <= N1end1 && I >= N1start1) || (I <= N1end2 && I >= N1start2)
signalphase(I) = 0;
else
signalphase(I) = pi / 2 * (I – Istart) + 1 * pi / Itotal / 2 * (I – Istart)^2;
end
end

signalF = 1/2 * signalFFTShow .* exp(1i * (signalphase + signalFFT_phase));
signalT = ifft(NFFT * signalF);
t1 = (0:NFFT-1) * T;

%% Phase spectrogram
x1 = 1:NFFT;
figure(3);
plot(x1, signalphase)
xlabel(‘frequency’);
ylabel(‘Phase’);
title(‘Phase spectrum’);

%% Signal waveform after modulation
signaltime = (0:NFFT-1) * T;
figure(4);
plot(signaltime, real(signalT));
xlabel(‘time’);
ylabel(‘amplitude’);
signalT_real = real(signalT);
savefile = ‘signalT_real.txt’;
x = Stime’;
y = signalT_real’;
save(savefile, ‘x’, ‘y’, ‘-ASCII’)

%% STFT Calculation
NN = 10;
win_size = round(datalength / NN);
overlap = 0;
fft_size = win_size * 2 + 1;
stft_result = [];

for i = 1:win_size*(1-overlap):length(signalT)-win_size
window = signalT(i:i+win_size-1);
fft_result = fft(window, fft_size+1);
stft_result = [stft_result; abs(fft_result(1:fft_size))];
end

time_axis = linspace(0, Gridlength, size(stft_result, 1));
freq_axis = linspace(fstart, fend, fft_size);

% Extract intensity for notch frequencies (116 to 119 kHz)
notch_indices = find(freq_axis >= 116000 & freq_axis <= 119000);
notch_intensity = mean(stft_result(:, notch_indices), 2);

% Plot the STFT spectrogram
figure(5);
imagesc(time_axis * 1000, freq_axis / 1000, stft_result.’);
axis([Gridlength * 950 Gridlength * 1050 110e3 / 1000 125e3 / 1000]);
colorbar;
ylabel(‘Frequency (kHz)’);
xlabel(‘Time (ms)’);

% % Plot the intensity over time for the notch frequencies
% figure(6);
% plot(time_axis * 1000, notch_intensity, ‘LineWidth’, 2);
% xlabel(‘Time (ms)’);
% ylabel(‘Intensity’);
% title(‘Intensity at Notch Frequencies (116 to 119 kHz)’);
% grid on;
end function [y, signalT, notch_intensity] = signalfunccycle(Gridlength, startfreq, endfreq, notchstart1, notchend1, datalength)
close all;
clear all;
PIx2 = 2 * pi; % used to calculate the frequency values of the output signal
fstart = 0;
Gridlength = 200e-3;
startfreq = 40e3;
endfreq = 200e3;
stimeStep = 5e-8;
datalength = round(Gridlength / stimeStep);

notchstart1 = 116000;
notchend1 = 119000;
dataN = datalength;
Stimelength = Gridlength;
fend = 10^6;
Fs = dataN / Gridlength;
NFFT = dataN; % Number of points used in the FFT analysis
fstep = (Fs – fstart) / NFFT; % Point step

T = 1 / Fs;
Stime = linspace(0, Stimelength, dataN);
t = (0:dataN-1) * T; % Length of time
signal = zeros(1, length(t));
f1 = notchstart1;
f2 = notchend1;
f3 = notchstart1;
f4 = notchend1;

for freq = startfreq:0.5*10^3:endfreq
if (f1 <= freq && freq <= f2) || (f3 <= freq && freq <= f4)
continue;
end
signal = sin(PIx2 * freq * t + pi) + signal;
end

%% Time-domain signal graphs
figure(1);
plot(Stime, signal);
grid on;
title(‘signal’);
xlabel(‘time’);
ylabel(‘amplitude’);

%% FFT
fMax = NFFT / 2 + 1; % The maximum frequency to plot
signalFFT0 = fft(signal, NFFT); % FFT of the signal
signalFFT_phase = angle(fft(signal, NFFT)); % Phase of the FFT result
signalFFT = abs(fft(signal, NFFT)); % Absolute value (modulus) of the FFT result
signalFFTShow = 2 * signalFFT / dataN; % Actual amplitude values
signalFFTShow(1) = signalFFTShow(1) / 2; % DC component
f = Fs / 2 * linspace(0, 1, fMax); % Actual frequencies

% Plot FFT of the signal
figure(2);
plot(f, signalFFTShow(1:fMax), ‘r’, ‘LineWidth’, 3); % Red color, line width set to 3
grid on; % Display grid lines
title(‘FFT Signal’);
xlabel(‘Frequency’);
ylabel(‘Amplitude’);

%% Phase function modulation
signalstart = startfreq;
signalend = endfreq;
Istart = round(signalstart / fstep);
Iend = round(signalend / fstep);
Itotal = Iend – Istart;

N1start1 = round(f1);
N1end1 = round(f2);
N1start2 = round(f3);
N1end2 = round(f4);
signalphase(1:NFFT) = 0;

for I = 1:NFFT
if (I <= N1end1 && I >= N1start1) || (I <= N1end2 && I >= N1start2)
signalphase(I) = 0;
else
signalphase(I) = pi / 2 * (I – Istart) + 1 * pi / Itotal / 2 * (I – Istart)^2;
end
end

signalF = 1/2 * signalFFTShow .* exp(1i * (signalphase + signalFFT_phase));
signalT = ifft(NFFT * signalF);
t1 = (0:NFFT-1) * T;

%% Phase spectrogram
x1 = 1:NFFT;
figure(3);
plot(x1, signalphase)
xlabel(‘frequency’);
ylabel(‘Phase’);
title(‘Phase spectrum’);

%% Signal waveform after modulation
signaltime = (0:NFFT-1) * T;
figure(4);
plot(signaltime, real(signalT));
xlabel(‘time’);
ylabel(‘amplitude’);
signalT_real = real(signalT);
savefile = ‘signalT_real.txt’;
x = Stime’;
y = signalT_real’;
save(savefile, ‘x’, ‘y’, ‘-ASCII’)

%% STFT Calculation
NN = 10;
win_size = round(datalength / NN);
overlap = 0;
fft_size = win_size * 2 + 1;
stft_result = [];

for i = 1:win_size*(1-overlap):length(signalT)-win_size
window = signalT(i:i+win_size-1);
fft_result = fft(window, fft_size+1);
stft_result = [stft_result; abs(fft_result(1:fft_size))];
end

time_axis = linspace(0, Gridlength, size(stft_result, 1));
freq_axis = linspace(fstart, fend, fft_size);

% Extract intensity for notch frequencies (116 to 119 kHz)
notch_indices = find(freq_axis >= 116000 & freq_axis <= 119000);
notch_intensity = mean(stft_result(:, notch_indices), 2);

% Plot the STFT spectrogram
figure(5);
imagesc(time_axis * 1000, freq_axis / 1000, stft_result.’);
axis([Gridlength * 950 Gridlength * 1050 110e3 / 1000 125e3 / 1000]);
colorbar;
ylabel(‘Frequency (kHz)’);
xlabel(‘Time (ms)’);

% % Plot the intensity over time for the notch frequencies
% figure(6);
% plot(time_axis * 1000, notch_intensity, ‘LineWidth’, 2);
% xlabel(‘Time (ms)’);
% ylabel(‘Intensity’);
% title(‘Intensity at Notch Frequencies (116 to 119 kHz)’);
% grid on;
end matlab code MATLAB Answers — New Questions

​

Implement the Gauss-Seidel method for solving a system of linear equations from scratch in MATLAB and walk me through your thought process in constructing the code. Additionally, demonstrate that your implementation works by applying it to the following system. 2x +3y −z =7 , x −2y+4z =1 , 3x +y+2z =8
%Gauss seidel
A=input(‘Enter a co-efficient matrix A: ‘);
B=input(‘Enter source vector B: ‘);
P=input(‘Enter initial guess vector: ‘);
n=input(‘Enter number of iterations: ‘);
e=input(‘Enter tolerance: ‘);
N=length(B);
X=zeros(N,1);
Y=zeros(N,1); %for stopping criteria
for j=1:n
for i=1:N
x(i)=(B(i)/A(i,i))-(A(i,[1:i-1,i+1:N])*P([1:i-1,i+1:N]))/A(i,i);
P(i)=X(i);
end
fprintf(‘Iterations no %d n’ ,j)
X
if abs (Y-X)<e
break
end
Y=X;
end
This is my code is this correct? please help meImplement the Gauss-Seidel method for solving a system of linear equations from scratch in MATLAB and walk me through your thought process in constructing the code. Additionally, demonstrate that your implementation works by applying it to the following system. 2x +3y −z =7 , x −2y+4z =1 , 3x +y+2z =8
%Gauss seidel
A=input(‘Enter a co-efficient matrix A: ‘);
B=input(‘Enter source vector B: ‘);
P=input(‘Enter initial guess vector: ‘);
n=input(‘Enter number of iterations: ‘);
e=input(‘Enter tolerance: ‘);
N=length(B);
X=zeros(N,1);
Y=zeros(N,1); %for stopping criteria
for j=1:n
for i=1:N
x(i)=(B(i)/A(i,i))-(A(i,[1:i-1,i+1:N])*P([1:i-1,i+1:N]))/A(i,i);
P(i)=X(i);
end
fprintf(‘Iterations no %d n’ ,j)
X
if abs (Y-X)<e
break
end
Y=X;
end
This is my code is this correct? please help me Implement the Gauss-Seidel method for solving a system of linear equations from scratch in MATLAB and walk me through your thought process in constructing the code. Additionally, demonstrate that your implementation works by applying it to the following system. 2x +3y −z =7 , x −2y+4z =1 , 3x +y+2z =8
%Gauss seidel
A=input(‘Enter a co-efficient matrix A: ‘);
B=input(‘Enter source vector B: ‘);
P=input(‘Enter initial guess vector: ‘);
n=input(‘Enter number of iterations: ‘);
e=input(‘Enter tolerance: ‘);
N=length(B);
X=zeros(N,1);
Y=zeros(N,1); %for stopping criteria
for j=1:n
for i=1:N
x(i)=(B(i)/A(i,i))-(A(i,[1:i-1,i+1:N])*P([1:i-1,i+1:N]))/A(i,i);
P(i)=X(i);
end
fprintf(‘Iterations no %d n’ ,j)
X
if abs (Y-X)<e
break
end
Y=X;
end
This is my code is this correct? please help me @wantquickhelp MATLAB Answers — New Questions

​