回答済み
Much slower valid convolution using complementary size of kernels.
I would suggest to do specific conv with MEX programing. Not sure the chance to beat MATLAB though.

5年以上前 | 0

回答済み
Is there any method to accelerate many small matrix and vector's "mldivide" (4*4)?
I have submitted fast method FEX file MultipleQRSolver (C-compiler required) that can carried out the job in 1.5 second, 15 time...

5年以上前 | 0

回答済み
How to create an N-ary array
N=3 K=5 [~,A] = ismember(dec2base(0:N^K-1,N),['0':'9' 'A':'Z']); A = A-1

5年以上前 | 0

回答済み
How to a='10101111' convert array b=[1,0,1,0,1,1,1,1]
>> a='10101111' a = '10101111' >> a-'0' ans = 1 0 1 0 1 1 1 1

5年以上前 | 3

回答済み
How to average along all permutations of a 4D array?
A=randn(10,10,10,10); p=num2cell(perms(1:4),2) B=cellfun(@(p) permute(A,p), p, 'unif', 0); B=mean(cat(5,B{:}),5);

5年以上前 | 1

回答済み
set order of elseif
prefered_order = 'xzy'; b = [x y z]; % logical conditions (in if/elseif) code_to_be_exec = { @xcode, @ycode, @zcode }; % cod...

5年以上前 | 2

回答済み
Creating a non-square diagonal matrix
B = diag(A(:)); B(94037,1) = 0;

5年以上前 | 1

回答済み
Generating matrix with ones and zeros
Look at KRON >> kron(eye(4),ones(1,3)) ans = 1 1 1 0 0 0 0 0 0 0 0 0 ...

5年以上前 | 0

回答済み
Quadratic Spline Interpolation Code
My toolbox provide any order (1D) spline, including quadratic https://www.mathworks.com/matlabcentral/fileexchange/25872-free-k...

5年以上前 | 0

回答済み
is this repr.m redundant to existing standard approach?
A decend tool https://www.mathworks.com/matlabcentral/fileexchange/29457-serialize-deserialize A true matlab serialization is...

5年以上前 | 0

回答済み
3 equations 4 unknowns
>> M=[0 1 6; 4 0 4; 0 1 0]; >> [X,K]=eig(M); >> k=1./diag(K); % select 3rd eigen vector for example >> X(:,3) ans = ...

5年以上前 | 0

回答済み
Whats the difference between the two statements
The second creates 2-row matrix. The first creates 2-column matrix, since it make a transpose after reshape.

5年以上前 | 0

回答済み
Alternative to blkdiag and mat2cell functions
[I,J] = ndgrid(1:Nb,1:Nc); C = accumarray([I(:),J(:)+(I(:)-1)*Nc],A(:));

5年以上前 | 0

回答済み
Alternative to blkdiag and mat2cell functions
What about about a simple for-loop [Nb,Nc] = size(A); C = zeros(Nb,Nb*Nc); for r=1:Nb C(r,(r-1)*Nc+(1:Nc)) = A(r,:); en...

5年以上前 | 0

回答済み
Inexplicable different execution speeds when filling a matrix with complex entries
I also dig deeper into this "issue" of slowness. The strange issue I discover is that when filling with real values into a comp...

5年以上前 | 0

| 採用済み

送信済み

Set partition
List all partitions a set n elements

5年以上前 | ダウンロード 1 件 |

4.8 / 5

回答済み
Inexplicable different execution speeds when filling a matrix with complex entries
You should preallocate M as complex M = 1i+zeros(m,n); if you intend to populate it with complex number later in the for-loop....

5年以上前 | 0

回答済み
Sum of selected elements in Matrix
A=rand(1000); [m,n]=size(A); G=max((1:m)',1:n); s=cumsum(accumarray(G(:),A(:)))

5年以上前 | 0

回答済み
Trying to simplify code to see if 3 vectors are at right angle. Unsure of outcome.
>> dot(Q-P,R-P) ans = 0 So PQ is right angle with PR, etc...

5年以上前 | 0

回答済み
Excluding one vector from another vector with repetition
From Alexander Gallard above comment (I agree he should open a new thread) "This method didn't work for me if my C was, say, [...

5年以上前 | 0

回答済み
How can I create a set from elements and combined elements?
Not sure exactly what kind of "combination" you look for, but it looks like a subset of these b=logical(dec2bin(1:2^4-1,4)-'0')...

5年以上前 | 1

回答済み
matrix manupulation and reading diagonally
Matrix_test = [13,10,5,2;7,14,11,6;3,8,15,12;1,4,9,16] [m,n]=size(Matrix_test); [i,j]=ndgrid(1:m,1:n); jmi = j(:)-i(:); [~...

5年以上前 | 0

回答済み
Point or multiple points is/are in a triangle??
tf = inpolygon(Points(:,1),Points(:,2),Triangle(:,1),Triangle(:,2))

5年以上前 | 0

| 採用済み

回答済み
Does using Functions make your code run faster?
MATLAB makes optimization when the code is in function, it does not in script (why it's also mysterious to me, there might be th...

5年以上前 | 1

回答済み
How to create N+1 dimensional array by taking exterior product of 1st dimension of two N dimensional arrays?
Solution without expansion as requested sizeA = size(A); reshapeB = reshape(B,[1 size(B)]); reshapeA = reshape(A,[sizeA(1) 1 ...

5年以上前 | 0

| 採用済み

回答済み
When is minimum p-norm solution independent of p?
The optimizers fail to initialize find the first feasible point with its gradient and return the same guess, which are identical...

5年以上前 | 0

回答済み
Curve fit or spline fit for a wavy function
You can use my tool https://www.mathworks.com/matlabcentral/fileexchange/25872-free-knot-spline-approximation load('RPM.mat') ...

5年以上前 | 0

回答済み
How to check 2D rectangle and 3D rectangular prism intersect with each other?
Quick and dirty code load('prism.mat') load('rectangle.mat') % Normalize the coordinates M = prism; minM = min(M,[],1); ...

5年以上前 | 0

回答済み
how to write special matrices
>> dec2bin(0:2^4-1)-'0' ans = 0 0 0 0 0 0 0 1 0 0 1 0 0 0 ...

5年以上前 | 2

回答済み
Finding shortest path on calculated trajectories via A*
Check out this of you nsist on using A* https://www.mathworks.com/matlabcentral/fileexchange/10922-matlabbgl Otherwise Matlab ...

5年以上前 | 0