yangyyyy_ 2023-03-03 20:31 采纳率: 70%
浏览 33
已结题

关于fortran矩阵运算的

现在就是有一个(6,60)(这应该是6列60行)的矩阵,然后一个一维行向量也是六个元素,让他两每一行对应元素相乘,我看了MATLAB里直接用.*乘就可以,不知道fortran里是不是也这样。同理一个(60,6)的矩阵和一个列向量,每一列对应相乘该怎么写程序呢

  • 写回答

3条回答 默认 最新

  • 地球屋里老师 2023-03-04 11:23
    关注
    
    program Test
    implicit none
    integer i
    !虽然fortran数组存储是列优先
    !但矩阵a是符合数学定义上的,2行3列
    !    1 3 5
    !    2 4 6
    integer:: a(2,3) =reshape([1,2,3,4,5,6],[2,3]), d(2,3)
    !一维数组既可以表示行向量,也可表示列向量
    !具体看你怎么用, 以下是矩阵乘法的例子
    integer:: b1(2)=1,b2(3)=1
    integer c1(3),c2(2)
    ! b1行向量,c1列向量
    c1 = matmul(b1,a)
    !结果 3,7,11
    print*,c1
    ! b2列向量,c2行向量
    c2 = matmul(a,b2)
    !结果 9,12
    print*,c2
    
    ! matlab中 .* 是对应元素相乘
    ! fortran中用 *
    ! 以下是示例
    do i=1, 3
      d(:,i) = a(:,i) * b1(:)
    end do
    do i=1, 2
      d(i,:) = a(i,:) * b2(:)
    end do
    end program
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论 编辑记录
查看更多回答(2条)

报告相同问题?

问题事件

  • 系统已结题 3月21日
  • 已采纳回答 3月13日
  • 创建了问题 3月3日

悬赏问题

  • ¥15 依据报错在原代吗格式的基础上解决问题
  • ¥15 在虚拟机中安装flash code
  • ¥15 单片机stm32f10x编写光敏电阻调节3.3伏大功率灯亮度(光强越大灯越暗,白天正常光强灯不亮,使用ADC,PWM等模块)望各位找一下错误或者提供一个可实现功能的代码
  • ¥20 verilog状态机方法流水灯
  • ¥15 pandas代码实现不了意图
  • ¥15 GD32H7 从存储器到外设SPI传输数据无法重复启用DMA
  • ¥25 LT码在高斯信道下的误码率仿真
  • ¥45 渲染完成之后将物体的材质贴图改变,自动化进行这个操作
  • ¥15 yolov5目标检测并显示目标出现的时间或视频帧
  • ¥15 电视版的优酷可以设置电影连续播放吗?