Community Profile

# Fabio Freschi

### Politecnico di Torino

Last seen: 25 días ago Active since 2014

All

#### Content Feed

View by

Minimal-code to get axes limits of log-log plots, but using gscatter
You can remove the negative part of the axis before setting the log scale set(gca,'XLim',[0 max(xlim)],'YLim',[0 max(ylim)]); ...

2 meses ago | 0

How to plot Temperature vs Efficiency curve using For loop and Linspace?
I have commented your plots to avoid moltiple figures % numebr of temperatures nT = 100; % temperature vector T = linspace(2...

2 meses ago | 1

| accepted

Question

Calculation of consistent initial conditions for ODE
(minor edit) Hi everyone, I have written wrapper functions calling MATLAB's ODE solvers. Specifically, my problems are DAEs wi...

2 meses ago | 1 answer | 0

### 1

How to find same unique in two different size arrays?
Use ismember % data x=[1 22 22 3 8 19]; y=[22 10 19 35]; % engine [iii,jjj] = ismember(y(:),x(:)); countCommon = su...

2 meses ago | 0

| accepted

Having trouble using for loops to create a matrix
I reproduce the result of the picture, because the question is not clear. a = 25 ; b = 5 ; c = (2*a)/b - 1 ; % do not use mi...

3 meses ago | 0

| accepted

How to plot a function
You may use fplot figure, hold on fplot(@(x)x.^-3,[1 1e5]) fplot(@(x)x.^-2.8,[1 1e5]) % use bi-logaritmic scale set(gca,'...

3 meses ago | 0

Ax = 0 with 2 values of x known
The post is old but the question arises many times, so I post here my solution. It is a generalization of @Matt J solution. Sup...

3 meses ago | 0

Looking for visualization method to present finite difference numerical solution in 2D
You may use the low-level patch function. Since I don't have your data, please find below a self-consistent example showing the ...

5 meses ago | 0

smooth step function with rise time
My colleague Sofia have found a paper [1] where authors report the analytic approximation of the piece-wise linear step function...

10 meses ago | 1

| accepted

Question

smooth step function with rise time
Hi everyone, I know this question is more related to math than Matlab, but I know there are many gurus out there who can easily...

10 meses ago | 3 answers | 1

### 3

How can I calculate the sum of the absolute error?
Please check the code below. clear all, close all % model F = @(c,x)c(1)+c(2)*x+c(3)*(x.^2-1); % data xk = [-2 -1 0 1 ...

11 meses ago | 1

| accepted

Drawing a rotating circle on a sphere
The code should be self-explicative [X,Y,Z] = sphere; % figure figure, hold on axis equal, view([1 1 1]); % sphere wit...

11 meses ago | 0

| accepted

Showing me wrong Matrix Size.
multiply without the 'dot' OO= N*kk;

11 meses ago | 0

How can you solve a system of second order differential equation, which is not time, but length and radius dependent?
This is not an ODE but a PDE. you should have a look at pdepe

11 meses ago | 2

| accepted

We have a column vector and a row vector both of same lengths. How to use sum function to get element wise sum of two vectors and final answer display in row vector
vect3 = vect1.'+vect2 If you want to be sure of the element by element sum, regardles the shape of the two vectors (row or colu...

11 meses ago | 0

| accepted

How to randomly change element of matrix?
You can use linear indexing % your matrix A = rand(10,100); % create the folder mkdir('storage'); % loop for i = 1:100...

11 meses ago | 2

How to create a series of circle by using a specific radius in cartesian?
Try something like this % anonymous function to plot a circle given the center coordinates and the % radius plotcircle = @(x,...

11 meses ago | 0

| accepted

Converting Euler ODE to the MATLAB code
f is an anonymous function f = @(t,y)[-2*y(1)+y(2)+15*cos(t); 2*y(1)-y(2)]; note that your script is not correct, since t is t...

11 meses ago | 1

| accepted

How to call a variable matrix from one script to another script
I am not sure I have understtod correctly your problem. I give it a try. I would add an extra layer using an anonymous function...

11 meses ago | 1

| accepted

How to use ~= operator as constraint function?
You may have a look here or here

11 meses ago | 0

How to compare two time signals with different amount of samples?
Use interp1 to resample your vector(s) in order to have the same length clear variables, close all % time vector t1 = 0:1/1...

11 meses ago | 1

solve efficiently sparse large linear system of equations
I would let the backslash operator work on this matrix sol = L\f; Can you share the original matrix to make other attempts?

11 meses ago | 0

how to use lsqnonlin to solve conditional equation
The following code should be self explainatory. In the opposite case, simply ask clear all, close all % some dummy params va...

11 meses ago | 0

| accepted

how to ascribe two values to a parameter such that it will be time dependent
You can pass an extra parameter to RiceG1 function using an anonmymous function. In your driver function add % extra anonimous ...

11 meses ago | 1

Can someone help me implement my Kuramoto Model with time delay please?
Your function requires additional parameters. One way to pass them is to use an anonymous function as extra layer % params N =...

12 meses ago | 0

| accepted

Euler's method in matlab
Any ODE solver gives the unknown function(s) as result, not its derivative. So, yes, your results are x and y. Note that in you...

12 meses ago | 0

| accepted

for loop/periodic implementation for number or arrays?
If you want to force the max and min values of your vectors, you can preprocess them x0(x0 < 1 & x0 > nx+1) = 1; If the bounda...

12 meses ago | 1

| accepted

How to get this while cycle?
if B and C are matrices, make them vectors using : while any(abs(a_0-a) > 1e-6) || any(abs(B_0(:)-B(:)) > 1e-6) || any(abs(C_0(...

12 meses ago | 1

| accepted