Как узнать, какое линейное соответствие лучше?

Допустим, у нас есть простой гармонический осциллятор масса/струна. Измерив период движения, я могу вычислить угловую частоту. Знаю это

ю "=" к м и Т "=" 2 π ю

мы можем прийти к двум функциям:

ю 2 "=" 1 м к и м "=" Т 2 4 π 2 к

которая в основном представляет собой ту же функцию, выражающую ю по-другому. Оба соотношения являются линейными с к .

Чтобы рассчитать к , я использовал линейную подгонку. Я сделал два сюжета, с ю 2 на у ось и 1 / м на Икс ось на первом графике и м по оси Y и Т 2 / 4 π 2 на Икс ось на втором графике. Но эти два графика дали мне совершенно разные результаты, с разницей почти в 50%. Как узнать, какой из них правильный и почему это произошло?

Используемые данные:

T (s)  0.8283 0.9622 1.0912 1.1195 1.2896
m (kg) 0.02   0.03   0.04   0.05.  0.06
Ваши необработанные данные измеряются с точки зрения Т и м , поэтому можно было бы подумать, что ваша линейная подгонка будет также с точки зрения этих переменных. Почему вы хотите «массировать» данные, чтобы построить график с точки зрения ю 2 и 1 м ?

Ответы (2)

Проблема в том, что ваши данные плохо моделируются прямой пропорциональностью между м и Т 2 / 4 π 2 . Обе посадки должны проходить через начало координат; но если вы наивно построите линейную подгонку по методу наименьших квадратов в каждом случае, вы обнаружите существенное смещение для ваших данных:

введите описание изображения здесь

Вы можете получить лучшее согласие между двумя методами, заставив кривую соответствия проходить через начало координат (детали того, как это сделать, различаются в зависимости от используемого вами программного обеспечения). Согласие все еще невелико, но оно лучше.

В более широком смысле, не обязательно, что наклон линии, которая минимизирует расстояние по методу наименьших квадратов для функции у "=" м Икс будет таким же, как наклон линии, которая минимизирует расстояние по методу наименьших квадратов для функции ( 1 / Икс ) "=" м ( 1 / у ) . Это связано с тем, что «наивная» подгонка методом наименьших квадратов предполагает равномерно распределенные ошибки во всех ваших точках данных и незначительные ошибки в независимой переменной; но когда вы меняете местами два графика, вы меняете роли зависимой и независимой переменных и меняете размеры ошибок на них. В общем, наилучшей практикой на вводном занятии было бы взять в качестве независимой переменной ту, которая измеряется более точно (в данном случае, вероятно, массу), использовать распространение ошибки для оценки планок погрешностей для каждого из ваших измерений зависимую переменную и выполнить подгонку данных с соответствующим взвешиванием.

Это не линеаризация, это просто выражение к двумя разными способами. Если вы выберете два случайных значения для Т и м вы получите два разных результата. Подвох в том, что Т и м связаны через к , т.е. к "=" ф ( Т , м ) или к "=" г ( ю , м ) и в обоих случаях вы должны получить тот же результат, если вы удовлетворяете ю "=" 2 π Т и вы используете то же самое м .

ф ( Т , м ) "=" 4 π 2 Т 2 м , г ( ю , м ) "=" ю 2 м

В общем, если к 0 "=" ф ( Т 0 , м 0 ) "=" г ( ю 0 , м 0 ) где ю 0 "=" 2 π Т 0 затем

ф ( Т 0 , м 0 + Δ ) "=" г ( ю 0 , м 0 + Δ )

ф ( Т 0 + Δ , м 0 ) г ( ю 0 + Δ , м 0 )

поскольку функции ф и г линейны в м для фиксированного Т 0 и ю 0 , и нелинейный по Т и ю для фиксированного м 0 .


Что мне кажется, так это то, что вы запутались, как появляются сюжеты ( Т , м , ф ) и ( ю , м , г ) не похожи друг на друга. Чтобы понять это, рассмотрим более простой случай:

час ( Икс ) "=" Икс , г ( у ) "=" 1 у , Икс "=" 1 у

Очевидно, что ( Икс , час ) и ( у , г ) графики не совпадают, хотя две функции одинаковы в Икс . Однако, ( Икс , час ) и ( 1 / у , г ) сюжеты одинаковые, и это именно то, что происходит в вашем случае.

конечно, с каждой точкой сам по себе я получаю тот же результат, но наклон линии тренда другой. Разве это значение не должно быть равным k и таким же?
Я не уверен, что понимаю, что вы имеете в виду. Если вы спросите, следует ли к разница одинакова при равных приращениях ю и ответ нет, так как ю и Т обратно пропорциональны.
Что я сделал, так это построил (две функции отдельно) вычисленные точки и вычислил наклон линии тренда. Я понимаю, что k должен быть наклоном, но результаты разные.
Это потому, что две функции для к нелинейны в Т и ю .
Чтобы быть более конкретным, график представляет собой ω2 по оси y и 1/m по оси x на первом графике, и m по оси y и T2/4π2 по оси x на втором графике.
Да, в том-то и дело! Если вы положите 4 π 2 / Т 2 вместо Т 2 / 4 π 2 вы получите то, что ожидаете, т.е. тот же наклон.
Прости, что я не понимаю. Отношение между y и x на графике является линейным с наклоном k. Я не рисую ни T, ни ω сами по себе.
Давайте упростим это, представьте, что у вас есть ф ( Икс ) "=" Икс и г ( у ) "=" у где у "=" 1 / Икс . Будут ли сюжеты ( Икс , ф ) и ( у , г ) быть таким же?
Нет, графики разные, но разве наклон обеих функций не должен быть равен 1, хотя точки разные?
Если графики разные, почему вы ожидаете, что наклон будет одинаковым? Однако в упрощенном примере графики ( Икс , ф ) и ( 1 / у , г ) одинаковы, что также означает один и тот же наклон.
Итак, в данном случае, какой наклон равен k? Я все еще немного в замешательстве..
Можете ли вы быть более конкретным, что вы подразумеваете под «какой склон к "? Вы имеете в виду, каков наклон ф ( Т , м ) в отношении Т и наклон г ( ю , м ) в отношении ю ?
Может поможет, если вы покажете графики.