R怎么画函数y=tanx+siny+x的图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求在R中做一些超越函数以及比较复杂图的方法~谢谢

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：Tan 超越

这个 模拟就ok了
比如
x<-seq(-3,3,by=0.001)
  z<- sin(x)+tan(x)+x
  plot(x,z,type="l",ylim=range(-50,50))

看错了x,y 了， 再想想看

y=tan(x) + sin(y) +x
先想到的是 @nuomin 以前介绍的一个软件math...,呼唤 nuomin 看看能否用那个软件计算出 y.

另外就是迭代的形式看看能否接近。
y   = tan(x)          + sin(y)             +x
                         .......sin(tan(x) +sin(y) +x)
                         ..........................sin(tan(x) +sin(y) +x)
                        ...................................................

1. y=function(x,n=600){
   
2.     if(n==1) 0
   
3.     else tan(x)+sin(y(x,n-1))+x
   
4.   }
   
5.   x=seq(-2,2,by=0.01)
   
6.   y(x,600)
   
7.   plot(x,y(x,600),type="l",ylim=range(-100,100))
   
8.   

复制代码

x=0 ，y=0
x=1,   y(1,600)=2.847387 和 tan(1)+sin(2.847387) +1=2.847387
x=2,   y(2,600)=-1.054882 和 tan(2) +sin(-1.054882)+2= -1.054882
.......所以觉得能很好的接近。

欢迎讨论与补充

y=function(x,n=600){
    if(n==1) 0
    else tan(x)+sin(y(x,n-1))+x
  }
  x=seq(-50,50,by=0.01)
  y(x,600)
  plot(x,y(x,600),type="l",ylim=range(-500,500))
搞个图

看到jmpamao念叨我，就过来了

1. require(rootSolve)
   
2. x <- seq(-100,100,0.1)        #设置一个数组，避开pi/2的整数倍
   
3. xx <- x+tan(x)                #为求y-sin(y)的数值解准备方程右边的值
   
4. def_func <- function(x){
   
5.     oper <-ifelse(x>0,'-','+')
   
6.     eval(parse(text=paste("func_NS<- function(y) y-sin(y)",oper,abs(x),sep="")))
   
7.     uniroot(func_NS,c(-1000,1000))$root
   
8. }
   
9. 
   
10. yy <- sapply(xx,def_func)   #在每一个x+tan（x）值求y的数值解
    
11. plot(x,yy,type="l")         #画折线图
    
12. 
    
13. curve(x+tan(x),from=-100,to=100,n=2001,add=T,col="red",type="p")

复制代码

目测在|y|>100的时候图形可以画成y=tanx+x，误差不超过1%，这两者的误差如下图：
那些极值点都已经捕捉到了。

nuomin 发表于 2014-2-24 23:43
看到jmpamao念叨我，就过来了

哈哈 把你召唤来了
要像你学习这包。

jmpamao 发表于 2014-2-25 00:08
哈哈 把你召唤来了
要像你学习这包。

你用的递归思路很好，可以不用借助外部函数
我也发现了错误，已经修正好了，现在图形的趋势是向右上方倾斜。
如果想数据的点能取得密，y的取值范围要相应增大。

An alternative is to view the equation as a 2-variable function z=y-tan(x)-sin(y)-x and plot the solution of z=0:

1. my.fun <- function(x,y){y - sin(y) - x - tan(x)}
   
2. x<-seq(-50, 50, by = 0.02)
   
3. y<-seq(-500, 500, by = 0.2)
   
4. z<-outer(x,y,my.fun)
   
5. options(max.contour.segments=50000)
   
6. contour(x, y, z, levels=0, labels="", labcex=0, xlab = "x", ylab = "y")

复制代码

Results:
A big drawback is that a lot more data have to be generated than actually used in the plot.
@jmpamao The recursion method is brilliant.

zhangyangsmith 发表于 2014-2-25 15:35
An alternative is to view the equation as a 2-variable function z=y-tan(x)-sin(y)-x and plot the sol ...

请教
对contour用法不了解，是怎么判断 z为0

另外运行此程序后， 我的显示中间为 白板。怪~