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