菜菜路过，google查了下，实在不知道cos(c*P)应该怎么处理，那么我觉得至少我们可以折中的采用下面的方法
1：人为地给C赋值
2：做变量V1=s, V2=s*s, V3=cos(c*P)，V4 = v1*v3, V5= v2*v3
模型变为：Y=an*v3 + (bn-a)*v1*v3- b*v2*v3, 也就是Y=an*v3 + (bn-a)*v4 - b*v5,
然后跑跑模型proc logisitc ; Y= v3 v4 v5/dis=neg link=log;run;看看能得到什么系数，并将系数计算得到一组a, b, n c, 及方程的拟合程度。
3： 重复1，2, 将得到的a, b, n c, 及方程的拟合程度整理，比较，结合实际情况判断哪一组数据比较合适。

等待大神降临！