JPH0555637B2 - - Google Patents
Info
- Publication number
- JPH0555637B2 JPH0555637B2 JP62263093A JP26309387A JPH0555637B2 JP H0555637 B2 JPH0555637 B2 JP H0555637B2 JP 62263093 A JP62263093 A JP 62263093A JP 26309387 A JP26309387 A JP 26309387A JP H0555637 B2 JPH0555637 B2 JP H0555637B2
- Authority
- JP
- Japan
- Prior art keywords
- basis weight
- dry basis
- weight profile
- absolute dry
- profile
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 claims description 27
- 238000005259 measurement Methods 0.000 claims description 23
- 210000000988 bone and bone Anatomy 0.000 claims description 7
- 238000011156 evaluation Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 10
- 239000002994 raw material Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 240000001931 Ludwigia octovalvis Species 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 238000004260 weight control Methods 0.000 description 1
Landscapes
- Paper (AREA)
Description
【発明の詳細な説明】
＜産業上の利用分野＞
本発明は、抄紙機の制御システムにおいて、操
作端であるスライスボルトの位置に対応する絶乾
坪量プロフイル測定点の位置を精度よく決定する
方法に関する。 ＜従来技術＞ 第５図は、抄紙プロセスの概要図である。原料
パルプＰは種箱１、種口弁２を経て白水サイロ３
からのリターン原料と合流してポンプ４によりヘ
ツドボツクス５に供給される。ヘツドボツクス内
の原料は、スライスリツプ６の間隙よりワイヤー
パート８上にシート状に吐出され、プレスパート
９で搾水され、ドライパート１０で蒸気により加
熱乾燥された後、カレンダー１１を経てリール１
２に巻き取られて製品となる。７はスライスリツ
プの間隙開度を操作するための複数（例えば40
本）のスライスボルトである。 １３は製品を紙幅方向にスキヤンして紙の坪量
（１m2の重さ）と水分率を測定するＢ／Ｍ計であ
り、この測定値に基づいて紙の絶乾坪量プロフイ
ルが紙の横幅方向を複数点、例えば60点又は360
点に分割して計算される。 １４は、計算で求められた各点の絶乾坪量プロ
フイル測定値PVとプロフイル設定値SVの差を制
御演算する制御装置であり、その操作出力MVに
より測定点に対応するスライスボルト７が駆動さ
れ、スライスリツプ６の開度が操作され、製品の
絶乾坪量プロフイルが設定絶乾坪量プロフイルに
一致するようなフイードバツク制御が実行され
る。 ＜発明が解決しようとする問題点＞ 初期の絶乾坪量プロフイル制御システムでは、
ヘツドボツクス５よりスライスリツプ６を介して
吐出した原料は、第６図Ａに示す様にヘツドボツ
クスのｉ番スライスボルト位置よりまつすぐに飛
び出し、ワイヤー上をまつすぐ流れ、乾燥段階で
一様に収縮するとしてｊ番プロフイル測定点位置
に対応するものとし、Ｂに示すようにプロフイル
の測定値PVと設定値SVの偏差に基づいて対応す
るスライスボルトを操作して絶乾坪量プロフイル
の制御を実行している。Ｃはこのように制御され
た場合のスライスリツプ開度プロフイルの例であ
る。 所が、実際には、原料はヘツドボツクスからま
つすぐに吐出するとは限らないし、流れもまつす
ぐではない。又、収縮率も紙端近傍のほうが中央
部より大きく、現実の流れによる対応点は、ｊ番
目からずれてｌ番のプロフイル測定点となつてい
る場合が多い。 このように、実際の測定点と操作端の位置対応
のずれが大きいと、絶乾坪量プロフイル制御にお
いて好ましくないプロフイルの鋸歯状波を生ずる
ことがある。 出願人は、このような問題点を解消する対応点
決定手法として４次以上の回帰モデルによる最小
自乗法による対応点の決定手法を特願昭62−
264917号で提案している。 第７図によりその概要を説明する。ｉ番スライ
スボルトを操作したときに、抄紙機制御装置１４
からは例えば60点の幅方向のプロフイル測定点の
平滑化されたＮスキヤン分の複数の時系列デー
タ、1 1〜N 1、1 2〜N 2…1 60〜N 60がオンライ
ン
出力され、Ｉ／Ｏコントローラ１５を介してデー
タ記録装置１６にBD（絶乾坪量）データとして
フロツピーデイスクなどの記憶手段に記録保持さ
れる。 このBDデータは、オンラインでＩ／Ｏコント
ローラ１７を介してデータ処理装置１８のBDデ
ータ読出し部１９により読み出され、定常状態と
見なされるＮスキヤンのCD方向（幅方向）のBD
データの収束値N 1、N 2…N 60を20で４次又は６
次の回帰モデル、 N n(r)＝4 〓i=0 airi または、 N n(r)＝6 〓i=0 airi により最小２乗処理する。ここでｒはCO方向の
距離、ｂはBD出力の平滑化された値である。 第８図及び第９図は、この処理結果得られた最
小２乗４次及び６次の近似カーブであり、そのピ
ーク値位置がピーク値決定部２１により求めら
れ、この位置（ｊ＋１）番をｉ番スライスボルト
の対応測定点として出力部２２より出力し、オフ
ライン的に抄紙機制御装置１４の対応点データを
修正する。 ＜発明が解決しようとする問題点＞ このような最小２乗法４次又は６次の回帰モデ
ルによるピーク値位置を求める手法は、固定的な
対応関係の制御又は２次の回帰モデルによる近似
カーブのピーク点による対応の決定手法に比較し
て制御性は格段に向上するが、なお第８図、第９
図に示すように、真の対応点ｊ点からのずれが発
生する場合がある。 このずれの原因は、例えば第８図の４次の近似
の場合では、両端の真ん中の周辺では、実測デー
タをよく追いかけているが、両端特に左端では生
データをあまり追いかけていない。ここで重要な
ことは、この両端、つまスライス操作の影響の小
さいところの変動が外乱による場合は、ここを忠
実に追いかけないほうが良い結果が得られる点で
ある。 第９図の６次の最小２乗近似カーブでは、両端
を４次の場合よりは忠実に追いかけているために
かえつてこれが災いして全体的に右にずれ、誤差
が＋１ポイント発生し、４次の場合と変化がない
結果となつている。 本発明は、このような両端のデータに関する問
題点を解消し、対応点決定の精度をより向上させ
たプロフイル測定点決定方法の提供を目的とする
ものである。 ＜問題点を解決するための手段＞ このような目的を達成する本発明は、ｉ番目の
スリイスボルトを操作した後絶乾坪量プロフイル
変化が収束した時点での絶乾坪量プロフイルデー
タを４次以上の重み付け最小２乗法により近似
し、この近似式のピーク点に対応する測定点ｊを
上記ｉ番目のスライスボルトの対応測定点とみな
す絶乾坪量プロフイル測定点決定方法において、
次のようにしたものである。 即ち、前記最小２乗近似の対称となる絶乾坪量
プロフイルデータを、当該スライスボルト操作の
影響を受ける中央ゾーンと、当該スライスボルト
操作の影響に比べて外乱の影響を大きく受ける左
右ゾーンに区分し、前記絶乾坪量プロフイルデー
タの重み付け最小２乗法における重み係数を、中
央ゾーンでの重み係数に比較して、左右ゾーンの
評価を小さくする重み係数としたことを特徴とし
ている。 ＜作用＞ 絶乾坪量プロフイルの収束値データが４次以上
の重み付けされた最小２乗法により近似され、こ
の近似式のピーク点に対する測定点ｊが操作され
たスライスボルトの対応測定点と見なされる。重
み付けは、中央ゾーンを左右ゾーンに比較して重
くしたので、外乱の影響を受ける左右ゾーンを小
さく評価してスライスボルト操作の影響のみを抽
出している。そこで、外乱の影響を受けず、正確
な対応点が得られる。 ＜実施例＞ 第１図は本発明の測定点決定方法を適用した装
置の実施例を示す構成図であり、第７図との相違
点は、データ処理装置１８におけるCD方向の距
離ｒのプロフイルの出力部２３と、このプロフイ
ル出力に対して重み付けを実行する重み決定部２
４が追加されており、他の構成要素並びに作用は
同一である。 一般にに関するｎ次の回帰モデルは、 N n(r)＝o 〓i=0 airi (1) である。これより、次式が得られる。 〓N n（rn′）ΛoФo (2) ここで、 Λo＝１、 １、 〓 １、r1、 r2、 ro′、r1 2、 r2 2、 ro′2、r1 3、 r2 3、 ro′3、…r1 n …r2 n 〓 …ro′n （n′≧ｎ） (3) ФT o＝［a0、a1、a2、……ao］ (4) また、実測データｂ（ro′）より、 ［〓N （ro′）］T＝［bN （r1）、bN （r2）… bN （ri）、bN （ri+1）…bN （rj-1）、 bN （rj）…bN （rj+1）］ (5) が得られる。 本発明では、重み関数Ｒを導入して、 Jo′＝（〓N （ro′）−〓N n（ro′））T ×R-1（〓N （ro′）−〓N n（ro′）） (6) Joを最小にするようにΦを決定する。 Фo＝［ΛT oR-1Λo］-1 ×ΛT oR-1 〓N （ro′） (7) ここで、重みＲは、第２図に示すようにｉ個の
右端ゾーン、ｊ個の左端ゾーンと中央ゾーンと
で、次のように重みを異ならせて設定する。 【表】 ここで、rR、rL≧１に選定され、rRの値は、bN
（r1）〜bN （ri）の最大値を｜bN (k)｜とすると、こ
れの例えば10倍に選択する。同様に、rLの値は、
bN（rj）〜bN （rj＋ｌ）中の最大値を｜bN (m)｜と
すると、これの例えば10倍に選択する。 さらに(8)式より、Ｒの逆行列は、 【表】 となるから、両端の影響は中央ゾーンの1/10とな
る。 さらに(7)式より、 bN（ri）＝n 〓i=0 airi (10) が得られ、これよりCD方向の実測データｂの指
定値（CD方向の実測データのフイルタ処理と同
等）が得られる。(7)式の計算は、更に演算時間の
短縮と演算精度の向上のために以下のような平行
根演算形式に修正される。 ΛT oR-1Λoは実対称行列ゆえ、平方根行列で表
すことができる。そこで、 Ｖ＝ΛT oR-1Λo (11) とすると、Ｖは次のように分解される。 Ｖ＝SST (12) ここでＳは、行列Ｖの平方根行列であり、下３
角行列となり、Ｓは上３角行列である。(12)式より
(7)式は、 Φo＝（S-TS-1）（ΛT oR-1 〓N （ro）） (13) となる。両辺にＳを乗ずると、 STΦo ＝S-1［ΛT oR-1 〓N （ro）］ (14) となる。 ここで、 Zo′＝S-1［Λo TR-1 〓N （ro）］ (15) とすると、(13)式は、 STΦo＝Zo′ (16) ここで、 ZT o′＝［Z1、Z2、……Zo′］ (17) STは、上３角行列故に、(15)式は、元の(7)式のよ
うに［ΛT oR-1Λo］-1の行列逆演算が不要であり、
単に代数演算処理ができるために、計算時間を著
しく短縮することができ、精度も向上する。次
に、 ST Φo＝S11、 ０、 〓 ０、S12、 S22、 ０、…S1o …S2o 〓 Sooa1 a2 〓 ao (18) から、 Sooao＝Zo′ (19) So-1o-1ao-1＋So-1oao ＝Zo-1′ (20) So-2o-2ao-2＋So-2o-1 ×ao-1＋So-2o-1ao ＝Zo-2′ （21） … S11a1＋S12a2＋…＋S1oao ＝Z1′ （22） (19)式より最初にaoを、 ao＝Zo／Soo （23） で求め、次に(20)式よりao-1を、 ao-1＝（Zo-1−So-1oao ／So-1o-1 （24） で求め、この順でai、……a3、a2、a1を次々に求
めて行けばよい。 以上の手順でΦoを求めるのは、第１図の２０
で実行され、bN (r)の計算は２１で実行されること
になる。第３図は、本発明による重み付けを実行
した場合（点線で示す）と重みをすべて１に設定
した場合（実践で示す）の比較結果を示すもので
あり、重み付けを実行しない場合の対応誤差＋１
ポイントが改善され、誤差がした場合のほうが対
応誤差がゼロとなり、精度の向上を期待すること
が可能である。 第４図は以上説明した最小２乗法ｎ次の計算手
順をフローチヤートで示したものである。 本発明では、最小２乗法の次数は、要求される
測定精度に応じて４以上の任意の偶数次数に選択
されるが、計算の繁雑さなどを考慮すれば、４次
又は６次程度が実用的である。 上記実施例では、データは平滑処理実測データ
ｂ(r)を用いたが、未平滑の実測データｂ(r)を用い
てもよい。 さらに、Ｖ＝ΛT oR-1Λoの平方根行列Ｓを求め
ずに、 Φo＝［ΛT oR-1Λo］-1 ×ΛoR-1 〓N （ro′） の計算から直接Φoを求めるようにすることも可
能である。 ＜発明の効果＞ 以上説明したように、本発明によれば、重みＲ
を幅方向プロフイル或いは幅方向各BD、Ｂ（r1）
〜ｂ（rn′）の外乱によるノイズの影響を、中央ゾ
ーンと左右ゾーンとの重み付けの差により除去す
ることで、スライス対応即定点を正確に決定でき
るという効果がある。 また、坪量制御装置においては、均一なプロフ
イルを得るためにスライス操作を行つているが、
この場合プロフイル制御の分配係数により制御を
行つている。ここで、本発明により得られる最小
２乗近似カーブから正しい干渉パターン（スライ
ス、BD間）が得られ、プロフイル制御の分配係
数が正確に決定できる。これにより、高精度の
BDプロフイール制御を行うことも可能になると
いう効果がある。
作端であるスライスボルトの位置に対応する絶乾
坪量プロフイル測定点の位置を精度よく決定する
方法に関する。 ＜従来技術＞ 第５図は、抄紙プロセスの概要図である。原料
パルプＰは種箱１、種口弁２を経て白水サイロ３
からのリターン原料と合流してポンプ４によりヘ
ツドボツクス５に供給される。ヘツドボツクス内
の原料は、スライスリツプ６の間隙よりワイヤー
パート８上にシート状に吐出され、プレスパート
９で搾水され、ドライパート１０で蒸気により加
熱乾燥された後、カレンダー１１を経てリール１
２に巻き取られて製品となる。７はスライスリツ
プの間隙開度を操作するための複数（例えば40
本）のスライスボルトである。 １３は製品を紙幅方向にスキヤンして紙の坪量
（１m2の重さ）と水分率を測定するＢ／Ｍ計であ
り、この測定値に基づいて紙の絶乾坪量プロフイ
ルが紙の横幅方向を複数点、例えば60点又は360
点に分割して計算される。 １４は、計算で求められた各点の絶乾坪量プロ
フイル測定値PVとプロフイル設定値SVの差を制
御演算する制御装置であり、その操作出力MVに
より測定点に対応するスライスボルト７が駆動さ
れ、スライスリツプ６の開度が操作され、製品の
絶乾坪量プロフイルが設定絶乾坪量プロフイルに
一致するようなフイードバツク制御が実行され
る。 ＜発明が解決しようとする問題点＞ 初期の絶乾坪量プロフイル制御システムでは、
ヘツドボツクス５よりスライスリツプ６を介して
吐出した原料は、第６図Ａに示す様にヘツドボツ
クスのｉ番スライスボルト位置よりまつすぐに飛
び出し、ワイヤー上をまつすぐ流れ、乾燥段階で
一様に収縮するとしてｊ番プロフイル測定点位置
に対応するものとし、Ｂに示すようにプロフイル
の測定値PVと設定値SVの偏差に基づいて対応す
るスライスボルトを操作して絶乾坪量プロフイル
の制御を実行している。Ｃはこのように制御され
た場合のスライスリツプ開度プロフイルの例であ
る。 所が、実際には、原料はヘツドボツクスからま
つすぐに吐出するとは限らないし、流れもまつす
ぐではない。又、収縮率も紙端近傍のほうが中央
部より大きく、現実の流れによる対応点は、ｊ番
目からずれてｌ番のプロフイル測定点となつてい
る場合が多い。 このように、実際の測定点と操作端の位置対応
のずれが大きいと、絶乾坪量プロフイル制御にお
いて好ましくないプロフイルの鋸歯状波を生ずる
ことがある。 出願人は、このような問題点を解消する対応点
決定手法として４次以上の回帰モデルによる最小
自乗法による対応点の決定手法を特願昭62−
264917号で提案している。 第７図によりその概要を説明する。ｉ番スライ
スボルトを操作したときに、抄紙機制御装置１４
からは例えば60点の幅方向のプロフイル測定点の
平滑化されたＮスキヤン分の複数の時系列デー
タ、1 1〜N 1、1 2〜N 2…1 60〜N 60がオンライ
ン
出力され、Ｉ／Ｏコントローラ１５を介してデー
タ記録装置１６にBD（絶乾坪量）データとして
フロツピーデイスクなどの記憶手段に記録保持さ
れる。 このBDデータは、オンラインでＩ／Ｏコント
ローラ１７を介してデータ処理装置１８のBDデ
ータ読出し部１９により読み出され、定常状態と
見なされるＮスキヤンのCD方向（幅方向）のBD
データの収束値N 1、N 2…N 60を20で４次又は６
次の回帰モデル、 N n(r)＝4 〓i=0 airi または、 N n(r)＝6 〓i=0 airi により最小２乗処理する。ここでｒはCO方向の
距離、ｂはBD出力の平滑化された値である。 第８図及び第９図は、この処理結果得られた最
小２乗４次及び６次の近似カーブであり、そのピ
ーク値位置がピーク値決定部２１により求めら
れ、この位置（ｊ＋１）番をｉ番スライスボルト
の対応測定点として出力部２２より出力し、オフ
ライン的に抄紙機制御装置１４の対応点データを
修正する。 ＜発明が解決しようとする問題点＞ このような最小２乗法４次又は６次の回帰モデ
ルによるピーク値位置を求める手法は、固定的な
対応関係の制御又は２次の回帰モデルによる近似
カーブのピーク点による対応の決定手法に比較し
て制御性は格段に向上するが、なお第８図、第９
図に示すように、真の対応点ｊ点からのずれが発
生する場合がある。 このずれの原因は、例えば第８図の４次の近似
の場合では、両端の真ん中の周辺では、実測デー
タをよく追いかけているが、両端特に左端では生
データをあまり追いかけていない。ここで重要な
ことは、この両端、つまスライス操作の影響の小
さいところの変動が外乱による場合は、ここを忠
実に追いかけないほうが良い結果が得られる点で
ある。 第９図の６次の最小２乗近似カーブでは、両端
を４次の場合よりは忠実に追いかけているために
かえつてこれが災いして全体的に右にずれ、誤差
が＋１ポイント発生し、４次の場合と変化がない
結果となつている。 本発明は、このような両端のデータに関する問
題点を解消し、対応点決定の精度をより向上させ
たプロフイル測定点決定方法の提供を目的とする
ものである。 ＜問題点を解決するための手段＞ このような目的を達成する本発明は、ｉ番目の
スリイスボルトを操作した後絶乾坪量プロフイル
変化が収束した時点での絶乾坪量プロフイルデー
タを４次以上の重み付け最小２乗法により近似
し、この近似式のピーク点に対応する測定点ｊを
上記ｉ番目のスライスボルトの対応測定点とみな
す絶乾坪量プロフイル測定点決定方法において、
次のようにしたものである。 即ち、前記最小２乗近似の対称となる絶乾坪量
プロフイルデータを、当該スライスボルト操作の
影響を受ける中央ゾーンと、当該スライスボルト
操作の影響に比べて外乱の影響を大きく受ける左
右ゾーンに区分し、前記絶乾坪量プロフイルデー
タの重み付け最小２乗法における重み係数を、中
央ゾーンでの重み係数に比較して、左右ゾーンの
評価を小さくする重み係数としたことを特徴とし
ている。 ＜作用＞ 絶乾坪量プロフイルの収束値データが４次以上
の重み付けされた最小２乗法により近似され、こ
の近似式のピーク点に対する測定点ｊが操作され
たスライスボルトの対応測定点と見なされる。重
み付けは、中央ゾーンを左右ゾーンに比較して重
くしたので、外乱の影響を受ける左右ゾーンを小
さく評価してスライスボルト操作の影響のみを抽
出している。そこで、外乱の影響を受けず、正確
な対応点が得られる。 ＜実施例＞ 第１図は本発明の測定点決定方法を適用した装
置の実施例を示す構成図であり、第７図との相違
点は、データ処理装置１８におけるCD方向の距
離ｒのプロフイルの出力部２３と、このプロフイ
ル出力に対して重み付けを実行する重み決定部２
４が追加されており、他の構成要素並びに作用は
同一である。 一般にに関するｎ次の回帰モデルは、 N n(r)＝o 〓i=0 airi (1) である。これより、次式が得られる。 〓N n（rn′）ΛoФo (2) ここで、 Λo＝１、 １、 〓 １、r1、 r2、 ro′、r1 2、 r2 2、 ro′2、r1 3、 r2 3、 ro′3、…r1 n …r2 n 〓 …ro′n （n′≧ｎ） (3) ФT o＝［a0、a1、a2、……ao］ (4) また、実測データｂ（ro′）より、 ［〓N （ro′）］T＝［bN （r1）、bN （r2）… bN （ri）、bN （ri+1）…bN （rj-1）、 bN （rj）…bN （rj+1）］ (5) が得られる。 本発明では、重み関数Ｒを導入して、 Jo′＝（〓N （ro′）−〓N n（ro′））T ×R-1（〓N （ro′）−〓N n（ro′）） (6) Joを最小にするようにΦを決定する。 Фo＝［ΛT oR-1Λo］-1 ×ΛT oR-1 〓N （ro′） (7) ここで、重みＲは、第２図に示すようにｉ個の
右端ゾーン、ｊ個の左端ゾーンと中央ゾーンと
で、次のように重みを異ならせて設定する。 【表】 ここで、rR、rL≧１に選定され、rRの値は、bN
（r1）〜bN （ri）の最大値を｜bN (k)｜とすると、こ
れの例えば10倍に選択する。同様に、rLの値は、
bN（rj）〜bN （rj＋ｌ）中の最大値を｜bN (m)｜と
すると、これの例えば10倍に選択する。 さらに(8)式より、Ｒの逆行列は、 【表】 となるから、両端の影響は中央ゾーンの1/10とな
る。 さらに(7)式より、 bN（ri）＝n 〓i=0 airi (10) が得られ、これよりCD方向の実測データｂの指
定値（CD方向の実測データのフイルタ処理と同
等）が得られる。(7)式の計算は、更に演算時間の
短縮と演算精度の向上のために以下のような平行
根演算形式に修正される。 ΛT oR-1Λoは実対称行列ゆえ、平方根行列で表
すことができる。そこで、 Ｖ＝ΛT oR-1Λo (11) とすると、Ｖは次のように分解される。 Ｖ＝SST (12) ここでＳは、行列Ｖの平方根行列であり、下３
角行列となり、Ｓは上３角行列である。(12)式より
(7)式は、 Φo＝（S-TS-1）（ΛT oR-1 〓N （ro）） (13) となる。両辺にＳを乗ずると、 STΦo ＝S-1［ΛT oR-1 〓N （ro）］ (14) となる。 ここで、 Zo′＝S-1［Λo TR-1 〓N （ro）］ (15) とすると、(13)式は、 STΦo＝Zo′ (16) ここで、 ZT o′＝［Z1、Z2、……Zo′］ (17) STは、上３角行列故に、(15)式は、元の(7)式のよ
うに［ΛT oR-1Λo］-1の行列逆演算が不要であり、
単に代数演算処理ができるために、計算時間を著
しく短縮することができ、精度も向上する。次
に、 ST Φo＝S11、 ０、 〓 ０、S12、 S22、 ０、…S1o …S2o 〓 Sooa1 a2 〓 ao (18) から、 Sooao＝Zo′ (19) So-1o-1ao-1＋So-1oao ＝Zo-1′ (20) So-2o-2ao-2＋So-2o-1 ×ao-1＋So-2o-1ao ＝Zo-2′ （21） … S11a1＋S12a2＋…＋S1oao ＝Z1′ （22） (19)式より最初にaoを、 ao＝Zo／Soo （23） で求め、次に(20)式よりao-1を、 ao-1＝（Zo-1−So-1oao ／So-1o-1 （24） で求め、この順でai、……a3、a2、a1を次々に求
めて行けばよい。 以上の手順でΦoを求めるのは、第１図の２０
で実行され、bN (r)の計算は２１で実行されること
になる。第３図は、本発明による重み付けを実行
した場合（点線で示す）と重みをすべて１に設定
した場合（実践で示す）の比較結果を示すもので
あり、重み付けを実行しない場合の対応誤差＋１
ポイントが改善され、誤差がした場合のほうが対
応誤差がゼロとなり、精度の向上を期待すること
が可能である。 第４図は以上説明した最小２乗法ｎ次の計算手
順をフローチヤートで示したものである。 本発明では、最小２乗法の次数は、要求される
測定精度に応じて４以上の任意の偶数次数に選択
されるが、計算の繁雑さなどを考慮すれば、４次
又は６次程度が実用的である。 上記実施例では、データは平滑処理実測データ
ｂ(r)を用いたが、未平滑の実測データｂ(r)を用い
てもよい。 さらに、Ｖ＝ΛT oR-1Λoの平方根行列Ｓを求め
ずに、 Φo＝［ΛT oR-1Λo］-1 ×ΛoR-1 〓N （ro′） の計算から直接Φoを求めるようにすることも可
能である。 ＜発明の効果＞ 以上説明したように、本発明によれば、重みＲ
を幅方向プロフイル或いは幅方向各BD、Ｂ（r1）
〜ｂ（rn′）の外乱によるノイズの影響を、中央ゾ
ーンと左右ゾーンとの重み付けの差により除去す
ることで、スライス対応即定点を正確に決定でき
るという効果がある。 また、坪量制御装置においては、均一なプロフ
イルを得るためにスライス操作を行つているが、
この場合プロフイル制御の分配係数により制御を
行つている。ここで、本発明により得られる最小
２乗近似カーブから正しい干渉パターン（スライ
ス、BD間）が得られ、プロフイル制御の分配係
数が正確に決定できる。これにより、高精度の
BDプロフイール制御を行うことも可能になると
いう効果がある。
第１図は本発明方法を適用した測定装置の実施
例を示す構成図、第２図は重み付け設定の説明
図、第３図は本発明方法による重み付けをした場
合と重み付けをしない場合の精度の比較説明図、
第４図はｎ次の重み付け最小２乗法の計算手順を
示すフローチヤート図、第５図は抄紙機における
プロフイル制御の概要を説明する構成図、第６図
はスライスボルトと測定点の対応関係に関する説
明図、第７図は従来技術の一例を示す構成図、第
８図、第９図は従来技術による対応点決定の説明
図である。 １４……抄紙機制御装置、１５，１７……Ｉ／
Ｏコントローラ、１６……データ記録装置、１８
……データ処理装置、１９……データ読出し部、
２０……最小２乗４次、６次処理部、２１……ピ
ーク値決定部、２２……測定点出力部、２３……
プロフイル出力部、２４……重み決定部。
例を示す構成図、第２図は重み付け設定の説明
図、第３図は本発明方法による重み付けをした場
合と重み付けをしない場合の精度の比較説明図、
第４図はｎ次の重み付け最小２乗法の計算手順を
示すフローチヤート図、第５図は抄紙機における
プロフイル制御の概要を説明する構成図、第６図
はスライスボルトと測定点の対応関係に関する説
明図、第７図は従来技術の一例を示す構成図、第
８図、第９図は従来技術による対応点決定の説明
図である。 １４……抄紙機制御装置、１５，１７……Ｉ／
Ｏコントローラ、１６……データ記録装置、１８
……データ処理装置、１９……データ読出し部、
２０……最小２乗４次、６次処理部、２１……ピ
ーク値決定部、２２……測定点出力部、２３……
プロフイル出力部、２４……重み決定部。
Claims (1)
- 【特許請求の範囲】 １ ｉ番目のスライスボルトを操作した後絶乾坪
量プロフイル変化が収束した時点での絶乾坪量プ
ロフイルデータを４次以上の重み付け最小２乗法
により近似し、 この近似式のピーク点に対応する測定点ｊを上
記ｉ番目をスライスボルトの対応測定点とみなす
絶乾坪量プロフイル測定点決定方法において、 前記最小２乗近似の対称となる絶乾坪量プロフ
イルデータを、当該スライスボルト操作の影響を
受ける中央ゾーンと、当該スライスボルト操作の
影響に比べて外乱の影響を大きく受ける左右ゾー
ンに区分し、 前記絶乾坪量プロフイルデータの重み付け最小
２乗法における重み係数を、中央ゾーンでの重み
係数に比較して、左右ゾーンの評価を小さくする
重み係数としたことを特徴とする絶乾坪量プロフ
イル測定点決定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP26309387A JPH01111087A (ja) | 1987-10-19 | 1987-10-19 | 絶乾坪量プロフィル測定点決定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP26309387A JPH01111087A (ja) | 1987-10-19 | 1987-10-19 | 絶乾坪量プロフィル測定点決定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH01111087A JPH01111087A (ja) | 1989-04-27 |
| JPH0555637B2 true JPH0555637B2 (ja) | 1993-08-17 |
Family
ID=17384732
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP26309387A Granted JPH01111087A (ja) | 1987-10-19 | 1987-10-19 | 絶乾坪量プロフィル測定点決定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH01111087A (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109215019B (zh) * | 2018-08-24 | 2021-04-23 | 华南农业大学 | 一种基于原木弯曲度的木材切分点定位方法和装置 |
| US11673420B2 (en) | 2020-11-30 | 2023-06-13 | Jessem Tool Corporation | Layout tool with grips |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6026877A (ja) * | 1983-07-20 | 1985-02-09 | 住友金属工業株式会社 | 油井管用管継手 |
| JPH0159394U (ja) * | 1987-10-08 | 1989-04-13 |
-
1987
- 1987-10-19 JP JP26309387A patent/JPH01111087A/ja active Granted
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6026877A (ja) * | 1983-07-20 | 1985-02-09 | 住友金属工業株式会社 | 油井管用管継手 |
| JPH0159394U (ja) * | 1987-10-08 | 1989-04-13 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH01111087A (ja) | 1989-04-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US3767900A (en) | Adaptive controller having optimal filtering | |
| US4374703A (en) | Control system for papermaking machine headbox | |
| US3711688A (en) | Computer control of paper machine in which basis weight is controlled through control of stock flow | |
| WO1999004090A1 (en) | Method and apparatus for executing grade change in paper machine | |
| US3676295A (en) | Noninteracting control of moisture and fiber content of fibrous sheet during manufacture | |
| KR101100660B1 (ko) | 시뮬레이션 방법, 섬유 배향 제어 방법, 및 섬유 배향 제어장치 | |
| US5781440A (en) | Process and device for monitoring the cross-sectional profile of a continuously produced web of material | |
| US4939929A (en) | Measurement of the properties of a web in paper production | |
| JPH0555637B2 (ja) | ||
| US3586601A (en) | Basis weight control system for a paper making machine | |
| JPH0547674B2 (ja) | ||
| EP1259674B1 (en) | Method and apparatus for identifying mapping of paper machine actuator | |
| JP2006526714A (ja) | パーシャル最小二乗法に基づく紙のカール及び捩れのモデル化、予測及びコントロール | |
| JPS583075B2 (ja) | シヨウシキニオケルプロフアイルコントロ−ルノホウホウ | |
| JPH0547672B2 (ja) | ||
| JP2512318B2 (ja) | 坪量プロフィル測定方法 | |
| JP2884696B2 (ja) | 抄紙機の抄速変更・抄替制御装置及びその方法 | |
| JPS63227885A (ja) | 絶乾坪量プロフイル測定点決定方法 | |
| JPH0159394B2 (ja) | ||
| JPH063000B2 (ja) | プロフィル制御装置 | |
| JPS63182491A (ja) | 厚さのプロフイ−ル制御方法 | |
| JPH0583677B2 (ja) | ||
| JPS629405A (ja) | プロセス制御装置 | |
| JPS6343513B2 (ja) | ||
| JPS6071793A (ja) | 抄紙機のグレ−ドチエンジ制御装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| LAPS | Cancellation because of no payment of annual fees |