网页资讯视频图片知道文库贴吧地图采购
进入贴吧全吧搜索

 
 
 
日一二三四五六
       
       
       
       
       
       

签到排名:今日本吧第个签到,

本吧因你更精彩,明天继续来努力!

本吧签到人数:0

一键签到
成为超级会员,使用一键签到
一键签到
本月漏签0次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行补签。
连续签到:天  累计签到:天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
06月29日漏签0天
matlab吧 关注:291,752贴子:1,658,221
  • 看贴

  • 图片

  • 吧主推荐

  • 游戏

  • 2回复贴,共1页
<<返回matlab吧
>0< 加载中...

等号左侧的表达式不是用于赋值的有效目标。 有大神能帮忙看看吗

  • 只看楼主
  • 收藏

  • 回复
  • wl56198273
  • 自带板凳
    3
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
>> Untitled
错误: 文件:Untitled.m 行:6 列:7
等号左侧的表达式不是用于赋值的有效目标。


  • wl56198273
  • 自带板凳
    3
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
import numpy as np
import numpy.linalg as la
def soft(x,t):
y=np.abs(x)
y-=t
np.maximum(y,0,y)
y/=(y+t)
y*=x
return y
def fista(A,b,tau,L0=1e-4,eta=2,tol=1e-6):
H=np.dot(A.T,A)
c=np.dot(A.T,b)
L=L0
x=c
y=x
t=1
cost=la.norm(b-np.dot(A,x))**2/2+tau*la.norm(x,1)
while True:
x0=x
cost0=cost
fy=la.norm(b-np.dot(A,y))**2
fpy=np.dot(H,y)-c
while True:
x=soft(y-fpy/L,tau/L)
cost=la.norm(b-np.dot(A,x))**2/2+tau*la.norm(x,1)
if cost<=fy+np.dot(x-y,fpy)+L*la.norm(x-y)**2/2+tau*la.norm(x,1):
break
L*=eta
if abs((cost-cost0)/cost)<tol:
break
t0=t
t=(1+np.sqrt(1+4*t0**2))/2
y=x+(t0-1)/t*(x-x0)
return x


2025-06-29 04:00:11
广告
  • wl56198273
  • 自带板凳
    3
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
显示的是y-=t这一行错了


登录百度账号

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!
  • 贴吧页面意见反馈
  • 违规贴吧举报反馈通道
  • 贴吧违规信息处理公示
  • 2回复贴,共1页
<<返回matlab吧
分享到:
©2025 Baidu贴吧协议|隐私政策|吧主制度|意见反馈|网络谣言警示