本文整理汇总了C++中dec2phase函数的典型用法代码示例。如果您正苦于以下问题：C++ dec2phase函数的具体用法？C++ dec2phase怎么用？C++ dec2phase使用的例子？那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了dec2phase函数的15个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: pulsesequence
//.........这里部分代码省略.........
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
obspower(tpwr); /* Set power for pulses */
dec2power(pwNlvl); /* Set decoupler2 power to pwNlvl */
decpower(pwClvl);
decpwrf(rf0);
obsoffset(tof);
decoffset(dofCHn);
delay(d1);
status(B);
rcvroff();
txphase(t1);
delay(9.0e-5);
/* spin-echo filter and 15N/13C double filters */
rgpulse(pw,t1,0.0,0.0);
if(C13filter[A]=='y')
{
decphase(t5);
zgradpulse(gzlvl5,gt5);
delay(tauCH1-gt5);
decrgpulse(pwC,t5,0.0,0.0);
txphase(t3); dec2phase(t5);
delay(tauNH1-tauCH1-pwN*0.5);
sim3pulse(2.0*pw,2*pwC,pwN,t3,zero,t5,0.0,0.0);
decphase(t6);
delay(tauCH2+tauCH1-tauNH1-pwN*0.5);
decrgpulse(pwC,t6,0.0,0.0);
delay(tauNH1+tauNH2-tauCH1-tauCH2-gt5-gstab);
zgradpulse(gzlvl5,gt5);
delay(gstab);
}
else
{
if (gzlvld1>0.0)
{
txphase(t3); dec2phase(t5);
delay(2.0e-6);
zgradpulse(gzlvld1,0.5*fact*tauNH1-pwN*0.25-11.0e-6);
delay(2.0e-5);
zgradpulse(-gzlvld1,0.5*fact*tauNH1-pwN*0.25-11.0e-6);
delay(2.0e-6);
sim3pulse(2.0*pw,0.0e-6,pwN,t3,zero,t5,0.0,0.0);
delay(2.0e-6);
zgradpulse(gzlvld1,0.5*fact*tauNH2-pwN*0.25-11.0e-6);
delay(2.0e-5);
zgradpulse(-gzlvld1,0.5*fact*tauNH2-pwN*0.25-11.0e-6);
delay(2.0e-6);
}示例2: pulsesequence
//.........这里部分代码省略.........
if( ix == 1) d2_init = d2 ;
t1_counter = (int)((d2-d2_init)*sw1 + 0.5);
if((t1_counter % 2))
{
tsadd(t1,2,4);
tsadd(t16,2,4);
}
if( ix == 1) d3_init = d3 ;
t2_counter = (int)((d3-d3_init)*sw2 + 0.5);
if((t2_counter % 2))
{
tsadd(t2,2,4);
tsadd(t16,2,4);
}
decstepsize(1.0);
initval(phi_CO, v1);
initval(phi_Ca, v2);
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
delay(d1-1.0e-3);
obsoffset(tof);
decoffset(dof);
obspower(tpwr);
decpower(pwClvl);
decpwrf(4095.0);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
rcvroff();
if(gt6 > 0.2e-6)
{
delay(10.0e-6);
decrgpulse(pwC, zero, 1.0e-6, 1.0e-6);
delay(0.2e-6);
zgradpulse(gzlvl6, gt6);
}
decpwrf(dfCa180);
delay(1.0e-3);
rgpulse(pw,zero,1.0e-6,1.0e-6);
delay(2.0e-6);
zgradpulse(gzlvl0,gt0);
delay(taua - gt0 - 2.0e-6 - WFG_START_DELAY);
simshaped_pulse("","offC17",2.0*pw,pwCa180,zero,zero,1.0e-6,1.0e-6);
/* c13 offset on CO, slp 180 on Ca */
delay(taua - gt0 - 500.0e-6 - WFG_STOP_DELAY);
zgradpulse(gzlvl0,gt0);
txphase(one);
delay(500.0e-6);
rgpulse(pw, one, 1.0e-6, 1.0e-6);
decphase(zero);
delay(2.0e-6);
zgradpulse(gzlvl3,gt3);
obspower(tpwrHd);
decpwrf(dfCa90);
delay(200.0e-6);示例3: pulsesequence
//.........这里部分代码省略.........
tau3 = 0.5*(d4_index/swC+0.5/swC)-pw -rof1 -pwC*2.0/M_PI ; /* increment corresponds to 13C increment */
tau3p = 0.5*(d4_index/swN+0.5/swN) -pw -rof1 -pwC -pwN*2.0/M_PI -tau3;
if(d4_index % 2) { tsadd(t7,2,4); tsadd(t8,2,4); tsadd(t31,2,4); }
tau2 = d3;
tau2 += 0.0*(-pw*4.0/M_PI-rof1*2.0);
if((f2180[A] == 'y') && (ni2 > 0.0)) {tau2 += ( 1.0 / (2.0*sw2) ); }
if(tau2 < 0.2e-6) {tau2 = 0.0;}
tau2 = tau2/2.0;
}
/* simultaneous Ntrosy-ChsqcSE, last part */
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (phase1 == 1) {icosel = 1; }
else { tsadd(t21,2,4); tsadd(t22,2,4); tsadd(t23,2,4); icosel = -1; }
if(d2_index % 2) { tsadd(t4,2,4); tsadd(t5,2,4); tsadd(t31,2,4); }
/* ECHO-ANTIECHO + STATES-TPPI t1, t1' in TROSY/HSQC last step */
tau1 = 1.0*d2_index/swC; /* increment corresponds to 13C increment */
tau1p = 1.0*d2_index*(1.0/swN-1.0/swC);
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
delay(d1);
/* Destroy 13C magnetization*/
decrgpulse(pwC*1.0, zero, 0.0, 0.0);
zgradpulse(-gzlvl0, gt0);
delay(gstab);
/* NOESY */
if(exp_mode[A]!='2')
{ /* 3-4D */
if(exp_mode[A]=='4')
{ /* full 4D */
/* t3 evolution, the very first HSQC */
/* Hz -> HzXz INEPT */
rgpulse(pw,two,rof1,rof1); /* 1H pulse excitation */
zgradpulse(gzlvl7, gt0); delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0); delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab); zgradpulse(gzlvl7, gt0); delay(gstab);
rgpulse(pw, one, rof1, rof1);
/* water defoc-refoc */
delay(gstab); zgradpulse(gzlvl8, gt8); delay(gstab);
/* t3 time */示例4: pulsesequence
//.........这里部分代码省略.........
decrgpulse(7.2*p_d,two,0.0,0.0);
decrgpulse(4.9*p_d,zero,0.0,0.0);
decrgpulse(7.4*p_d,two,0.0,0.0);
decrgpulse(6.8*p_d,zero,0.0,0.0);
decrgpulse(7.0*p_d,two,0.0,0.0);
decrgpulse(5.2*p_d,zero,0.0,0.0);
decrgpulse(5.4*p_d,two,0.0,0.0);
decrgpulse(0.6*p_d,zero,0.0,0.0);
decrgpulse(4.5*p_d,two,0.0,0.0);
decrgpulse(7.3*p_d,zero,0.0,0.0);
decrgpulse(5.1*p_d,two,0.0,0.0);
decrgpulse(7.9*p_d,zero,0.0,0.0);
decrgpulse(4.9*p_d,zero,0.0,0.0);
decrgpulse(7.9*p_d,two,0.0,0.0);
decrgpulse(5.0*p_d,zero,0.0,0.0);
decrgpulse(5.5*p_d,two,0.0,0.0);
decrgpulse(0.6*p_d,zero,0.0,0.0);
decrgpulse(4.6*p_d,two,0.0,0.0);
decrgpulse(7.2*p_d,zero,0.0,0.0);
decrgpulse(4.9*p_d,two,0.0,0.0);
decrgpulse(7.4*p_d,zero,0.0,0.0);
decrgpulse(6.8*p_d,two,0.0,0.0);
decrgpulse(7.0*p_d,zero,0.0,0.0);
decrgpulse(5.2*p_d,two,0.0,0.0);
decrgpulse(5.4*p_d,zero,0.0,0.0);
decrgpulse(0.6*p_d,two,0.0,0.0);
decrgpulse(4.5*p_d,zero,0.0,0.0);
decrgpulse(7.3*p_d,two,0.0,0.0);
decrgpulse(5.1*p_d,zero,0.0,0.0);
decrgpulse(7.9*p_d,two,0.0,0.0);
endhardloop();
dec2phase(zero);
decphase(zero);
txphase(zero);
decpwrf(rf10);
delay(tau2);
/* WFG3_START_DELAY */
sim3shaped_pulse("", "offC10", "", 2.0*pw, pwC10, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
if(pwC10>2.0*pwN) pwZ=0.0; else pwZ=2.0*pwN - pwC10;
delay(tau2);
decpwrf(rf0);
if (mag_flg[A] == 'y')
{
magradpulse(-icosel2*gzcal*gzlvl2, 1.8*gt1);
}
else
{
zgradpulse(-icosel2*gzlvl2, 1.8*gt1);
}
delay(2.02e-4);
decrgpulse(2.0*pwC, zero, 0.0, 0.0);
decpwrf(rf10);
if (mag_flg[A] == 'y')
{
magradpulse(icosel2*gzcal*gzlvl2, 1.8*gt1);
}
else
{
zgradpulse(icosel2*gzlvl2, 1.8*gt1);
}
delay(2.0e-4 + WFG3_START_DELAY + pwZ);示例5: pulsesequence
//.........这里部分代码省略.........
if (phase2 == 2) { tsadd(t3,1,4); }
if (d3_index % 2) { tsadd(t3,2,4); tsadd(t31,2,4); }
/* setup constant time in t2 (ni2) */
tau2 = d3;
t2max=2.0*(timeCN - CO180offCA.pw);
if((f2180[A] == 'y') && (ni2 > 0.0))
{tau2 += 0.5/sw2 ; t2max += 0.5/sw2 ;}
if(tau2 < 0.2e-6) {tau2 = 0.0;}
if ( (ni2-1.0)/sw2 > t2max)
{ text_error("too many ni2 increments in t2 ! "); psg_abort(1); }
if(FIRST_FID)
{
printf("t1max is %f\n",t1max);
printf("t2max is %f\n",t2max);
};
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
delay(d1);
zgradpulse(gzlvl2, gt2);
delay(gstab*3.0);
if (exp_mode[B]=='n') dec2rgpulse(2.0*pwN, zero, 0.0, 0.0); /* test for steady-state 15N */
/* Hz -> HzXz INEPT */
rgpulse(pw,zero,rof1,rof1); /* 1H pulse excitation */
zgradpulse(gzlvl0, gt0);
delay(tauNH -gt0);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab);
zgradpulse(gzlvl0, gt0);
delay(gstab);
rgpulse(pw, t6, rof1, rof1);
/* on HzNz now */
/* water flipback*/
obspower(tpwrs); obspwrf(tpwrsf_u);
shaped_pulse("H2Osinc",pwHs,zero,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
/* purge */
zgradpulse(gzlvl3, gt3);示例6: pulsesequence
//.........这里部分代码省略.........
*/
/* BEGIN PULSE SEQUENCE */
status(A);
delay(d1);
if ( dm3[B] == 'y' )
{ lk_hold(); lk_sampling_off();} /*freezes z0 correction, stops lock pulsing*/
rcvroff();
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
decpwrf(rf0);
obsoffset(tof);
decoffset(dofCO);
txphase(zero);
delay(1.0e-5);
if (TROSY[A] == 'n')
dec2rgpulse(pwN, zero, 0.0, 0.0); /*destroy N15 and C13 magnetization*/
decrgpulse(pwC, zero, 0.0, 0.0);
zgradpulse(-gzlvl0, 0.5e-3);
delay(1.0e-4);
if (TROSY[A] == 'n')
dec2rgpulse(pwN, one, 0.0, 0.0);
decrgpulse(pwC, zero, 0.0, 0.0);
zgradpulse(-0.7*gzlvl0, 0.5e-3);
delay(5.0e-4);
rgpulse(pw, zero, 0.0, 0.0); /* 1H pulse excitation */
dec2phase(zero);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
txphase(one);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0);
rgpulse(pw, one, 0.0, 0.0);
if (tpwrsf < 4095.0)
{obspwrf(tpwrsf);
obspower(tpwrs+6.0);} /* increases tpwrs by 6dB, now need */
else
obspower(tpwrs);
/* tpwrsf to be ~ 2048 for equivalence */
if (TROSY[A]=='y')
{txphase(two);
shaped_pulse("H2Osinc", pwHs, two, 5.0e-4, 0.0);
obspower(tpwr); obspwrf(4095.0);
zgradpulse(gzlvl3, gt3);
delay(2.0e-4);
dec2rgpulse(pwN, zero, 0.0, 0.0);
delay(0.5*kappa - 2.0*pw);
rgpulse(2.0*pw, two, 0.0, 0.0);
obspower(tpwrd); /* POWER_DELAY */
decphase(zero);
dec2phase(zero);示例7: pulsesequence
//.........这里部分代码省略.........
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
obsoffset(tof);
decoffset(dof); /* set Dec1 carrier at Co */
obspower(tsatpwr); /* Set transmitter power for 1H presaturation */
decpower(d_chirp); /* Set Dec1 power for hard 13C pulses */
dec2power(pwNlvl); /* Set Dec2 power for 15N hard pulses */
/* Presaturation Period */
if (fsat[0] == 'y')
{
delay(2.0e-5);
rgpulse(d1,zero,2.0e-6,2.0e-6); /* presaturation */
obspower(tpwr); /* Set transmitter power for hard 1H pulses */
delay(2.0e-5);
if (fscuba[0] == 'y')
{
delay(2.2e-2);
rgpulse(pw,zero,2.0e-6,0.0);
rgpulse(2*pw,one,2.0e-6,0.0);
rgpulse(pw,zero,2.0e-6,0.0);
delay(2.2e-2);
}
}
else
{
delay(d1);
}
obspower(tpwr); /* Set transmitter power for hard 1H pulses */
txphase(zero);
dec2phase(zero);
delay(1.0e-5);
/* Begin Pulses */
status(B);
rcvroff();
lk_hold();
delay(20.0e-6);
initval(1.0,v2);
obsstepsize(phase_sl);
xmtrphase(v2);
/* shaped pulse */
obspower(tpwrsl);
shaped_pulse(shp_sl,pw_sl,one,4.0e-6,0.0);
xmtrphase(zero);
obspower(tpwr); txphase(zero);
delay(4.0e-6);
/* shaped pulse */
rgpulse(pw,zero,0.0,0.0); /* 90 deg 1H pulse */
delay(0.2e-6);
zgradpulse(gzlvl5,gt5);
delay(2.0e-6);
delay(taua - gt5 - 2.2e-6); /* taua <= 1/4JNH */
sim3pulse(2*pw,0.0e-6,2*pwN,zero,zero,zero,0.0,0.0);
示例8: pulsesequence
//.........这里部分代码省略.........
if(t1_counter % 2) {
tsadd(t1,2,4);
tsadd(t5,2,4);
}
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
rlpower(tsatpwr,TODEV); /* Set transmitter power for 1H presaturation */
rlpower(dhpwr,DODEV); /* Set Dec1 power for 13C pulses */
rlpower(dhpwr2,DO2DEV); /* Set Dec2 power for 15N pulses */
obsoffset(tof);
/* Presaturation Period */
status(B);
if (fsat[0] == 'y')
{
delay(2.0e-5);
rgpulse(d1,zero,2.0e-6,2.0e-6); /* presaturation */
rlpower(tpwr,TODEV); /* Set transmitter power for hard 1H pulses */
delay(2.0e-5);
if(fscuba[0] == 'y')
{
delay(2.2e-2);
rgpulse(pw,zero,2.0e-6,0.0);
rgpulse(2*pw,one,2.0e-6,0.0);
rgpulse(pw,zero,2.0e-6,0.0);
delay(2.2e-2);
}
}
else
{
delay(d1);
}
rlpower(tpwr,TODEV); /* Set transmitter power for 1H CPMG pulses */
txphase(zero);
dec2phase(zero);
decphase(zero);
delay(1.0e-5);
/* Begin Pulses */
status(C);
rcvroff();
delay(20.0e-6);
decrgpulse(pwc,zero,4.0e-6,0.0);
delay(2.0e-6);
rgradient('z',gzlvl1);
delay(gt1);
rgradient('z',0.0);
delay(250.0e-6);
rgpulse(pw,zero,0.0,0.0);
decpower(d_reb);
delay(2.0e-6);
rgradient('z',gzlvl2);
delay(gt2);
rgradient('z',0.0);
delay(150.0e-6);
if(filter_flg[A] == 'y')
delay(taua - POWER_DELAY - gt2 - 152e-6
- WFG2_START_DELAY - 0.5*pwc_reb - 4.0/PI*pw);
else
delay(taua - POWER_DELAY - gt2 - 152e-6
- WFG2_START_DELAY - 0.5*pwc_reb);
simshaped_pulse("hard",sh_reb,2.0*pw,pwc_reb,zero,zero,0.0,0.0);
txphase(one);
decpower(dhpwr);
decphase(t4);
delay(taua - 0.5*pwc_reb - WFG2_STOP_DELAY - POWER_DELAY - gt2 - 152e-6 );
delay(2.0e-6);
rgradient('z',gzlvl2);示例9: pulsesequence
//.........这里部分代码省略.........
}
/* Set variables */
settable(t1, 2, phi1);
settable(t2, 1, phi2);
settable(t3, 2, rec);
/* PHASES AND INCREMENTED TIMES */
/* Phase incrementation for two fids to add */
if (psi2 == 2) tsadd(t2,2,4);
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (phase == 2) tsadd(t1,1,4);
/* Calculate modifications to phases for States-TPPI acquisition */
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t3,2,4); }
status(A);
delay(d0);
obsoffset(tof);
obspower(tpwrs);
txphase(zero);
decphase(zero);
decpower(pwClvl);
dec2phase(zero);
dec2power(pwNlvl);
delay(d1);
/* for a selective H20 one-lobe sinc pulse */
tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69)); /*needs 1.69 times more*/
tpwrs = (int) (tpwrs); /*power than a square pulse */
status(B);
dec2rgpulse(pwN,zero,rof1,rof1);
zgradpulse(gzlvl0,gt0);
delay(gstab);
obsstepsize(1.0);
xmtrphase(v8);
if (tpwrsf_d<4095.0)
{
obspower(tpwrs+6.0);
obspwrf(tpwrsf_d);
shaped_pulse("H2Osinc_d", pwHs, two, rof1, 0.0);
}
else
{
obspower(tpwrs);
shaped_pulse("H2Osinc_d", pwHs, two, rof1, 0.0);
}
xmtrphase(zero);
obspower(tpwr); obspwrf(4095.0);
示例10: pulsesequence
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pw1Hhytrap") + getval("pw2Hhytrap") + getval("tHX");
d.dutyoff = d1 + 4.0e-6 + getval("t1HYtrap") + getval("t2HYtrap");
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1Hhytrap,4,table1);
settable(phYhytrap,4,table2);
settable(ph2Hhytrap,4,table3);
settable(phXhx,4,table4);
settable(phHhx,4,table5);
settable(phRec,4,table6);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(ph1Hhytrap); dec2phase(phYhytrap);
obspwrf(getval("aXhx")); decpwrf(getval("aHhytrap")); dec2pwrf(getval("aYhytrap"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// TRAPDOR on H with Y Modulation
decrgpulse(getval("pw1Hhytrap"),ph1Hhytrap,0.0,0.0);
decphase(ph2Hhytrap);
decunblank();
dec2on();
delay(getval("t1HYtrap"));
dec2off();
decrgpulse(getval("pw2Hhytrap"),ph2Hhytrap,0.0,0.0);
decphase(phHhx);
decunblank();
decphase(phHhx);
decpwrf(getval("aHhx"));
delay(getval("t2HYtrap"));
// H to X Cross Polarization
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}示例11: pulsesequence
//.........这里部分代码省略.........
/* Set up f1180 */
tau1 = d2;
if((f1180[A] == 'y') && (ni > 1.0))
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1/2.0;
/* Calculate modifications to phases for States-TPPI acquisition */
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t4,2,4); }
kappa = tCN>ni/sw1? 1:tCN/(ni/sw1);
kappa = ((double)((int)(kappa*1000)))/1000;
/* 3 digits after the point, may cause a small decrease in value for safe boundary */
printf("kappa = %f\n", kappa);
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
decpwrf(rf0);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
delay(d1);
/* xxxxxxxxxxxxxxxxx CONSTANT SAMPLE HEATING FROM N15 RF xxxxxxxxxxxxxxxxx */
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx */
rcvroff();
/*dec2rgpulse(pwN, zero, 0.0, 0.0); */ /*destroy N15 magnetization*/
zgradpulse(gzlvl0, 0.5e-3);
delay(1.0e-4);
/*dec2rgpulse(pwN, one, 0.0, 0.0); */
zgradpulse(0.7*gzlvl0, 0.5e-3);
txphase(t1);
delay(5.0e-4);
rgpulse(pw,zero,0.0,0.0); /* 1H pulse excitation */
txphase(zero);
dec2phase(zero);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero,zero, zero, 0.0, 0.0);
txphase(one);
zgradpulse(gzlvl0, gt0);示例12: pulsesequence
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double aYxy8 = getval("aYxy8");
double pwYxy8 = getval("pwYxy8");
double nYxy8 = getval("nYxy8");
int cycles = (int) nYxy8/2.0;
nYxy8 = 2.0*cycles;
int counter = (int) (nYxy8 - 1.0);
initval((nYxy8 - 1.0),v8);
double onYxy8 = getval("onYxy8");
double srate = getval("srate");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
DSEQ mix = getdseq("Hmix");
strncpy(mix.t.ch,"dec",3);
putCmd("chHmixtppm='mix'\n");
strncpy(mix.s.ch,"dec",3);
putCmd("chHmixspinal='mix'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwX90") + 4.0*nYxy8*pwYxy8 + getval("pwX180");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d.c3 = d.c3 + (!strcmp(mix.seq,"tppm"));
d.c3 = d.c3 + ((!strcmp(mix.seq,"tppm")) && (mix.t.a > 0.0));
d.t3 = 2.0*nYxy8*(1.0/srate - 2.0*pwYxy8) + 1.0/srate - getval("pwX180");
d.c4 = d.c4 + (!strcmp(mix.seq,"spinal"));
d.c4 = d.c4 + ((!strcmp(mix.seq,"spinal")) && (mix.s.a > 0.0));
d.t4 = 2.0*nYxy8*(1.0/srate - 2.0*pwYxy8) + 1.0/srate - getval("pwX180");
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phX90,4,table1);
settable(ph1Yxy8,8,table2);
settable(ph2Yxy8,4,table3);
settable(phX180,4,table4);
settable(phRec,4,table5);
if (counter < 0) tsadd(phRec,2,4);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// X Single Pulse
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// xy8Y Period One
obspwrf(getval("aX180"));
txphase(phX180);
if (counter >= 0) {
_dseqon(mix);
delay(pwYxy8/2.0);
dec2pwrf(aYxy8);
sub(v1,v1,v1);
if (counter >= 1) {
if (counter > 1) loop(v8,v9);
getelem(ph1Yxy8,v1,v4);
incr(v1);
getelem(ph2Yxy8,ct,v2);
add(v4,v2,v2);
dec2phase(v2);
delay(0.5/srate - pwYxy8);
if (onYxy8 == 2)
dec2rgpulse(pwYxy8,v2,0.0,0.0);
else
delay(pwYxy8);
if (counter > 1) endloop(v9);
}
// X Refocussing Pulse
//.........这里部分代码省略.........
示例13: pulsesequence
//.........这里部分代码省略.........
tau3 = tau3/2.0;
flippwr = tpwr - 20.0*log10(flippw/(compH*pw*1.69));
flippwr = (int)(flippwr + 0.4);
if (fliphase < 0.2e-6) fliphase = fliphase + 360.0;
initval(fliphase, v10);
initval(7.0, v9);
obsstepsize(45.0);
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(dpwr);
dec2power(pwNlvl);
obsoffset(tofhn);
xmtrphase(v9);
delay(d1);
status(B); /* for 13C decoupling during t2, use dm='nynn' or 'nyyn' */
dec2rgpulse(pwN,zero,0.0,2.0e-6);
zgradpulse(gzlvl0,gt0);
rgpulse(pw,t1,1.0e-6,2.0e-6);
zgradpulse(gzlvl6,gt6);
delay(1.9*tauxh - gt1 - 2.0e-6);
dec2rgpulse(pwN, t2, 0.0, 0.0);
delay(tau2);
rgpulse(2.0*pw, t1, 0.0, 0.0);
dec2phase(zero);
delay(tau2);
dec2rgpulse(pwN, zero, 0.0, 0.0);
xmtrphase(zero);
delay(1.9*tauxh - gt6 - 600.0e-6 - SAPS_DELAY);
txphase(zero);
zgradpulse(gzlvl6,gt6);
txphase(zero);
delay(0.6e-3);
if (tau1 > pwN)
{
delay(tau1 - pwN);
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(tau1 - pwN);
}
else
{
delay(2.0*tau1);
}
rgpulse(pw, zero, 0.0, 1.0e-6);
status(A); /* no decoupling during mix period */
delay(mix - pwN - 1.5*gt7 - 2.0e-3);
zgradpulse(gzlvl7,gt7);
delay(1.0e-3);
dec2rgpulse(pwN, zero, 0.0, 2.0e-6);
zgradpulse(gzlvl7,gt7/2.0);
delay(1.0e-3 - 2.0e-6);示例14: pulsesequence
//.........这里部分代码省略.........
if( ix == 1) d3_init = d3;
t2_counter = (int) ( (d3-d3_init)*sw2 + 0.5 );
if(t2_counter % 2)
{ tsadd(t8,2,4); tsadd(t12,2,4); }
/* BEGIN PULSE SEQUENCE */
status(A);
delay(d1);
if (dm3[B]=='y') lk_hold();
rcvroff();
set_c13offset("ca");
obsoffset(tof);
obspower(tpwr);
obspwrf(4095.0);
decpower(pwClvl);
decpwrf(4095.0);
dec2power(pwNlvl);
txphase(three);
delay(1.0e-5);
dec2rgpulse(pwN, zero, 0.0, 0.0); /*destroy N15 and C13 magnetization*/
decrgpulse(pwC, zero, 0.0, 0.0);
zgradpulse(gzlvl0, 0.5e-3);
delay(1.0e-4);
dec2rgpulse(pwN, one, 0.0, 0.0);
decrgpulse(pwC, zero, 0.0, 0.0);
zgradpulse(0.7*gzlvl0, 0.5e-3);
delay(5.0e-4);
if(dm3[B] == 'y') /*optional 2H decoupling on */
{dec3unblank(); dec3rgpulse(1/dmf3, one, 0.0, 0.0);
dec3unblank(); setstatus(DEC3ch, TRUE, 'w', FALSE, dmf3);}
rgpulse(pw, zero, 0.0, 0.0); /* 1H pulse excitation */
txphase(zero);
decphase(zero);
zgradpulse(gzlvl0, gt0); /* 2.0*GRADIENT_DELAY */
decpwrf(rfst);
delay(tauCH - gt0 - WFG2_START_DELAY - 0.5e-3 + 70.0e-6);
simshaped_pulse("",stCshape, 2.0*pw, 1.0e-3, zero, zero, 0.0, 0.0);
delay(tauCH - gt0 - 0.5e-3 + 70.0e-6 - 150.0e-6);
decpwrf(rf0);
zgradpulse(gzlvl0, gt0); /* 2.0*GRADIENT_DELAY */
delay(150.0e-6);
rgpulse(pw, one, 0.0, 0.0);
zgradpulse(gzlvl3, gt3);
delay(2.0e-4);
decrgpulse(pwC, t1, 0.0, 0.0);
set_c13offset("co");
delay(zeta - 0.6*pwC - OFFSET_DELAY - POWER_DELAY - PWRF_DELAY - PRG_START_DELAY);
h1decon("DIPSI2", widthHd, 0.0); /*POWER_DELAY+PWRF_DELAY+PRG_START_DELAY */
delay(2.0*timeC - zeta);
c13pulse("co", "ca", "sinc", 90.0, t2, 0.0, 0.0); /* pwS1 */
delay(timeNCA - tau1);
c13pulse("ca", "co", "sinc", 180.0, zero, 2.0e-6, 2.0e-6);
sim3_c13pulse("", "co", "ca", "square", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0); /* pwS2 */
delay(timeNCA + tau1 + (60.0e-6));
initval(phshift, v3);
decstepsize(1.0);
dcplrphase(v3);
c13pulse("co", "ca", "sinc", 90.0, one, 0.0, 0.0); /* pwS1 */
delay(2.0*timeC + tauCC - OFFSET_DELAY - SAPS_DELAY - tau1);
c13pulse("ca", "co", "sinc", 180.0, zero, 0.0, 0.0);
delay(tauCC);
sim3_c13pulse("", "co", "ca", "square", "", 0.0, 180.0, 2.0*pwN, zero, zero, zero, 0.0, 60.0e-6);
delay(tau1);
set_c13offset("ca");
initval(phi7cal, v7);
decstepsize(1.0);
dcplrphase(v7); /* SAPS_DELAY */
dec2phase(t8);
nh_evol_se_train("ca", "co"); /* common part of sequence in bionmr.h */
if (dm3[B] == 'y') lk_sample();
} 示例15: pulsesequence
//.........这里部分代码省略.........
tau1 = tau1/2.0; tau2 = tau2/2.0; tau3 =tau3/2.0;
if (ix ==1 )
{
printf ("Current Spectral Width:\t\t%5.2f\n", swTilt);
printf ("Angle_H: %5.2f \n", angle_H);
printf ("Angle_C: %5.2f \n", angle_C);
printf ("Angle_N: %5.2f \n", angle_N);
printf ("\n\n\n\n\n");
}
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
delay(d1);
rcvroff();
obsoffset(tofali); obspower(tpwr); obspwrf(4095.0);
decoffset(dof); decpower(pwClvl); decpwrf(rf0);
dec2offset(dof2); dec2power(pwNlvl); dec2pwrf(4095.0);
if (gt0 > 0.2e-6)
{
decrgpulse(pwC,zero,10.0e-6,0.0);
dec2rgpulse(pwN,zero,2.0e-6,2.0e-6);
zgradpulse(gzlvl0,gt0);
delay(gstab);
}
txphase(t1); decphase(t2); dec2phase(zero);
status(B);
rgpulse(pw,t1,4.0e-6,2.0e-6);
zgradpulse(gzlvl3,gt3);
delay(2.0*tauCH - gt3 - 2.0*GRADIENT_DELAY -4.0e-6);
decrgpulse(pwC,t2,2.0e-6,2.0e-6);
decphase(zero);
/*======= Start of c13 evolution ==========*/
if ( ((tau2 -PRG_START_DELAY - POWER_DELAY -pwN - 2.0*pwC/PI -2.0e-6)> 0)
&& ((tau2 -PRG_STOP_DELAY - POWER_DELAY - pwN - 2.0*pwC/PI -2.0e-6)>0)
&& (codec[A] == 'y') )
{
decpower(copwr); decpwrf(copwrf);
decprgon(codecseq,1/codmf,cores);
decon();
delay(tau2 -PRG_START_DELAY - POWER_DELAY -pwN - 2.0*pwC/PI -2.0e-6);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, t1, zero, zero, 0.0, 0.0);
delay(tau2 -PRG_STOP_DELAY - POWER_DELAY - pwN - 2.0*pwC/PI -2.0e-6);
decoff();
decprgoff();
}
else if ( (tau2 -pwN - 2.0*pwC/PI -2.0e-6) > 0)
{
delay(tau2 -pwN - 2.0*pwC/PI -2.0e-6);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, t1, zero, zero, 0.0, 0.0);