ULprojectのホームページ © ULproject

 2022/3/25(金) N88-BASICで点と直線の距離 (1回目)   初めに 1. 原点(0,0)と直線ax+by+c=0との距離dを求める     y切片と傾きを書いたグラフの 2つの相似な直角三角形の 辺の比で求めています   直線ax + by + c = 0 [ y = -(a/b)x - c/b ]と 原点O(0, 0)との距離dは 上記図より d = |c| / √(a 2 +b 2 ) となる   次にグラフの平行移動を考える 2. 平行移動   ちなみに、平行移動を使って A(x 1 ,y 1 )を通る傾きmの直線の方程式 y-y 1 =m(x-x 1 ) を求める事ができる   最後に 3. 点P(x 1 ,y 1 )と直線ax+by+c=0との距離dを求める     点P(x 1 ,y 1 )と直線ax+by+c=0を x方向に-x 1 、y方向に-y 1 移動させると 点Pは原点(0,0)に、直線はax+by+(ax 1 +by 1 +c)=0 に移動するが、距離dは変わらないので 原点と直線の距離の公式のcを(ax 1 +by 1 +c)で 置き換えれば良い よって 直線L: ax + by + c = 0 [ y = -(a/b)x - c/b ]と 点P(x 1 , y 1 )との距離d = |ax 1 +by 1 +c|/√(a 2 +b 2 ) となる   また、距離d = |ax 1 +by 1 +c|/√(a 2 +b 2 )は 以下の図で求める事も出来る     点P(x 1 , y 1 )を-x 1 , -y 1  移動して原点に持ってくる ax + by + c = 0も-x 1 , -y 1  移動して a{x-(-x 1 )} + b{y-(-y 1 )} + c = 0 ax + by + ax 1 +by 1 +c = 0 この式と原点の距離がdと同じになる   上図のようにax + by + ax 1 +by 1 +c = 0 を変形して傾きと切片からグラフを書くと 相似な2つの直角三角形が書ける   三平方の定理より r' = √(a 2 /b 2  + 1 2 ) = √{(a 2 +b 2 )/b 2 } = √(a 2 +b 2 )/b   d : 1 = r : r' = (ax

2022/3/20(日) N88-BASICで同じ誕生日の確率   n人の誕生日が重複しない確率pを求める   (1年は365日で誕生日にバラツキがないとする) (誕生日が重複する確率は1-p)   1人目からn人目まで並べて 2人目が1人目と違う確率は(365-1)/365 3人目が1～2人目と違う確率は(365-2)/365 n人目が1～n-1人目と違う確率は{365-(n-1)}/365 なので、 n人の誕生日が重複しない確率pは p = (365-1)/365･(365-2)/365･ ... ･{365-(n-1)}/365 となります   別の求め方、   階乗(Factorial)、順列(Permutation)、組合せ(Combination) については https://ulprojectmail.blogspot.com/2021/11/n88-basicpc.html N88-BASICで順列組合せ   を参照して下さい   n人の誕生日が重複しない確率pを求める 別解   365種類の日付からn個を選び並べる方法が 365 P n  通りある つまり、n人が重複しない誕生日の選び方が 365 P n  通りある また、n人の誕生日の選び方は 365 n  通りあるので、   n人の誕生日が重複しない確率pは 365 P n  / 365 n  となる   p = 365 P n  / 365 n    = 365･364･ ... ･(365-n+1) / 365 n    = 364･363･ ... ･{365-(n-1)} / 365 n-1    = (365-1)/365･(365-2)/365･ ... ･{365-(n-1)}/365 と変形すると、前者の求め方と一致する   以上、2通りの考え方の式を導きました   プログラムでは前者の式を使用していますが、 計算精度の関係で確率は近似値になります (すごい誤差が出ているかもしれません)   また、10組分のシミュレートも表示します   同じ誕生日が2人以上の時 誕生日(人数)を表示しています   NL-BASICとblg～.zip(birth001.bas)は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Read

2022/3/18(金) N88-BASICで平均の落とし穴 (2回目)   1回目は平均(mean)が現実を表していない 例を示しました   2回目は中央値(median)が現実を表していない 例です ちなみに   平均は全員の合計年収を、人数で割った値   中央値は小さい順(大きい順でも可)に年収を 並べて真ん中の人が1人の時(全体が奇数人の時)は その人の年収になります 真ん中が2人の時(全体が偶数人の時)は その2人の平均年収になります   NL-BASICとblg～.zip(ave002.bas)は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい

2022/3/16(水) N88-BASICで平均の落とし穴 (1回目)   平均値を聞いた時、現実と違いすぎる と感じたことはありませんか   平均的なもの(average)を表す数値には 平均値(mean)、中央値(median)などがあります   それぞれ、数値によって、現実的な間隔と 近い場合と、遠い場合があります その様な例を表示しました   1回目は 平均年収が1000万円という村(人口10人)があり 高収入な人が多いと思い、村人に話を聞くと 生活が苦しいと言う人が多いと感じました   そんな事があるのでしょうか そんな村の年収を表示しました ちなみに   平均は全員の合計年収を、人数で割った値   中央値は小さい順(大きい順でも可)に年収を 並べて真ん中の人が1人の時(全体が奇数人の時)は その人の年収になります 真ん中が2人の時(全体が偶数人の時)は その 2人の平均年収になります NL-BASICとblg～.zip(ave001.bas)は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい    

2022/3/12(土) N88-BASICでCIRCLE (3回目)   前回の3点を通る円弧 を連続円弧 にします   入力した座標を記録し、 再生するプログラムを セーブ出来る様にしました   このプログラムは一筆書きのみの対応です   NL-BASICとblg～.zip(cir003.bas)は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

2022/3/10(木) N88-BASICでCIRCLE  (2回目)   前回の3点を通る円を円弧にし ます 3点と中心の位置ベクトルをA,B,C,O とすると、 ベクトルA-Oとx軸(正方向)のなす角θ を求めます   A-O=(x,y)とすると、 0°≦θ＜ 360°は、 -90°≦ Tan -1 (y/x) ≦ 90°として、 x ＝ 0 のとき、   不明  (if y ＝ 0)    90° (if y ＞ 0)   270° (if y ＜ 0) x ＞ 0 のとき、   Tan -1 (y/x)        (if y ≧ 0)   360°+ Tan -1 (y/x) (if y ＜ 0) x ＜ 0 のとき、   180°+ Tan -1 (y/x)   3点A,B,Cのどれを始点、間の点、終点 にするかによって、描かれる円弧が 変わります   BASICでは角度はラジアン(弧度法)なので 0≦θ＜2π です   今回は、A,B,Cの順に点を打つとして、 AからCへBを経由する円弧を描きます   A<B<C → A to C B<C<A → A to C C<A<B → A to C A<C<B → C to A B<A<C → C to A C<B<A → C to A   860 IF W > 0 THEN T = ATN(H / W) - (H < 0) * 2 * PI: RETURN この行の(H < 0)は、真なら-1、偽なら0になる   NL-BASICとblg～.zip(cir002.bas)は 以下のリンクからダウンロードできます https://ulprojectmail.blogspot.com Readme.txtを読んで遊んで下さい  

2022/3/8(火) N88-BASICでCIRCLE (1回目)   3点を通る円をCIRCLE命令で描画します   ベクトルは大文字、スカラーは小文字で 書くことにします   ベクトルの内積の行列表記 A=(a x ,a y ), B=(b x ,b y ) A･B = a x b x  + a y b y   |a x  a y ||b x | = a x b x  + a x b y  = A･B              |b y |   行列 M =|a b|の行列式det|M|=ad-bc              |c d| det|M|はベクトル(a b),(c d)または (a c),(b d)が作る平行四辺形の面積(方向有)   行列や、平行四辺形の面積は、このブログの   https://ulprojectmail.blogspot.com/2021/06/n88-basicmatrix-1.html N88-BASICで行列(matrix) (1回目)   https://ulprojectmail.blogspot.com/2022/03/n88-basicvector-2.html N88-BASICでベクトル (2回目)   で、それぞれ解説しています   逆行列 M -1 =   1      | d -b|  (det|M|= ad-bc ≠0)                      ad-bc |-c  a| 単位行列 I=|1 0| 　　　      |0 1| (ad-bc)M -1 M=|ad-bc -ab+ab|=|ad-bc   0  |                         |cd-cd -cb+ad|  |  0   ad-bc| M -1 M=|1 0| = I             |0 1|   -------------------------------------------- 円周上の 3点と中心の位置ベクトルを A,B,C,O、円の半径をrとする   A=(a x ,a y ), B=(b x ,b y ), C=(c x ,c y ), O=(x,y), r   中心から各 3点までの距離はrなので |A-O|=r |B-O|=r |C-O|=r   両辺 2乗して |A| 2