How to find the difference between two vectors

Question

vani 2022-10-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1824053-how-to-find-the-difference-between-two-vectors

评论： dpb 2022-10-14

采纳的回答： dpb

在 MATLAB Online 中打开

Hello all,

My doubt is,

I have to find the difference between two vectors( for eg: y(8) and y(4) or y(6) and y(2)) and have to append the difference output also get the appended output in plot. So could anyone please tell me how to do these functions?.

y_try=[0 1e-9 0 0 0 1e-9 0 0] in this
y(1) =[0] ; y(2) = [1e-9] ; y(3) =[0] ; y(4) =[0] 
y(5) =[0]; y(6) = [1e-9]; y(7) = [0] ; y(8) = [0]

3 个评论
显示 1更早的评论隐藏 1更早的评论

vani 2022-10-12

I have to find the difference of output for each vectors (y(8) and y(4) , y(6) and y(2), y(7) and y(3), y(5) and y(1)) and append all these difference output and plot as a figure.

vani 2022-10-13

Thank you for your response. The difference is working and now I have problem with appending.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

dpb 2022-10-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1824053-how-to-find-the-difference-between-two-vectors#answer_1072878

在 MATLAB Online 中打开

Well, this is a guess, but it appears that maybe you have x,y data strung together in a single vector and want to separate out the two and difference -- try

y_try=[0 1e-9 0 0 0 1e-9 0 0];

d=diff(reshape(y_try,[],2),[],2);

plot(d)

Of course, with the above dataset the differences are all identically zero; this may be just owing to the particular values used as example or, it could be a complete misunderstanding of the intent of the question...

10 个评论
显示 8更早的评论隐藏 8更早的评论

vani 2022-10-13

在 MATLAB Online 中打开

function cantilever_trialadd
global count
global max
global env
global damping1
global damping2
global omega_0_square1
global omega_0_square2
global freq_sound
global F_ext
global alpha
global beta
global gamma
global a_f1
global a_f0
global R
global u_dc1
global u_dc2
global kalib1
global kalib2
global freq_act
global u_ac
global tau
global a_flist
global T_list
global difference_list
global difference_aflist
global flag
global t0
global sampling
    R=15;
    freq_res1= 1000;
    freq_res2=1050;
    bandwidth1 = 100;
    Q1 = freq_res1/bandwidth1;
    bandwidth2 = 150;
    Q2 = freq_res2/bandwidth2;
    omega_0_1=freq_res1*2*pi;
    omega_0_2=freq_res2*2*pi;
    omega_0_square1=omega_0_1^2;
    omega_0_square2=omega_0_2^2;
    damp1=1/(2*Q1);
    damp2=1/(2*Q2);
    damping1=damp1*2*omega_0_1;
    damping2=damp2*2*omega_0_2;
    alpha=19.2;
    beta=1.0066e+03;
    gamma=1.6214e+07;
    kalib1=1e6; 
    kalib2=2e6; 
    u_dc1=-0.2;
    u_dc2=-0.1;
    a_f0 = 0;
    tau=0.001;
    F_ext=0.5;
    freq_sound=1000;
    t_final = 0.2;
    zeit = [0 t_final];
     y_start2=[0 1e-9 0 0 0 1e-9 0 0];
    options = odeset('MaxStep',1e-6);
       a_f1 = 0.3;
    [T1,Y1] = ode15s(@(zeit,y_start2) cantimodel_temp(zeit,y_start2),zeit,y_start2,options);
    figure;
    plot(T1,Y1(:,2),'r--');
    hold on;
    plot(T1,Y1(:,6),'b--');
    legend('canti1','canti2');
    hold off;
end
function dy = cantimodel_temp(t,y)
    global damping1
    global damping2
    global omega_0_square1
    global omega_0_square2
    global freq_sound
    global F_ext
    global alpha
    global beta
    global gamma
    global tau
    global a_f0
    global a_f1
    global R
    global u_dc1
    global u_dc2
    global kalib1
    global kalib2
    global env
    global max
    global count
    global t0
        dy(1) = -beta.*y(1)+gamma.*(tanh(a_f0.*kalib1.*y(4)+u_dc1)/R)^2;
        dy(2) = y(3);
        dy(3) = -damping1*y(3)-omega_0_square1*y(2)+alpha*y(1)+F_ext*sin(2*pi*freq_sound*t);
        dy(4)=(-y(4)/tau+y(3));
        dy(5) = -beta.*y(5)+gamma.*(tanh(a_f1.*kalib2.*y(8)+u_dc2)/R)^2;
        dy(6) = y(7);
        dy(7) = -damping2*y(7)-omega_0_square2*y(6)+alpha*y(5)+F_ext*sin(2*pi*freq_sound*t);
        dy(8)=(-y(8)/tau+y(7));        
      dy=dy';
      diff=(y(8)-y(4));
      d=diff(reshape(y_start2,[y(8)],2),[y(4)],2);  %(I have inserted a code here) 
    t(end);
end

dpb 2022-10-13

编辑：dpb 2022-10-13

在 MATLAB Online 中打开

Well, in

d=diff(reshape(y_start2,[y(8)],2),[y(4)],2); %(I have inserted a code here)

you didn't even bother to switch the variable name to the variable you used in your code instead of that of the original example.

But even before you got there, in

diff=(y(8)-y(4));

you aliased the MATLAB builtin diff() function with a local variable so broke the use of diff as a function. Get rid of that line entirely and don't use builtin function names as variables; "There be dragons!"

...
dy(1) = -beta.*y(1)+gamma.*(tanh(a_f0.*kalib1.*y(4)+u_dc1)/R)^2;
dy(2) = y(3);
dy(3) = -damping1*y(3)-omega_0_square1*y(2)+alpha*y(1)+F_ext*sin(2*pi*freq_sound*t);
dy(4)=(-y(4)/tau+y(3));
dy(5) = -beta.*y(5)+gamma.*(tanh(a_f1.*kalib2.*y(8)+u_dc2)/R)^2;
dy(6) = y(7);
dy(7) = -damping2*y(7)-omega_0_square2*y(6)+alpha*y(5)+F_ext*sin(2*pi*freq_sound*t);
dy(8)=(-y(8)/tau+y(7));        
dy=dy(:);                       % make it a column vector
d=diff(reshape(dy,[],2),[],2);  % reshape to 2 columns and subtract the two

will return the difference vector you're asking for, but it's only going to be for a single point in time and am not at all sure what you think it means to have done so.

You'll have to figure out from your application how to return it to the caller and what to try to do with it; I'm at a loss to have a klew beyond simply computing what it is you've asked for, sorry.

vani 2022-10-14

Thank you so much for your response. Your code helped a lot to rectify the error.

dpb 2022-10-14

Glad to help where can...even though sometimes we don't have enough information to have any idea what the result means, can still at least fix syntax errors or find the correct command for the job.

If the above resolved the problem, please go ahead and "Accept" the answer if or no other reason than to let others know it's not still an outstanding issue.

请先登录，再进行评论。

How to find the difference between two vectors

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

10 个评论
显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

How to find the difference between two vectors

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

10 个评论 显示 8更早的评论隐藏 8更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

3 个评论
显示 1更早的评论隐藏 1更早的评论

10 个评论
显示 8更早的评论隐藏 8更早的评论