本文整理汇总了C++中sim3pulse函数的典型用法代码示例。如果您正苦于以下问题：C++ sim3pulse函数的具体用法？C++ sim3pulse怎么用？C++ sim3pulse使用的例子？那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了sim3pulse函数的15个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: pulsesequence
//.........这里部分代码省略.........
delay(tauxh-gt3- WFG2_START_DELAY - 0.5e-3 + 70.0e-6 ); /* delay=1/4J(XH) */
if (TROSY[A] == 'y')
{
rgpulse(pw, v3, 0.0,0.0);
zgradpulse(-0.5*gzlvl3,gt3);
decpwrf(rf0); delay(200.0e-6);
decrgpulse(pwC, t2, 0.0, 0.0);
txphase(t1); decphase(zero);
if (tau1>0.0)
{
if ( (N15refoc[A]=='y') && (tau1 > (pwN/2.0 +2.0*pwC/PI) ) )
{delay(tau1 - pwN -2.0*pwC/PI);
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(tau1 - pwN -2.0*pwC/PI);}
else
{ if (tau1>2.0*pwC/PI)
delay(2.0*tau1-4.0*pwC/PI);
else
delay(2.0*tau1);
}
}
rgpulse(pw, t1, 0.0,0.0);
decpwrf(rfst);
zgradpulse(0.3*gzlvl3,gt3);
delay(tauxh-gt3-WFG2_START_DELAY - 0.5e-3 + 70.0e-6);
sim3shaped_pulse("","stC200","",2*pw,1.0e-3,0.0,zero,zero,zero,0.0,0.0); zgradpulse(0.3*gzlvl3,gt3);
delay(tauxh-gt3-WFG2_START_DELAY - 0.5e-3 + 70.0e-6);
sim3pulse(pw,pwC,0.0,two,zero,one,0.0,0.0);
zgradpulse(gzlvl3,gt3);
txphase(v4);
delay(tauxh-gt3-WFG2_START_DELAY - 0.5e-3 + 70.0e-6);
rgpulse(pw*0.231,v4,0.0,0.0);
delay(d3);
rgpulse(pw*0.692,v4,0.0,0.0);
delay(d3);
rgpulse(pw*1.462,v4,0.0,0.0);
delay(d3/2);
sim3shaped_pulse("","stC200","",0.0,1.0e-3,0.0,zero,zero,zero,0.0,0.0);
txphase(v1);
delay(d3/2);
rgpulse(pw*1.462,v1,0.0,0.0);
delay(d3);
rgpulse(pw*0.692,v1,0.0,0.0);
delay(d3);
rgpulse(pw*0.231,v1,0.0,0.0);
dec2phase(t3);
zgradpulse(gzlvl3,gt3); decpwrf(rf0);
delay(tauxh-gt3-WFG2_START_DELAY - 0.5e-3 + 70.0e-6);
decrgpulse(pwC, t3, 0.0,0.0);}
else
{
rgpulse(pw, t1, 0.0,0.0);
zgradpulse(0.5*gzlvl3,gt3); decpwrf(rf0);
delay(200.0e-6);示例2: pulsesequence
//.........这里部分代码省略.........
tsadd(t2,2,4);
tsadd(t14,2,4);
}
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
obspower(tpwr); /* Set power for pulses */
dec2power(pwNlvl); /* Set decoupler2 power to pwNlvl */
initval(ncyc+0.1,v10); /* for DIPSI-2 */
delay(d1);
status(B);
rcvroff();
/*destroy N15 magnetization*/
dec2rgpulse(pwN, zero, 0.0, 0.0);
zgradpulse(gzlvl1, gt1);
delay(9.0e-5);
/* 1H-15N INEPT */
rgpulse(pw,zero,1.0e-6,0.0);
txphase(zero); dec2phase(zero);
zgradpulse(gzlvl2,gt2);
delay(taua -pwN-0.5*pw -gt2 ); /* delay=1/4J(NH) */
sim3pulse(2.0*pw,0.0e-6,2*pwN,zero,zero,zero,0.0,0.0);
txphase(one); dec2phase(t1);
zgradpulse(gzlvl2,gt2);
delay(taua -1.5*pwN -gt2); /* delay=1/4J(NH) */
sim3pulse(pw,0.0e-6,pwN,one,zero,t1,0.0,0.0);
if (tpwrsf < 4095.0)
{obspwrf(tpwrsf); tpwrs=tpwrs+6.0;}
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, two, 2.0e-6, 0.0);
obspower(tpwr); obspwrf(4095.0); tpwrs=tpwrs-6.0;
txphase(zero); dec2phase(zero);
zgradpulse(gzlvl3,gt3);
delay(taub -1.5*pwN -gt3 -pwHs-2.0e-6-2.0*POWER_DELAY-WFG_START_DELAY); /* delay=1/4J(NH) */
sim3pulse(2.0*pw,0.0e-6,2*pwN,zero,zero,zero,0.0,0.0);
dec2phase(one);
zgradpulse(gzlvl3,gt3);
delay(taub -1.5*pwN -gt3 ); /* delay=1/4J(NH) */
dec2phase(one);
dec2rgpulse(pwN,one,0.0,0.0);
/* relaxation recovery */
if (ncyc>0.6)
{示例3: pulsesequence
//.........这里部分代码省略.........
delay(1.0e-4);
dec2rgpulse(pwN, one, 0.0, 0.0);
decrgpulse(pwC, one, 0.0, 0.0);
zgradpulse(0.7*gzlvl0, 0.5e-3);
txphase(t1);
decphase(zero);
dec2phase(zero);
initval(135.0,v1);
obsstepsize(1.0);
xmtrphase(v1);
delay(5.0e-4);
rcvroff();
rgpulse(pw, t1, 50.0e-6, 0.0); /* 1H pulse excitation */
xmtrphase(zero); /* SAPS_DELAY */
txphase(zero);
if (tau1 > (2.0*GRADIENT_DELAY + pwN + 0.64*pw + 5.0*SAPS_DELAY))
{
if (tau1>0.002)
{
zgradpulse(gzlvl6, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw) - SAPS_DELAY);
}
else
{
delay(tau1-pwN-0.64*pw);
}
if (C13refoc[A]=='y')
sim3pulse(0.0, 2.0*pwC, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
else
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
if (tau1>0.002)
{
zgradpulse(-1.0*gzlvl6, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw) - SAPS_DELAY);
}
else
{
delay(tau1-pwN-0.64*pw);
}
}
else if (tau1 > (0.64*pw + 0.5*SAPS_DELAY))
delay(2.0*tau1 - 2.0*0.64*pw - SAPS_DELAY );
rgpulse(pw, zero, 0.0, 0.0);
delay(mix - gt4 - gt5 -gstab -200.0e-6);
dec2rgpulse(pwN, zero, 0.0, 0.0);
zgradpulse(gzlvl4, gt4);
delay(gstab);
rgpulse(pw, zero, 200.0e-6,0.0); /* HSQC begins */
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);示例4: pulsesequence
pulsesequence()
{
double gzlvl1,
gt1,
gzlvl3,
gt3,
gstab,
pwx2lvl,
pwx2,
tau,
gtau,
hsglvl,
hsgt;
char nullflg[MAXSTR],
sspul[MAXSTR];
int iphase,
icosel;
gzlvl1 = getval("gzlvl1");
gt1 = getval("gt1");
gzlvl3 = getval("gzlvl3");
gt3 = getval("gt3");
pwx2lvl = getval("pwx2lvl");
pwx2 = getval("pwx2");
hsglvl = getval("hsglvl");
hsgt = getval("hsgt");
gstab = getval("gstab");
getstr("nullflg",nullflg);
getstr("sspul",sspul);
iphase = (int)(getval("phase")+0.5);
icosel = 1;
tau = 1/(2*(getval("j1xh")));
gtau = 2*gstab + GRADIENT_DELAY;
if (tau < (gt3+gstab))
{
text_error("tau must be greater than gt3+gstab\n");
psg_abort(1);
}
settable(t1,2,ph1);
settable(t2,4,ph2);
settable(t3,4,ph3);
assign(zero,v4);
getelem(t1,ct,v1);
getelem(t3,ct,oph);
if (iphase == 2)
icosel = -1;
initval(2.0*(double)((int)(d2*getval("sw1")+0.5)%2),v10);
add(v1,v10,v1);
add(v4,v10,v4);
add(oph,v10,oph);
status(A);
dec2power(pwx2lvl);
if (sspul[0] == 'y')
{
zgradpulse(hsglvl,hsgt);
rgpulse(pw,zero,rof1,rof1);
zgradpulse(hsglvl,hsgt);
}
delay(d1);
status(B);
rcvroff();
if (nullflg[0] == 'y')
{
rgpulse(0.5*pw,zero,rof1,rof1);
delay(tau);
sim3pulse(2.0*pw,0.0,2.0*pwx2,zero,zero,zero,rof1,rof1);
delay(tau);
rgpulse(1.5*pw,two,rof1,rof1);
zgradpulse(hsglvl,hsgt);
delay(1e-3);
}
rgpulse(pw,zero,rof1,rof1);
delay(tau - rof1 - (2*pw/PI));
dec2rgpulse(pwx2,v1,rof1,1.0e-6);
delay(gt1+gtau - (2*pwx2/PI) - pwx2 - 1.0e-6 - rof1);
dec2rgpulse(2*pwx2,v4,rof1,1.0e-6);
delay(gstab - pwx2 - 1.0e-6);
zgradpulse(gzlvl1,gt1);
delay(gstab - rof1 - pw);
delay(d2/2);
rgpulse(2.0*pw,zero,rof1,rof1);
delay(d2/2);
delay(gstab - rof1 - pw);
zgradpulse(gzlvl1,gt1);
delay(gstab - pwx2 - rof1);
dec2rgpulse(2*pwx2,zero,rof1,1.0e-6);
delay(gt1+gtau - (2*pwx2/PI) - pwx2 - 1.0e-6 - rof1);
//.........这里部分代码省略.........
示例5: pulsesequence
//.........这里部分代码省略.........
tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69)); /*needs 1.69 times more*/
tpwrs = (int) (tpwrs); /*power than a square pulse */
if (tpwrsf_d<4095.0) tpwrs=tpwrs+6; /* nominal tpwrsf_d ~ 2048 */
/* tpwrsf_d,tpwrsf_u can be used to correct for radiation damping */
/* BEGIN ACTUAL PULSE SEQUENCE */
status(A);
obspower(satpwr); /* Set power for presaturation */
decpower(pwClvl); /* Set decoupler1 power to pwClvl */
decpwrf(rf90onco);
dec2power(pwNlvl); /* Set decoupler2 power to pwNlvl */
/* Presaturation Period */
if (satmode[0] == 'y')
{
rgpulse(d1,zero,rof1,0.0);
obspower(tpwr); /* Set power for hard pulses */
}
else
{
obspower(tpwr); /* Set power for hard pulses */
delay(d1);
}
status(B);
rcvroff();
decpwrf(rf90onco); /* Set decoupler1 power to rf90onco */
sim3pulse(0.0,pw90onco,pwN,zero,zero,zero,rof1,rof1); /* 90 for 15N and 13C' */
zgradpulse(gzlvl0,gt0);
delay(gstab);
/* transfer from HN to N by INEPT */
/* shaped pulse for water flip-back */
obspower(tpwrs); obspwrf(tpwrsf_d);
shaped_pulse("H2Osinc_d",pwHs,one,2.0e-6,0.0);
obspower(tpwr); obspwrf(4095.0);
/* shaped pulse */
rgpulse(pw,zero,rof1,0.0);
zgradpulse(gzlvl0*1.3,gt0);
delay(gstab);
delay(tauhn - gt0 - gstab); /* 1/(4JHN) */
sim3pulse(2.0*pw,0.0,2.0*pwN,zero,zero,zero,rof1,rof1);
delay(tauhn - gt0 - gstab); /* 1/(4JHN) */
zgradpulse(gzlvl0*1.3,gt0);
delay(gstab);
rgpulse(pw,three,rof1,0.0); /* 90 1H */
zgradpulse(gzlvl3,gt3);
delay(gstab);
decpwrf(rf180onco); /* Set decoupler power to rf180onco */示例6: pulsesequence
//.........这里部分代码省略.........
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);
delay(tau2 -pwN - 2.0*pwC/PI -2.0e-6);
}
else
{
delay(2.0*tau2);
decphase(t14);
delay(4.0e-6);
sim3pulse(2.0*pw, 2.0*pwC, 2.0*pwN, t1, t14, zero, 0.0, 0.0);
delay(4.0e-6);
}
decpower(pwClvl);
decpwrf(rf0);
decphase(zero);
/*======= End of c13 evolution ==========*/
decrgpulse(pwC,zero, 2.0e-6,2.0e-6);
txphase(zero);
delay(2.0*tauCH + tau1 - gt3 - 4.0*pwC
- gstab -4.0e-6 - 2.0*GRADIENT_DELAY);
zgradpulse(gzlvl3,gt3);
delay(gstab);
示例7: pulsesequence
//.........这里部分代码省略.........
}
}
else
{
delay(d1);
}
obspower(tpwr); /* Set transmitter power for hard 1H pulses */
txphase(one);
dec2phase(zero);
delay(1.0e-5);
/* Begin Pulses */
status(B);
rcvroff();
lk_hold();
shiftedpulse("sinc", pwHs, 90.0, 0.0, one, 2.0e-6, 0.0);
txphase(zero);
delay(2.0e-6);
/* xxxxxxxxxxxxxxxxxxxxxx 1HN to 15N TRANSFER xxxxxxxxxxxxxxxxxx */
rgpulse(pw,zero,0.0,0.0); /* 90 deg 1H pulse */
delay(0.2e-6);
zgradpulse(gzlvl1, gt1);
delay(2.0e-6);
delay(taua - gt1 - 2.2e-6); /* taua <= 1/4JNH */
sim3pulse(2*pw,0.0,2*pwN,zero,zero,zero,0.0,0.0);
txphase(three); dec2phase(zero); decphase(zero);
delay(taua - gt1 - gstab -0.2e-6 - 2.0e-6);
delay(0.2e-6);
zgradpulse(gzlvl1, gt1);
delay(gstab);
/* xxxxxxxxxxxxxxxxxxxxxx 15N to 13CO TRANSFER xxxxxxxxxxxxxxxxxx */
if(sel_flg[A] == 'n') {
rgpulse(pw,three,2.0e-6,0.0);
delay(0.2e-6);
zgradpulse(gzlvl2, gt2);
delay(gstab);
dec2rgpulse(pwN,zero,0.0,0.0);
delay( zeta + pwS4 );
dec2rgpulse(2*pwN,zero,0.0,0.0);
c13pulse("co", "ca", "sinc", 180.0, zero, 0.0, 0.0);
dec2phase(one);
delay(zeta - 2.0e-6);
dec2rgpulse(pwN,one,2.0e-6,0.0);
}示例8: pulsesequence
//.........这里部分代码省略.........
status(B);
decpower(dpwr);
delay(tau2);
rgpulse(2.0*pw,t1,0.0,0.0);
decphase(zero);
if (tau2 > 2.0*pwN)
{
dec2rgpulse(2.0*pwN,zero,0.0,0.0);
delay(tau2 - 2.0*pwN);
}
else
delay(tau2);
status(A);
decpower(pwClvl);
decrgpulse(pwC,zero, 2.0e-6,2.0e-6);
txphase(zero);
delay(tauch + tau1 + SAPS_DELAY - gt3 - 4.0*pwC - 500.0e-6);
zgradpulse(gzlvl3,gt3);
delay(500.0e-6);
decrgpulse(pwC,zero,0.0, 0.0);
decphase(one);
decrgpulse(2.0*pwC, one, 0.2e-6, 0.0);
decphase(zero);
decrgpulse(pwC, zero, 0.2e-6, 0.0);
delay(tau1);
rgpulse(pw,zero,0.0,0.0);
delay(mix - gt4 - gt5 - pwN - 2.0e-3);
obsoffset(tof);
zgradpulse(gzlvl4,gt4);
delay(1.0e-3);
sim3pulse((double)0.0,pwC,pwN,zero,zero,zero,0.0,2.0e-6);
zgradpulse(gzlvl5,gt5);
delay(1.0e-3);
rgpulse(pw,zero,0.0,2.0e-6);
zgradpulse(gzlvl6,gt6);
delay(taunh - gt6 - 2.0e-6);
sim3pulse(2.0*pw,(double)0.0,2*pwN,zero,zero,zero,0.0,0.0);
delay(taunh - gt6 - 500.0e-6);
zgradpulse(gzlvl6,gt6);
txphase(one);
delay(500.0e-6);
rgpulse(pw,one,0.0,0.0);
txphase(two);
obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, two, 2.0e-6, 2.0e-6);
obspower(tpwr);
zgradpulse(gzlvl7,gt7);
dec2phase(t3);
decoffset(dofcaco);
decpwrf(rfst);
delay(200.0e-6);
dec2rgpulse(pwN,t3,0.0,0.0);
dec2phase(t4);
delay(tau3);
rgpulse(2.0*pw, zero, 0.0, 0.0);
decshaped_pulse("stC200", 1.0e-3, zero, 0.0, 0.0);
delay(tau3);
delay(gt1 + 202.0e-6 - 1.0e-3 - 2.0*pw);
dec2rgpulse(2.0*pwN, t4, 0.0, 2.0e-6);示例9: pulsesequence
//.........这里部分代码省略.........
icosel = -1;
if (phase1 == 1) {tsadd(t10,2,4); icosel = +1;}
/* 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(d2_index % 2) { tsadd(t1,2,4); tsadd(t31,2,4); }
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
delay(d1);
status(B);
dec2rgpulse(pwN, zero, 0.0, 0.0); /*destroy N15 magnetization*/
zgradpulse(gzlvl0, 0.5e-3);
delay(1.0e-4);
rgpulse(pw,zero,rof1,rof1); /* 1H pulse excitation */
delay(gstab);
zgradpulse(gzlvl0, gt0);
delay(lambda - gt0 -gstab);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
txphase(one);
delay(lambda - gt0 -gstab);
zgradpulse(gzlvl0, gt0);
delay(gstab);
rgpulse(pw, one, rof1, rof1); /* on NzHz now */
if(water_sat[A]=='n') /* water to -Z */
{
obspower(tpwrs); obspwrf(tpwrsf_d);
shaped_pulse("H2osinc",pwHs,two,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
}
/* purge */
zgradpulse(gzlvl3, gt3);
delay(gstab);
/* HzNz-> Nz */
dec2rgpulse( pwN, zero, 0.0, 0.0);
delay(gstab);
zgradpulse(gzlvl6, gt0);
delay(lambda - gt0 -gstab);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
dec2phase(t3);
delay(lambda - gt0 -gstab);
zgradpulse(gzlvl6, gt0);
delay(gstab);
dec2rgpulse( pwN, t3, 0.0, 0.0);
/* purge */
zgradpulse(gzlvl7, gt7);
dec2phase(t3);
delay(gstab);示例10: pulsesequence
pulsesequence()
{
double pwx2lvl,
pwx2,
hsglvl,
hsgt,
satfrq,
satdly,
satpwr,
tau,
j1xh,
null,
phase;
int iphase;
char sspul[MAXSTR],
nullflg[MAXSTR],
PFGflg[MAXSTR],
satmode[MAXSTR];
pwx2lvl = getval("pwx2lvl");
pwx2 = getval("pwx2");
hsglvl = getval("hsglvl");
hsgt = getval("hsgt");
getstr("PFGflg",PFGflg);
getstr("nullflg",nullflg);
satfrq = getval("satfrq");
satdly = getval("satdly");
satpwr = getval("satpwr");
null = getval("null");
j1xh = getval("j1xh");
tau = 1/(4*j1xh);
phase = getval("phase");
getstr("satmode",satmode);
getstr("sspul",sspul);
iphase = (int) (phase + 0.5);
settable(t1,4,ph1);
settable(t2,2,ph2);
settable(t3,8,ph3);
settable(t4,16,ph4);
settable(t5,16,ph5);
getelem(t2,ct,v2);
getelem(t5,ct,oph);
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v14);
if (iphase == 2)
incr(v2);
add(v2,v14,v2);
add(oph,v14,oph);
status(A);
dec2power(pwx2lvl);
obspower(tpwr);
if (sspul[0] == 'y')
{
if (PFGflg[0] == 'y')
{
zgradpulse(hsglvl,hsgt);
rgpulse(pw,zero,rof1,rof1);
zgradpulse(hsglvl,hsgt);
}
else
{
obspower(tpwr-12);
rgpulse(500*pw,zero,rof1,rof1);
rgpulse(500*pw,one,rof1,rof1);
obspower(tpwr);
}
}
delay(d1);
if (satmode[0] == 'y')
{
obspower(satpwr);
if (satfrq != tof)
obsoffset(satfrq);
rgpulse(satdly,zero,rof1,rof1);
if (satfrq != tof)
obsoffset(tof);
obspower(tpwr);
delay(1.0e-5);
}
status(B);
if (PFGflg[0] == 'y')
{
if (nullflg[0] == 'y')
{
rgpulse(0.5*pw,zero,rof1,rof1);
delay(2*tau);
sim3pulse(2.0*pw,0.0,2.0*pwx2,zero,zero,zero,rof1,rof1);
delay(2*tau);
//.........这里部分代码省略.........
示例11: pulsesequence
//.........这里部分代码省略.........
initval(ncyc_nc,v11);
txphase(t1);
decphase(zero);
dec2phase(zero);
delay(5.0e-4);
rcvroff();
if(pw11[A] == 'y')
{
rgpulse(pw/2, t1, 50.0e-6, 0.0);
if (URA[A] == 'y')
delay(1/(2*(tof_12-tof)));
else if (CYT[A] == 'y')
delay(1/(2*(tof_75-tof)));
else
delay(1/(2*(tof_75-tof)));
rgpulse(pw/2, t1, 0.0, 0.0);
}
else
{
rgpulse(pw, t1, 50.0e-6, 0.0);
}
txphase(zero);
if (C13refoc[A]=='y')
{
if (tau1 > (0.001-(2.0*GRADIENT_DELAY + pwN + 0.64*pw )))
{
zgradpulse(gzlvlr, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw) );
sim3pulse(0.0, 2.0*pwC, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(-1.0*gzlvlr, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
}
else
{
if (tau1 > (0.001-(pwN + 0.64*pw )))
{
delay(tau1 - pwN - 0.64*pw );
sim3pulse(0.0, 2.0*pwC, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tau1 - pwN - 0.64*pw);
}
else
{
if (tau1 > (0.64*pw ))
delay(2.0*tau1 - 2.0*0.64*pw );
}
}
}
else
{
if (tau1 > (0.001-(2.0*GRADIENT_DELAY + pwN + 0.64*pw )))
{
zgradpulse(gzlvlr, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw) );
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
zgradpulse(-1.0*gzlvlr, 0.8*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
delay(0.2*(tau1 - 2.0*GRADIENT_DELAY - pwN - 0.64*pw));
}
else
{
if (tau1 > (0.001-(pwN + 0.64*pw )))
{示例12: pulsesequence
//.........这里部分代码省略.........
dec2power(pwNlvl); /* Set decoupler2 power to pwNlvl */
obsoffset(tof);
decoffset(dof);
dec2offset(dof2);
/* Presaturation Period */
if (satmode[0] == 'y')
{
rgpulse(d1,zero,rof1,rof1);
obspower(tpwr); /* Set power for hard pulses */
}
else
{
obspower(tpwr); /* Set power for hard pulses */
delay(d1);
}
if (tpwrsf_d<4095.0) tpwrs=tpwrs+6.0; /* allow for fine power control via tpwrsf_d */
rcvroff();
dec2phase(zero);
obspwrf(tpwrsf_d); obspower(tpwrs);
shaped_pulse("H2Osinc_d", pwHs, one, rof1, rof1);
obspower(tpwr); obspwrf(4095.0);
rgpulse(pw,zero,rof1,0.0);
zgradpulse(gzlvl0,gt0);
delay(gstab);
delay(tauhn - gt0 - gstab); /* delay for 1/4JHN coupling */
sim3pulse(2.0*pw,0.0,2.0*pwN,zero,zero,zero,rof1,rof1);
delay(tauhn - gt0 - gstab); /* delay for 1/4JHN coupling */
zgradpulse(gzlvl0,gt0);
delay(gstab -rof1);
dec2phase(zero);
rgpulse(pw,three,rof1,0.0);
zgradpulse(gzlvl3,gt3);
delay(gstab);
dec2rgpulse(pwN,zero,0.0,0.0);
if (abfilter[0] == 'b')
{
zgradpulse(gzlvl5*1.2,gt5);
delay(gstab);
delay(0.5*taunco - gt5 - gstab); /* 1/8J(NCO) */
sim3pulse(0.0*pw,0.0*pw180onco,0.0,zero,zero,zero,rof1,rof1);
delay(0.5*taunco - gt5 - gstab); /* 1/8J(NCO) */
zgradpulse(gzlvl5*1.2,gt5);
dec2phase(t4);
delay(gstab);
sim3pulse(0.0,pw180onco,2.0*pwN,zero,zero,t4,rof1,rof1);示例13: pulsesequence
//.........这里部分代码省略.........
/* H1 EVOLUTION BEGINS */
txphase(t3); dec2phase(zero);
if (C13refoc[A]=='y' && (tau1 -2.0*pwC -3.0*SAPS_DELAY) > 0.2e-6)
{
delay(tau1 -2.0*pwC -3.0*SAPS_DELAY);
decrgpulse(pwC, zero, 0.0, 0.0);
decphase(one);
decrgpulse(2.0*pwC, one, 0.0, 0.0);
decphase(zero);
decrgpulse(pwC, zero, 0.0, 0.0);
delay(tau1 -2.0*pwC -SAPS_DELAY);
}
else if (2.0*tau1 -2.0*SAPS_DELAY > 0.2e-6)
delay(2.0*tau1 -2.0*SAPS_DELAY);
/* H1 EVOLUTION ENDS */
rgpulse(pw, t3, 0.0, 0.0);
initval(phincr_a,v14);
if (phincr_a != 0) xmtrphase(v14);
if (tpwrsf_a < 4095.0) { obspwrf(tpwrsf_a); obspower(tpwrs+6.0); }
else obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, t3, 2.0e-6, 0.0);
obspower(tpwr);
if (tpwrsf_a < 4095.0) obspwrf(4095.0);
txphase(zero); if (phincr_a != 0) xmtrphase(zero);
zgradpulse(gzlvl4,gt4);
delay(taua -pwN -0.5*pw -gt4 -2.0*GRADIENT_DELAY -pwHs_dly); /* delay=1/4J(NH) */
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
dec2phase(t3);
zgradpulse(gzlvl4,gt4);
delay(taua -1.5*pwN -gt4 -2.0*GRADIENT_DELAY -pwHs_dly); /* delay=1/4J(NH) */
initval(phincr_d,v14);
txphase(two); if (phincr_d != 0) xmtrphase(v14);
if (tpwrsf_d < 4095.0) { obspwrf(tpwrsf_d); obspower(tpwrs+6.0); }
else obspower(tpwrs);
shaped_pulse("H2Osinc", pwHs, two, 2.0e-6, 0.0);
obspower(tpwr);
if (tpwrsf_d < 4095.0) obspwrf(4095.0);
txphase(zero); if (phincr_d != 0) xmtrphase(zero);
sim3pulse(pw, 0.0, pwN, zero, zero, t3, rof1, rof1);
dec2phase(zero);
zgradpulse(gzlvl5,gt5);
delay(taua -1.5*pwN -wtg_dly -gt5 -2.0*GRADIENT_DELAY);
if(wtg3919[0] == 'y') /* 3919 watergate */
{
txphase(v11);
rgpulse(pw*0.231,v11,rof1,rof1);
delay(d3919);
rgpulse(pw*0.692,v11,rof1,rof1);
delay(d3919);
rgpulse(pw*1.462,v11,rof1,rof1);
delay(d3919/2.0-pwN);
dec2rgpulse(2.0*pwN, zero, rof1, rof1);
txphase(v12);示例14: 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);
txphase(zero);
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);
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 (TROSY[A]=='y')
{
txphase(two);
obspower(tpwrs);
shaped_pulse("H2Osinc",pwHs,two,5.0e-4,0.0);
obspower(tpwr);
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);
dec2phase(zero);
decpwrf(rf2);
delay(timeTN - 0.5*kappa);
}
else
{
txphase(zero);
obspower(tpwrs);
shaped_pulse("H2Osinc",pwHs,zero,5.0e-4,0.0);
obspower(tpwrd);
zgradpulse(gzlvl3, gt3);示例15: pulsesequence
//.........这里部分代码省略.........
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(rfC);
delay(tau2);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tau2);
decpwrf(rfC);
zgradpulse(-icosel2*gzlvl2, 1.8*gt1);
delay(grecov+2.0e-6);
decrgpulse(2.0*pwC, zero, 0.0, 0.0);
decpwrf(rfC);
zgradpulse(icosel2*gzlvl2, 1.8*gt1);
delay(grecov + pwN);
decrgpulse(pwC, zero, 0.0, 0.0);
decpwrf(rfC);
decrgpulse(pwC, zero, 2.0e-6, 0.0);
zgradpulse(gzlvl5, gt5);
delay(0.5*del1 - gt5);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl5, gt5);
txphase(one);
decphase(t10);
delay(0.5*del1 - gt5);
simpulse(pw, pwC, one, t10, 0.0, 0.0);
zgradpulse(gzlvl6, gt5);
txphase(zero);
decphase(zero);
delay(0.5*del2 - gt5);