SAS解方程怎么解，多谢了

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

用SAS解方程怎么解，比如解

w^2-u/s=0，其中

w=(1+b^2)/(1+(b-1)^2),

s= pi-sin(t)-2*t , u=pi-sin(v)-2*v,

v=atan(b-1) , t=atan(b)， pi 是圆周率

多谢了啊！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：解方程 Tan 圆周率 圆周率 解方程

Single variable equations are solvable in PROC FCMP

1. proc fcmp outlib=work.funcs.test;
   
2.         function huixingfu(b);
   
3.                 v = atan(b-1);
   
4.                 t = atan(b);
   
5.                 s = constant('pi')-sin(t)-2*t;
   
6.                 u = constant('pi')-sin(v)-2*v;
   
7.                 w = (1+b**2)/(1+(b-1)**2);
   
8.                 return (w**2-u/s);
   
9.         endsub;
   
10.         function s(b);
    
11.                 t = atan(b);
    
12.                 return (constant('pi')-sin(t)-2*t);
    
13.         endsub;
    
14.         function solveh(b0, status);
    
15.                 outargs status;
    
16.                 return (solve('huixingfu', {b0,.,.,.,status}, 0, .));
    
17.         endsub;
    
18.         function solves(b0, status);
    
19.                 outargs status;
    
20.                 return (solve('s', {b0,.,.,.,status}, 0, .));
    
21.         endsub;
    
22.         x = solve('huixingfu', {1,1e-12,1e-10,100,status}, 0, .);
    
23.         put x= status=;
    
24.         y = huixingfu(x);
    
25.         put y=;       
    
26.         x = solve('huixingfu', {2,1e-12,1e-10,100,status}, 0, .);
    
27.         put x= status=;
    
28.         y = huixingfu(x);
    
29.         put y=;       
    
30. run;
    
31. options cmplib=(work.funcs);
    
32. data huixingfu;
    
33.         do b=-5 to 5 by 0.01;
    
34.                 y = huixingfu(b);
    
35.                 output;
    
36.         end;
    
37.         status = .;
    
38.         b = solves(2, status);
    
39.         y = .;
    
40.         output;
    
41.         stop;
    
42. proc sort data=huixingfu;
    
43.         by b;
    
44. proc sgplot data=huixingfu;
    
45.         series x=b y=y / break;
    
46.         refline 0 / axis=y lineattrs=(color=blue);
    
47.         yaxis max=10 min=-10;
    
48.         xaxis values=(-5 to 5 by 1);
    
49. run;

复制代码

邓贵大 发表于 2013-10-21 22:43
Single variable equations are solvable in PROC FCMP

顺便问下邓兄
解多元方程的话用什么过程？

我刚查了下，线性方程组的话可以用IML
若是非线性的，只能求数值角的话，不知道用什么过程方便，谢谢！

playmore 发表于 2013-10-22 09:36
顺便问下邓兄
解多元方程的话用什么过程？

I think the official method is
http://support.sas.com/kb/22/779.html

If the functions are differentiable, you may write your own Newton-Raphson method in PROC FCMP or IML to solve an equation system.

好好学习一下，还没学过这个解方程的代码

wow值得学习下