MATLAB: Can I multiply a Vector with a Matrix

arrayMATLABmatrixmatrix arrayvector

Hi;

Sorry for my english. I've a big problem, I want multiply a vector by a matrix:

A part of my code is:

[...]
x = -5000:10.:5000;
y = -5000:10.:5000;
[X,Y]=meshgrid(x);
U=   []; %is an array <1x15000 double>

teta=[]; %is an array <1x15000 double>
%condizioni iniziali
if U <= 5500
    rho = 0.1;
else U> 5500;
    rho = 0.16;
end
u = ((DeltaH/Hsor)^rho).*U ; 
%velocità del vento corretta
h = Hsor+DeltaH-((vg*2.1./u));
%altezza considerando la sedimentazione in atmosfera
dd=((270-teta)*pi)/180;
%Modifica direzione vento
x11 = cos(dd);
x111 = x11*X;
[...]

I've a problem in the last operation ??? Error using ==> mtimes Inner matrix dimensions must agree.

 Error in ==> SCORMAT at 90
 x111 = x11*X;

[EDITED, Jan, code formatted]

Best Answer

Change x11*X to x11.*X

Related Solutions

MATLAB: How to index matrices inside a loop

Hello,

You can change your code a bit like this

%% Quantum Control
format long
clear all
clc
%% Variables
time = -20000:.1:20000; % Time

n = length(time);
U = zeros(1,n);
        for i = 1:n % Time
        ft = time(i);    
       
        alpha = 0.001; % Transition parameter

        
        beta = 0.001; % Transition parameter
        
        ai = 0.2; % Initial population |1>
        
        ai2 = 0.8; % Initial population |2>
        
        af = 1; % Final population |1>
        
        w0 = 0.02; % w0 = (E2 - E1)/h
        
        mi = 6; % dipole operator
        
        phii = pi/24; % initial relative phase
        
        phif = pi/4; % final relative phase
        
        E1 = 1.98; % eigen-energy 1
        E2 = 2; % eigen-energy 2
        H = [0 1;1 0]; % Interation hamiltonian
        H0 = [E1 0;0 E2]; % hamiltonian
%% Control function
    
        g = 1./(1+exp(-alpha.*ft));
           
        f = ai*(1-g)+af*g;
         
        p = 1./(1+exp(-beta.*ft));
        
        h = phii*(1-p)+phif*p;
        
%% Electric field 
        
        E = alpha.*(af-ai).*exp(alpha.*ft).*sin(w0.*ft+h)./(mi.*(1+exp(alpha.*ft)).*sqrt((1-ai+(1-af).*exp(alpha.*ft)).*(ai+af.*exp(alpha.*ft))))+2.*beta.*(phif-phii).*exp(beta.*ft).*sqrt(f.*(1-f)).*cos(w0.*ft+h)/(mi.*(1-2.*f).*(1+exp(beta.*ft).^2));
    
%% Time evolution operator
        [V, D]=eig(H);
        
        Z = expm(H0);
        Z1 = expm(D);
        
        U(i) = exp(-1i*ft/2)*Z*exp(1i*mi*E*ft)*V*Z1*inv(V)*exp(-1i*ft/2)*Z
        
        end
    
%% End

Now, U will be a vector.

MATLAB: Find(vector==value) not working

You have encountered floating-point approximation error. See the documentation on Floating-Point Numbers for a full explanation.

Try this:

vector=0:0.01:0.3;
Result = find(vector<=0.16, 1, 'last');

producing:

Result =
    17

Best Answer

Related Solutions

MATLAB: How to index matrices inside a loop

MATLAB: Find(vector==value) not working

Related Question