Fabio Freschi

Politecnico di Torino

180 total contributions since 2014

Fabio Freschi's Badges

Answered
n! permutation matrices
Is this what you wish? % identity matrix A = eye(4); % permuatations idx = perms(1:4); % all matrices in a cell array B = ...

4 jours ago | 0

Answered
Solve the system of the 1-st order differential equations
Have a look at the ode suite

4 jours ago | 0

Answered
How to make two vectors the same size?
Try fkld=sum(dot(p,log(p./q))); with the element-wise division ./

4 jours ago | 1

| accepted

Answered
Inner matrix dimensions must agree
fkj= sum(.5*((p.^2-q.^2).^2)./((p.*q).^1.5)); Note that .^ .* and ./ operators are used for element-wise operation (see also St...

4 jours ago | 1

| accepted

Answered
Inverse matrix in Matlab very different to Excel
The answer is not very different. As Jeremy commented, the results are similar a part from numerical tolerances. You can verify ...

4 jours ago | 1

| accepted

Answered
How to call a multi output function without repeating the run
function [c, ceq] = nonconst(x) [total_value, Delta_value, energy_value] = digital(x); c = total_value-u*v;... ` `%fi...

5 jours ago | 0

| accepted

Answered
How to color a mesh by temperature(or other value) preserving the elevation. Using two matrix(same length)
% dummy data [X,Y] = meshgrid(1:0.5:10,1:20); Z = sin(X).*cos(Y); T = X+Y; % plot surf(X,Y,Z,T) colormap('hot') colorba...

10 jours ago | 1

| accepted

Answered
How can I get an array of velocity if I have an array of position and another array of time
First, you need a element-wise division pos = [0 2 4 8 16 32 64]; t = [0 0.1 0.2 0.3 0.4 0.5 0.6]; V = diff(pos)./diff(t); ...

10 jours ago | 0

| accepted

Answered
Numerical Solution for a System of THREE Differential Equations
This is an ODE, not PDE. Look for the ODE suite

11 jours ago | 0

| accepted

Answered
Area under curves with different length (trapz)
If you have the x axes for both waves, it is pretty easy to use trapz: N1 = 1036; N2 = 9958; % wave 1 x1 = linspace(0,2*pi...

11 jours ago | 1

| accepted

Answered
Plotting the integral of a graph
Suppose that x and y are the vectors defining your waveform. % dummy values x = linspace(0,2*pi,1000); y = sin(x); % integ...

11 jours ago | 1

| accepted

Answered
Does anybody know how to write a function which returns the coordinates of a convex hull of a list of certain points?
% your data N = 20; List = rand(N,2); % convex hull (k is the pointer to the rows in List with the convex hull k = convhul...

11 jours ago | 0

Answered
How to divide a part of array to the another part of it?
T = psi2(1:5)./psi(end-4:end);

11 jours ago | 0

| accepted

Answered
How can speed up the blow codes?
Edit: added attempt for ncount I don't understand what ncount is doing, it looks like a counter but I don't understand why you ...

12 jours ago | 0

Answered
defined n with the equation
n= size(P4,1)

27 jours ago | 0

Answered
Array indices must be positive integers or logical values
Your index k starts from 0 and you are indexing h with k, so you are asking at the first iteration h(0). Matlab has 1-indexing s...

27 jours ago | 0

Answered
while loop performance improvement
I don't know the starting value of k, assuming it is 1 idx = 1:rate:period mysum = sum(array(idx,2)); % note that sum is the n...

28 jours ago | 1

| accepted

Answered
Euler method error in code Index exceeds the number of array elements (1). in my my code
Checking the sintax of the code only I see 3 issues 1) In the for loop you wrote a(400-x(n)), that is you are addressing the po...

28 jours ago | 1

| accepted

Answered
Population growth in 7 years
% rate r = 0.0063; % years t = 7; % initial population a = 5e6; % population y = a*(1+r)^t

30 jours ago | 2

| accepted

Answered
Apply operation to each element of an array
You should use element-wise operations https://it.mathworks.com/help/matlab/matlab_prog/array-vs-matrix-operations.html In you...

30 jours ago | 0

| accepted

Answered
How do you solve for a double sum
There are some things that are vgue, like the x-axis, in any case, this code shows something similar to the desired result clea...

environ un mois ago | 1

| accepted

Answered
Solving System of Equations
It seems you are solving a numeric system of equations. Why don't you simply put the coefficient matrix in A and the right-hand-...

environ un mois ago | 0

Answered
Finding values within a matrix
It is difficult to answer because we don't have the data files 'areatemp1' and 'machtemp1'. In any case in your if statement yo...

environ un mois ago | 0

Answered
Making a 3D Matrix and Adding Matrices to make a bigger matrix
In the loop, you wrote T(i+k-1,j+k-1) = d(i,j,k-1); k starts from one, so you are trying to access the d(i,j,0) elemen...

environ un mois ago | 0

Answered
Computational Electromagnetics and MATLAB
This could be a good start https://it.mathworks.com/academia/books/computational-electromagnetics-with-matlab-sadiku.html

environ un mois ago | 1

| accepted

Answered
Why am i getting the : errir using inv ,Matrix must be square?
As James Tursa said, inv(H) is not defined when H is not square. If you have the matrix H and the rhs g, simply use f = H\g;

environ un mois ago | 0

Answered
Struggling to Vectorize the following Code.
Almost all execution time is due to the plot function as highlighted by the profiler: You can call the plot only once at the ...

environ un mois ago | 0

Answered
How can I pick certain element from each row of an array ?
I see that the edges shared by two triangles are not duplicated. To check for this case I sorted the nodes and the edges. Let me...

environ un mois ago | 0

| accepted

Answered
Combining my newton's method function for nonlinear equations and my lu decomposition function.
Is this what you are asking for? delta = LUgauss(J(x0),-F_value); Note that the buil-in lu Matlab function can do the work for...

environ un mois ago | 0

Answered
How to sum over rows within specific ranges without for loop?
Not sure if it is the best way, but this seems to work index = [0; cumsum(B)]; C = cell2mat(arrayfun(@(i)sum(A(index(i)+1:in...

environ un mois ago | 0

| accepted