Community Profile

photo

Davide Masiello


Paul Scherrer Institut (PSI)

Last seen: Today Attivo dal 2017

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

Statistics

All
  • Knowledgeable Level 5
  • Pro
  • Introduction to MATLAB Master
  • 3 Month Streak
  • Thankful Level 2
  • Knowledgeable Level 3
  • Solver
  • First Answer
  • Thankful Level 1

Visualizza badge

Content Feed

Visto da

Risposto
Solving a differential equation
Since you have data for V and C you could fit them with the analytical solution of the differential equation. Let's assume this...

circa 2 ore fa | 0

Risposto
solving a dynamic equation
Hi @sam - you did not provide enough info for us to try and give a complete answer, but I did some guessing and came up with a s...

circa 18 ore fa | 0

Risposto
Forward Euler solution plotting for dy/dt=y^2-y^3
There are a couple of issues with your code 1) Indexes in MatLab start at 1, not 0, so y(0) is not valid syntax and must be rep...

13 giorni fa | 0

| accettato

Risposto
Dot indexing is not supported for variables of this type.
You are assuming that sh_16_data is a structure. Instead, do sh_16_freq = sh_16_data(:,2)

18 giorni fa | 0

Risposto
find columns adjacent to logical values
Reductive example. Assume this is your matrix A = rand(10,5) and let's assume your logical criteria is that you want to extra...

19 giorni fa | 0

| accettato

Risposto
Loop plotting for individual years
There's no legend in the example below, because I am not sure what you wanted to do with the cellfun function. Nevertheless, it...

20 giorni fa | 0

| accettato

Risposto
How to vertically concatenate many tables produced from a for loop
I would try something like this M = []; for k=1:numel(A) % automatically brings in the correct amount of tables M = [M;re...

20 giorni fa | 0

| accettato

Risposto
Duplicating an entire vector
A=[1 2 3 4]; repmat(A,1,3)

20 giorni fa | 1

Risposto
plotting graph with 2 different size of vector
Applying diff to a Nx1 array produces an (N-1)x1 array, hence your problem. In truth, you are not really calculating the deriva...

20 giorni fa | 0

| accettato

Risposto
Build a matrix with a gradual value of intensity and an horizontal line with a constant value
A = repmat(linspace(0,1,1000),500,1); A(220:280,:) = 0.5; imshow(A)

20 giorni fa | 0

| accettato

Risposto
removing specified data from variable
I think the following code is a simpler way of achieving your task, but it does not implement the "pulling a number from a rando...

21 giorni fa | 0

Risposto
Call different functions inside a function
The way you have explained it in your question is already correct. See the example below. %% Main script mainfunction %% Mai...

21 giorni fa | 1

Risposto
Create an airfoil in MATLAB
From the wikipedia page https://en.wikipedia.org/wiki/NACA_airfoil t = 0.25; % maximum thickness as a fraction of the ch...

24 giorni fa | 0

Risposto
Saving Multiple .dat files
This should work A = rand(2000,1000); for col = 1:size(A,2) filename = ['Column_',num2str(col),'.dat']; writematri...

24 giorni fa | 0

| accettato

Risposto
MATLAB double addition format long
I suggest reading the following thread https://ch.mathworks.com/matlabcentral/answers/1838113-wrong-sum-50-0-01

25 giorni fa | 0

| accettato

Risposto
How do I find slope for large dataset?
Just do A = diff(data)./diff(time);

25 giorni fa | 0

| accettato

Risposto
How do i set up a legend for a graph with variable parameters?
t = -2:0.1:4; d = 1; c = 1; A = 1:0.2:2; figure hold on for i = 1:length(A) y= A(i) * (t - d).*(t-c) ; plot (...

25 giorni fa | 0

Risposto
Generate Alternating Tones Based On Number Of Cycles
I am still quite unsure if this is what you are asking, but I will give it a go. [t,w] = wave(465); plot(t,w) function [t,f] ...

26 giorni fa | 0

Risposto
Interpolate data to present with limited size of data
x=rand(5,7) figure contourf(x,'LineColor','none') shading interp axis equal off

26 giorni fa | 1

| accettato

Risposto
Getting A blank Graph
See below for a correct use of indexing with a for loop x = 0:0.01:1; for i = 1:length(x) y1(i) = (0.5*x(i)^5)-(0.6*x(i)^4)...

26 giorni fa | 0

Risposto
trying to plot 3 variables in 3d
surf can handle only matrices, so this is what you want to do theta= linspace(.01, .55); %for thetamax = 32 degrees M=linspace...

26 giorni fa | 0

| accettato

Risposto
While loop with all elements meeting the conditions
t = -100:200; F = t; indices = 1:100; while all(F(indices)<0) indices = indices + 1; end

circa un mese fa | 1

| accettato

Risposto
Wrong sum 50 + 0.01
I don't think this is an issue at all, but it rather has to do with the max accuracy MatLab uses to store numbers. Such accurac...

circa un mese fa | 0

| accettato

Risposto
solve multiple equation in one variable
I think the problem is that you have 5 equations and 4 variables, because p1 does not appear anywhere in the system. Therefore ...

circa un mese fa | 0

| accettato

Risposto
How can i get a logical vector for all rows when compare each row of a column matrix with a another row column matrix
You can simply do A = randi(10,5) B = randi(10,5) l = A == B

circa un mese fa | 0

| accettato

Risposto
How to delete previous values within a certain difference?
You can do this A = [1,2,3,4,20,21,22,35,36,37,38,50,51,52,53,54,57,70]; A = A(diff(A)>10)

circa un mese fa | 0

| accettato

Risposto
How do I set some variables of a function to a constant and plot the others?
This is a possibility E = -3:.001:3; plot(E,Theta(E,1)) function out = Theta(E,Delta) out = zeros(size(E)); out(abs(E)<=D...

circa un mese fa | 1

| accettato

Risposto
How to delete adjacent values
A = [1 5 6 7 8 9 10 11 12 20 21 22 23]; A([false,diff(A) == 1]) = []

circa un mese fa | 0

| accettato

Risposto
IF statement not generating results?
data=load('artificial_slope_2.txt') %loading data % defining input values based of input txt file nx = data(:,4); ny = data(...

circa un mese fa | 0

Risposto
a problem in rectangular function.
You are using plot incorrectly, see below. a=1; syms t f=piecewise(a / 0.2e1 < abs(t), 0, abs(t) == a / 0.2e1, 0.1e1 / 0.2e1,...

circa un mese fa | 1

| accettato

Carica altro