Bo

Dear Harry, Thank you for your help. The first question is answered by setting PTMX. For the second question, I just rewrote the function coco_add_func. Best Regards, Bo

Dear Harry, Thank you for your help. The first question is answered by setting PTMX. For the senod question, I just rewrote the function coco_add_func. Best Regards, Bo

Dear Harry, Thank you for your help. The first question is answered by setting PTMX. For the senod question, I wil try again. Best Regards, Bo

Dear Harry, Thank you for your reply. I am not sure I understand what you mean. why can we not calculate the entire frequency response curve at conce? Reason 1: If I take the excitation period T=2pi/5 as the initial parameter, so that in the jump region of the frequency response curve, the actual excitation period T may be not continue with the frequency response curve (the excitation period has changed, not continued). Instead, we need to constrain the curve from the other point, such as T=2pi/20....

Dear Harry, Thank you for your reply. I have changed the codes based on your comments. I have a few questions I would like to ask for your help. In the line 23 of demo_po.m, cont_args = {1, {'po.period' 'T'}, [2*pi/200 2*pi/2]}; it appears that the period of the system is limited between 2 pi/200 and 2 pi/2 (gen_sym_po.m and demo_po.m) . The result is shown in Figure 1. If I change the following to read p0=[2*pi/20; 1]; [~,x0]=ode45(@(t,x)funcs{1}(t,x,p0),[0 20*pi],[0.01; 0; 0; 0; 0; 0; 0; 0; 0;...

Dear Harry, Thank you for your reply. I have modified the codes based on your comments. I have a few questions to ask you for help. In the line 23 of demo_po.m, cont_args = {1, {'po.period' 'T'}, [2*pi/200 2*pi/2]}; it appears that the period of the system is limited between 2 pi/200 and 2 pi/2 (gen_sym_po.m and demo_po.m) . The result is shown in Figure 1. If I change the following to read p0=[2*pi/20; 1]; [~,x0]=ode45(@(t,x)funcs{1}(t,x,p0),[0 20*pi],[0.01; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0]); [t0,x0]=ode45(@(t,x)funcs{1}(t,x,p0),[0...

Dear Harry, Thank you for your reply. I have modified the codes based on your comments. I have a few questions to ask you for help. In the line 23 of demo_po.m, cont_args = {1, {'po.period' 'T'}, [2*pi/200 2*pi/2]}; it appears that the period of the system is limited between 2 pi/200 and 2 pi/2 (gen_sym_po.m and demo_po.m) . The result is shown in Figure 1. If I change the following to read p0=[2*pi/20; 1]; [~,x0]=ode45(@(t,x)funcs{1}(t,x,p0),[0 20*pi],[0.01; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0]); [t0,x0]=ode45(@(t,x)funcs{1}(t,x,p0),[0...

Dear Harry, I have a problem plotting the frequency amplitude curve. I am following the example for plotting the frequency amplitude curve of bistable, the following codes are constructed. Unfortunately, MX is displayed. I do not see in the documentation how to adjust the continuation settings to avoid MX. The other equation is about the frequency amplitude curve in relation to the excitation frequency Omega https://sourceforge.net/p/cocotools/tickets/51/#021b. Are there details on how to plot this?...