Community Profile

photo

Davide Masiello


Paul Scherrer Institut (PSI)

Last seen: 3 días hace Con actividad desde 2017

Programming Languages:
MATLAB
Spoken Languages:
English, Italian
Professional Interests:
Chemical Kinetics, Fluid Dynamics, Chemical Engineering, Heat and Mass Transfer

Estadísticas

All
  • 6 Month Streak
  • MATLAB Mini Hack 2022 Participant
  • Knowledgeable Level 5
  • Pro
  • Introduction to MATLAB Master
  • Thankful Level 2
  • Knowledgeable Level 3
  • Solver
  • First Answer
  • Thankful Level 1

Ver insignias

Content Feed

Ver por

Respondida
Insert zeros in random positions following another zero
This should work without the need for a for loop. X = [2 1 0 4 5 0 778 90 0 3 88 0 77 66 0 12 23 0 45 80 0 89 67 0 34 67 0 76 3...

4 meses hace | 0

Respondida
Modelling with 2 input variables and 1 output variable
Hello - I would recommend reading the documentation of the function fit.

4 meses hace | 0

Respondida
Array indices must be positive integers or logical values.
The problem is in the convergence check step, where you try to evaluate eta at n-1 even when n = 1. Adding n > 1 as additional ...

7 meses hace | 0

Respondida
Calculating the area of overlapping polgyons
How about this? % Clear the workspace and close all figures clear,clc % Define the vertices of the polygons polygon1 = [0,...

8 meses hace | 0

Respondida
How to morph a shape into another
I guess it boils down to what mapping you choose from the points of shape 1 to the points of shape 2. Let's take your case for ...

8 meses hace | 1

| aceptada

Respondida
How do I shade the area between two curves on a plot that is generated by arrays?
Use the patch function. x1 = [1 2 3]; y1 = [2 5 8]; x2 = [1 2.5 3.5]; y2 = [1 6 3]; hold on patch([x1 flip(x2)], [y1 fli...

8 meses hace | 1

| aceptada

Respondida
Combining legend data ?
Here's an example with made up functions (your code was too slow on the browser). Is this what you're looking for? y = linspa...

8 meses hace | 0

Respondida
Is it possible to programmatically suppress figure roll-over menu (zoom, rotate, etc)?
This should work, but it also suppresses the fixed toolbar above the figure f = figure; f.ToolBar = 'none';

8 meses hace | 0

Respondida
Solve for pressure swing adsorption cycle for two gases in one bed.
Hi @Samridh Sharma, please see the modifications I made to the code and the resulting output. Please let me know if you have fu...

8 meses hace | 0

| aceptada

Respondida
Solving simultaneous set of differential equations for range of values to get desired set output
I'd recommend taking a look at the material at the following link https://ch.mathworks.com/help/matlab/math/boundary-value-prob...

11 meses hace | 0

Respondida
Draw 1 hyperboloid with equation -0.2x^2 - 0.2y^2 +z^2 =1 with colormap jet
[x,y] = meshgrid(linspace(-10,10,100)); z = sqrt(1+0.2*x.^2+0.2*y.^2); Then you can plot it in many ways. Examples figure...

12 meses hace | 0

Respondida
ODE simulation: issue with the size of arrays
In addition to my comment above, here's a solution to your problem. It's based on the fact that you can specify tspan as a vect...

12 meses hace | 0

Respondida
How to do custom equation (non linear) regression?
Assume these are your experimental data x = linspace(0,20,30); y = rand(size(x))/3+(pi*x.^(sqrt(2)/2))./(26.1-x); figure(1) ...

12 meses hace | 0

Respondida
help with curve fitting with a data function
Just needed better initial guesses, I think data = readmatrix('ottica10.txt'); M = 2; myfun = @(p,x) 1/M*exp(-p(1)*0.5*x).*ab...

12 meses hace | 0

| aceptada

Respondida
Hello! How to define on which plane lies a given point with coordinates? I have coordinates of plane and normal
Four points can be verified to belong on the same plane if the determinant of the matrix is equal to zero. I'd translate tha...

12 meses hace | 0

Respondida
Interpolation of given data in 2d domain
Example xdata = linspace(1,10,10) ydata = linspace(0,1,11) zdata = rand(length(ydata),length(xdata)) x0 = [2.5 3.5 4.5 5.5];...

12 meses hace | 0

Respondida
How to include NaN elements in an array based on a condition?
a = [2 4 6 12 14 18]; A = a(1):2:a(end); new_a = zeros(size(A)); idx = any(a'== A,1); new_a(idx) = a; new_a(~idx) = nan; n...

12 meses hace | 0

| aceptada

Respondida
Minimize the sum of squared errors between the experimental and predicted data in order to calculate two parameters
You can use MatLab's fmincon. z = [0.006;0.012;0.018;0.024;0.03;0.046;0.069;0.091;0.122;0.137;0.152;0.162]; u_exp = [0.709...

casi 1 año hace | 0

| aceptada

Respondida
Why do I receive "Not enough input arguments"? or How can I better code this ODE?
Here's a correct version of your code. 1) Time span and initial conditions must be in the same file as the ode solver, not sure...

alrededor de 1 año hace | 1

Respondida
I have solved a set of three coupled odes in MATLAB. I have to plot another function which is a function of one of these odes. How this can be done?
If none of the differential equations you are solving depends on deln, then you can compute it after solving the ODE system. To...

alrededor de 1 año hace | 1

| aceptada

Respondida
In this code, I want to find out slope at every value of x?
You could use the gradient function provided by MatLab. x=[10.33 10.331 10.332 10.333 10.334 10.335 10.336 10.337 10.338 10.339...

alrededor de 1 año hace | 0

Respondida
How to plot scatter plot with conditions?
data = readmatrix('table.xlsx'); time = data(:,1); PW = data(:,2); You can identify all the indexes for which the differenc...

alrededor de 1 año hace | 0

| aceptada

Respondida
How do I exclude values in linspace?
Hi @lea.c, I have several comments about your code. 1) MatLab already implements a function called integral, so I would recomme...

alrededor de 1 año hace | 0

Respondida
Loop over initial guesses with fsolve
Try with func=@f; x0=[1:0.1:1.9;1:0.1:1.9;1:0.1:1.9;1:0.1:1.9;0.01:0.01:0.1;0.01:0.01:0.1;0.01:0.01:0.1]; x = zeros(10,7); ...

alrededor de 1 año hace | 1

| aceptada

Respondida
Help with numerical integral
I guess this is what you're looking for e=0.1; f=@(x)exp(x)./(1+x); a=0.5; b=1.1; n=2; h=(b-a)/n; i=1:1:n-1; S=f(a...

alrededor de 1 año hace | 0

Respondida
Legend with certain number of rows and columns
Use legendflex from file exchange. x = linspace(0,1,1000); n = (1:6)'; y = x.^n; p = plot(x,y); title('y = x^n') legendf...

alrededor de 1 año hace | 1

| aceptada

Respondida
Creating a row vector of combinations?
You could use nchoosek. The example below is from -3 to 3 for memory issues (you'll generate a humongous matrix). Just substitut...

alrededor de 1 año hace | 0

| aceptada

Respondida
How to remove the zeros from the array upon breaking out of the for loop?
Would this suit your purpose? RHo = 0.5; % Relative Humidity air [%] M= 0.018; % molar mass of water mol...

alrededor de 1 año hace | 0

Respondida
How o define two different data set of (t , T) for an ode solver that it's variables are (t , C)?
I think your code is fine, below you can see a slightly different approach that uses interpolation to define the T(t) dependence...

alrededor de 1 año hace | 0

| aceptada

Respondida
Create a row vector comparing strings in a for loop
index = []; s1 = { 'FMO1','FMO2','FMO3','FMO4', 'BAO1', 'BAO2', 'BAO3', 'BAO4'}; s2 = {'mean','FMO1','BAO4','FMO2'}; for k =...

alrededor de 1 año hace | 0

| aceptada

Cargar más