设为首页收藏本站

爱吱声

 找回密码
 注册
搜索
查看: 2553|回复: 4
打印 上一主题 下一主题

[信息技术] 继续请教问题:关于 Pytorch 的 Autograd

[复制链接]
  • TA的每日心情

    2025-9-22 22:19
  • 签到天数: 1183 天

    [LV.10]大乘

    跳转到指定楼层
    楼主
     楼主| 发表于 2023-2-14 13:09:28 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
    本帖最后由 雷达 于 2023-2-14 13:12 编辑
    8 k; _" d' N3 @- ?1 P9 B9 Q# H
    % @, D' b  E9 Y+ ~6 e" F* y为预防老年痴呆,时不时学点新东东玩一玩。0 a! d7 l2 o. }- E
    Pytorch 下面的代码做最简单的一元线性回归:6 W  d" g5 ^: e5 v) f, H
    ----------------------------------------------
    ; e& b% j! f: Simport torch
    : k- {: D0 z1 R: \8 r; o; G: r1 }" _import numpy as np
    7 s- }* T$ e; c$ S4 i! Uimport matplotlib.pyplot as plt
    + {' E' q0 S- ?import random
    & Z; m3 Z; a$ y% z& o
    ( \/ }0 e$ g  w8 i0 m/ Ux = torch.tensor(np.arange(1,100,1))
    & ^# A9 j9 L7 e: Hy = (x*27+15+random.randint(-2,3)).reshape(-1)  # y=wx+b, 真实的w0 =27, b0=15, P' b6 C* _/ i8 L6 N$ b0 h
    ! f+ e$ W+ F/ F
    w = torch.tensor(0.,requires_grad=True)  #设置随机初始 w,b
    , K0 W; n4 o$ p7 O" W! ?) q; h  Gb = torch.tensor(0.,requires_grad=True)
    . ]) `; o0 |% v2 n& r2 \! D, j  t
    ! ?  [: ^6 l" E) Pepochs = 100; M6 y. [9 o( G) a( r
    % d9 }$ a; R( n* q
    losses = []
    - \  B" _1 u6 Hfor i in range(epochs):2 ]3 V4 d7 P7 r* n$ M
      y_pred = (x*w+b)    # 预测' s- [& [; ^* Z
      y_pred.reshape(-1)
    / v0 N! h1 Z8 O( J" g/ ]+ b# i  X
    % w+ U! d& t8 `& B! R2 Q  loss = torch.square(y_pred - y).mean()   #计算 loss
    , H; ]3 ]! ?3 X& S" R# p  losses.append(loss)
      M* \' [! o, R% k  
    2 Y5 P! ], u8 G7 J  loss.backward() # autograd
    ( V) S7 A! q1 q# V. A3 ?0 ^  with torch.no_grad():
    / r( U' s. h3 V, M    w  -= w.grad*0.0001   # 回归 w
    3 k- c0 K3 `# B9 G8 O# l    b  -= b.grad*0.0001    # 回归 b % v* v5 j* ?3 j) a
      w.grad.zero_()  . ?! B; R+ `( S; _; Y
      b.grad.zero_()
    8 p9 s: y# M3 Y9 I+ |. E0 a
    3 C8 V/ r. ~* `print(w.item(),b.item()) #结果9 S5 A6 j/ t) Z/ h; r/ x4 @7 I9 D3 ?
    ) V, h$ S, R7 V8 _; g0 _: w
    Output: 27.26387596130371  0.49745178222656251 |+ B" S! {! f5 c- c# g
    ----------------------------------------------& Q) a& L( ]+ l: U7 s+ g, G
    最后的结果,w可以回到 w0 = 27 附近,b却回不去 b0=15。两处红字,损失函数是矢量计算后的均值,感觉 b 的回归表达有问题。0 ^6 s! f, D  N3 I
    高手们帮看看是神马原因?- k  \0 M" ~; D3 K! B# I

    评分

    参与人数 1爱元 +10 收起 理由
    老票 + 10 不明觉厉

    查看全部评分

    该用户从未签到

    沙发
    发表于 2023-2-14 19:23:02 | 只看该作者
    本帖最后由 老福 于 2023-2-14 21:58 编辑
    ( f0 b; V% e7 ~% ^, M7 X% y  ~* U1 [. E4 n" m
    没有用过pytorch,但你把随机噪音部分改成均值为0的正态分布再试试看是不是符合预期?: L) D( r$ D# j$ v
    -------
    7 H3 i$ |, P$ `+ Y2 q不好意思,再看一遍,好像你在自算回归而不是用现成的工具直接出结果,上面的评论只有一点用,就是确认是不是算法有问题。/ K6 D! e6 l0 T$ @, y6 |; f, F
    -------
    7 j8 x( V' `& u8 A8 k算法诊断部分,建议把循环次数改为1000, 再看看loss是不是收敛。有点怀疑你循环次数不够,因为你起点是0, 步长很小。只是直观建议。

    评分

    参与人数 1爱元 +10 收起 理由
    雷达 + 10 谢谢建议

    查看全部评分

    回复 支持 反对

    使用道具 举报

  • TA的每日心情

    2025-9-22 22:19
  • 签到天数: 1183 天

    [LV.10]大乘

    板凳
     楼主| 发表于 2023-2-14 21:52:57 | 只看该作者
    老福 发表于 2023-2-14 19:23
    # I5 r+ _6 e" H( K2 g4 W没有用过pytorch,但你把随机噪音部分改成均值为0的正态分布再试试看是不是符合预期?
    * m0 h% L! B, H; L! z2 \-------
    5 O# t4 F- N7 l$ n( ?不好意思, ...
    2 z1 z7 t! H- G" l% f4 [
    谢谢,算法应该没问题,就是最简单的线性回归。
    4 O" z  I- W' ?& L% h我特意没有用现成的工具,就是想从最基本的地方深入理解一下。
    回复 支持 反对

    使用道具 举报

    该用户从未签到

    地板
    发表于 2023-2-14 22:00:48 | 只看该作者
    本帖最后由 老福 于 2023-2-14 22:02 编辑
    9 ^8 z, v& A, u$ v, J9 w5 Z
    雷达 发表于 2023-2-14 21:524 K: M( Z$ C  V/ Q! J9 @5 I
    谢谢,算法应该没问题,就是最简单的线性回归。
    % X) _' @8 E! F! P' O- C我特意没有用现成的工具,就是想从最基本的地方深入理解 ...

    4 e/ B. [0 v$ N& C; w
    " i; y- G2 [3 t' ]% T5 k; g& ^刚才更新了一下,建议增加循环次数或调一下步长,查一下loss曲线。: G# Y' @: U, w6 L: ^7 U7 Z

    8 s( ~5 L& m% ?* S' }6 J& d: I) ]或者把b但的起点改为1试试。
    回复 支持 反对

    使用道具 举报

  • TA的每日心情

    2025-9-22 22:19
  • 签到天数: 1183 天

    [LV.10]大乘

    5#
     楼主| 发表于 2023-2-15 00:25:26 | 只看该作者
    本帖最后由 雷达 于 2023-2-15 00:31 编辑
    - c, m* P8 G) t. e$ r+ F; {
    老福 发表于 2023-2-14 22:00; d. v8 [' q# ]  z; ~7 z3 Z
    刚才更新了一下,建议增加循环次数或调一下步长,查一下loss曲线。5 w5 O" _/ ]7 H% V2 V, V1 f
    - S% F) f  Q4 h9 d
    或者把b但的起点改为1试试。 ...
    4 y/ `; U( M2 f5 y; w8 U
    $ }2 \7 H: C. |( [
    你是对的。" B4 B& h2 `7 @# t) H1 w" F5 g" N) U3 j
    去掉了随机部分* O% Q  `, c& D/ o# ^- e
    #y = (x*27+15+random.randint(-2,3)).reshape(-1)2 _( }2 J  P* _# F/ X8 C) q
    y = (x*27+15).reshape(-1)
    ) s/ D3 j1 Y" }7 o+ F6 r+ O# Z* P# I5 {" \) g7 d
    循环次数加成10倍,就看到 b 收敛了
    ; p1 J% [+ w3 g9 `w , b7 y$ X; V5 f# m% h) p7 t
    27.002620697021484 14.826167106628418$ k7 J5 k- K" z' {! i4 g. P

    9 H+ b  k( M* l& Q, d3 |* @* q" o: U和 b 的起始位置无关,但 labeled data 用 y = (x*27+15+random.randint(-2,3)).reshape(-1) ,收敛就很慢。
    回复 支持 反对

    使用道具 举报

    手机版|小黑屋|Archiver|网站错误报告|爱吱声   

    GMT+8, 2026-3-7 05:28 , Processed in 0.052370 second(s), 17 queries , Gzip On.

    Powered by Discuz! X3.2

    © 2001-2013 Comsenz Inc.

    快速回复 返回顶部 返回列表