注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

阿弥陀佛

街树飘影未见尘 潭月潜水了无声 般若观照心空静...

 
 
 

日志

 
 
关于我

一直从事气象预报、服务建模实践应用。 注重气象物理场、实况场、地理信息、本体知识库、分布式气象内容管理系统建立。 对Barnes客观分析, 小波,计算神经网络、信任传播、贝叶斯推理、专家系统、网络本体语言有一定体会。 一直使用Java、Delphi、Prolog、SQL编程。

网易考拉推荐

scala 二维二阶切比雪夫多项式  

2013-08-31 17:27:15|  分类: Scala |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
package chebList
/**
 * 二维2阶切比雪夫多项式
 * 何险峰
 * 2013-08-31
 */
import Array._
case class D2Cheb(f: Array[Array[Float]]) {
  val x = CbVal(f(0).length)
  val y = CbVal(f.length)
  //二维要素的多项式权重
  val A = ofDim[Float](y.K0 + 1, x.K0 + 1)
  for (
    k <- (0 to y.K0).par; s <- (0 to x.K0).par;
    i <- (0 until y.I0).par; j <- (0 until x.I0).par
  ) {
    A(k)(s) += f(i)(j) * y.v(k)(i) * x.v(s)(j)
  }

  //二维要素拟合
  val fit = ofDim[Float](y.I0, x.I0)
  for (
    i <- (0 until y.I0).par; j <- (0 until x.I0).par;
    k <- (0 to y.K0).par; s <- (0 to x.K0).par
  ) {
    fit(i)(j) += A(k)(s) * y.v(k)(i) * x.v(s)(j)

  }
  //二维要素拟合
  def inverse(cols: Int, rows: Int, AA: Array[Array[Float]]): Array[Array[Float]] = {
    val xx = CbVal(cols)
    val yy = CbVal(rows)
    val z = ofDim[Float](yy.I0, xx.I0)
    for (
      i <- (0 until yy.I0).par; j <- (0 until xx.I0).par;
      k <- (0 to yy.K0).par; s <- (0 to xx.K0).par
    ) {
      z(i)(j) += AA(k)(s) * yy.v(k)(i) * xx.v(s)(j)
    }
    z
  }
}
  评论这张
 
阅读(537)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017